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The  Jet  Propulsion  Laboratory  has  developed  a computer  program^for  designing  and 
analyzing  the  performance  of  solar  array /battery  power  systems  for  the  U.S.  Coast 
Guard  Navigational  Aids.  This  program  is  called  the  Design  Synthesis/Perforaance 
Analysis  (DSPA)  Computer  Program.  The  basic  function  of  the  Design  Synthesis  portion 
of  the  DSPA  program  is  to  evaluate  functional  and  economic  criteria  to  provide 
specifications  for  viable  solar  array /battery  power  systems.  The  basic  function  of 
the  Performance  Analysis  portion  of  the  DSPA  program  is  to  simulate  the  operation 
of  solar  array/battery  power  systems  under  specific  loads  and  environmental  conditions. 

This  document  establishes  the  software  requirements  for  the  DSPA  computer  program, 
discusses  the  processing  that  occurs  within  the  program,  and  defines  the  necessary 
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INTRODUCTION 


The  Jet  Propulsion  Laboratory  ( JPL } has  developed  a computer  program  for 
designing  and  analyzing  the  performance  of  solar  array/battery  power 
systems  for  the  U.S.  Coast  Guard  Navigational  Aids.  This  program  is 
called  the  Design  Synthesis/Performance  Analysis  (DSPA)  Computer  Program. 
The  basic  function  of  the  Design  Synthesis  portion  of  the  DSPA  program  is 
to  evaluate  functional  and  economic  criteria  to  provide  specifications 
for  viable  solar  array /battery  power  systems.  The  basic  function  of  the 
Performance  Analysis  portion  of  the  DSPA  program  is  to  simulate  the 
operation  of  solar  array/battery  power  systems  under  specific  loads  and 
environmental  conditions. 

This  document  establishes  the  software  requirements  for  the  DSPA  computer 
program,  discusses  the  processing  that  occurs  within  the  program,  and 
defines  the  necessary  interfaces  for  operation. 
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PROGRAM  OVERVIEW 


The  DSPA  computer  program  combines  the  elements  of  design  synthesis  and 
performance  analysis.  A functional  block  diagram  of  the  main  driver  pro- 
gram is  shown  in  Figure  2-1.  As  shown  in  this  figure  and  the  algorithms 
which  follow,  the  DSPA  computer  program  utilizes  the  following 
methodology . 

a.  The  program  user  selects  the  desired  power  system  arrangement. 

b.  If  a design  synthesis  is  not  required,  then  the  program  user  must 
provide  information  on  the  electrical  size  of  the  equipment. 

c.  If  a design  synthesis  _is  requested,  the  program  user  must  supply 
information  on  the  parameters  used  in  determining  the  various  pro- 
files. The  computer  program  then  calculates  the  load  and  environ- 
mental profiles  needed  for  & local  profile  analysis.  Based  on  a 
profile  ener®r  balance  determined  as  part  of  the  load  profile  analy- 
sis, the  computer  estimates  the  electrical  size  of  the  equipment 
required  and  then  determines  the  physical  characteristics  of  the 
selected  equipment.  The  calculated  data  along  with  significant 
input  data  is  printed  out  in  the  appropriate  output  data  format. 

d.  If  performance  analysis  is  not  required,  the  execution  of  the  DSPA 
program  is  terminated. 

e.  If  performance  analysis  i_s  required,  the  program  user  must  provide 
information  on  the  parameters  used  in  determining  the  various  pro- 
files. The  computer  program  then  calculates  the  values  of  the  load 
and  environment  at  the  start  of  the  selected  mission  period.  These 
stimuli  (load  and  environment)  are  used  to  calculate  the  response 
(operational  characteristics)  of  the  equipment  at  that  point  in 
time.  The  process  is  repeated  for  selected  time  increments  until 
the  power  system  operational  characteristics  for  the  entire  mission 
period  have  been  determined.  This  information  is  then  printed  in 
the  appropriate  output  data  formats,  and  execution  of  the  program  is 
terminated. 


Flow  charts  of  the  DSPA  subprograms  were  not  furnished  in  the  Program 

Documentation  volumes  since: 

• Host  computer  facilities  have  programs  which  automatically  produce 
subroutine  flow  charts.  If  such  charts  are  desired,  the  program 
user  can  easily  select  the  subroutine  of  interest  and  obtain  a copy 
of  the  latest  version  of  the  subroutine. 

e Preparation,  reproduction,  and  inclusion  of  all  of  the  present  ver- 
sions of  the  DSPA  subroutines  in  the  Program  Documentation  would  be 
more  costly  than  if  the  flow  charts  were  prepared  by  the  program 
user  automatically.  Additionally,  these  flow  charts  would  become 
obsolete  as  modifications  were  made  to  the  DSPA  computer  program. 
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FIGUPE  2-1.  DESIGN  SYNTHESIS  PERFORMANCE  ANALYSIS  COMPUTER  PROGRAM 
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PROGRAM  ALGORITHMS 

Step  1 Obtain  Program  Parameters 
IPRG  = Program  Selector: 

0 = Design  Synthesis  only 

1 ■ Performance  Analysis  only 

2 = Both  Design  Synthesis  and  Performance  Analysis 
ITAPE  a Weather  Data  Input  Selector: 

-1  * Statistical  Input  Tape 
0 * User  Input  Data 

YYDDD  = Merge  Tape  Input  beginning  at  year  * YY  and 
day  * DDD 

DEBUG  * Debug  Printout  Request  Flag: 

0 ' So  printout 

1 * Printout 

XLN  » Length  of  X-axis  (in  inches)  for  summary  plots 
YLN  * Length  of  Y-axis  (in  inches)  for  summary  plots 

Step  2 Execute  Design  Synthesis  program  if  requested 

If:  IPRG  i 1 

Then:  Call  DSDRVR 

Step  3 Execute  Performance  Analysis  program  if  requested 

If:  IPRG  f 0 

Then:  Call  PADRVR 

Step  1*  STOP  DS/PA 
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3.  DESIGN  SYNTHESIS 


The  Design  Synthesis  portion  of  the  DSPA  program  uses  load  and  environ- 
mental profiles  to  set  the  power  system  requirements.  Based  on  these 
requirements  and  on  the  electrical  characteristics  of  the  system  equipment 
the  computer  program  determines  the  electrical  size  (volts,  amperes,  watts 
arape re-hours,  watt-hours)  and  the  physical  characteristics  (weight,  area, 
cost)  of  the  power  system.  A functional  block  diagram  of  the  Design 
Synthesis  driver  program  is  shown  in  Figure  3-1.  As  shown,  the  selection 
of  lamp  and  flasher  combinations  as  well  as  the  day/night  load  durations 
enables  the  computer  program  to  estimate  a power  load  profile.  This 
profile,  after  modification  using  battery  charge-efficiency  and  a number 
of  power  system  cabling  and  diode  losses,  is  used  in  the  load  profile 
analysis.  The  object  of  the  load  prof  le  analysis  is  to  determine 
the  electrical  size  of  a balanced*  power  source  as  well  as  the  minimum 
theoretical  electrical  size  of  the  battery.  Once  this  information  is 
obtained,  it  is  a fairly  straight  forward  process  to  determine  the 
electrical  size  of  the  remaining  items  of  equipment. 


Balanced  power  source:  a power  source  which  provides  Just  enough  energy 
to  the  batteries  (during  recharge  periods)  to  offset  or  balance  the  energy 
loss  sustained  during  discharge  periods. 
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PROGRAM  ALGORITHMS 

Step  1 Obtain  Pertinent  Mission  and  Equipment  Information 

Q0N  = Solar  Insolation  Level  for  Lamp  Flasher 
Turn-On  - Watts/Meter^ 

Q0FF  = Solar  Insolation  Level  for  Lamp  Flasher 
Turn-Off  - Watts/Meter^ 

VBUS  = Nominal  value  of  Raw  Power  Bus  Operating  level  - 
VDC 

VBU3MN  * Minimum  allowable  Raw  Power  Bus  Operating 
voltage  - VDC 

TTAVE  « Average  yearly  temperature  in  location  selected  - 
°F 

INDFLS  = Lamp  Flasher  Condition  Indicator 

Where:  0 * Lamp  Flasher  is  Off  (Lamp  Not  Flashing) 
1 ■ Lamp  Flasher  is  On  (Lamp  Flashing) 

DTTESG  * Energy  Storage  Group  Temperature  Rise  - °F 

QBRES  ■ Battery  Reserve  (as  Stage-of-Charge) 

BRCEST  ■ Estimated  Normalized  Battery  Charge  Current  - 
Hours” 1 

BRDEST  « Estimated  Normalized  Battery  Discharge 
Current  - Hours"! 

IOHRT  « Battery  Charger  Type 

Where:  0 ■ No  Battery  Charger 

1 ■ Constant  Voltage  Charger  with  Current 
Limit 

DTTPSG  ■ Power  Source  Group  Equipment  Temperature  Rise  - °F 

DURAM  ■ Duration  of  Mission  - Years 

SARES  ■ Solar  Array  Reserve  (as  a fraction  of  total  area) 

NPREQ  ■ Number  of  Solar  Cells  in  Parallel  Required  for 
each  Solar  Cell  Array  Electrical  Section 
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CELPAC  = Solar  Cell  Packing  Factor  on  Solar  Array 

IJSAP  = Total  Number  of  Solar  Arrays  to  be  Procured 

BRCHMX  = Maximum  Allowable  Normalized  Battery  Charge 
Current  - Hours”-*- 


= Shunt  Limiter  Type 


0  = No  Shunt  Limiter 


1 = Ordinary  Zener  Diode 

2 = Temperature-Compensated  Zener  Diode 

3 = Active  Shunt  Limiter 

Biasing  Factor  for  Selecting  the  Number  of  Storage 
Cells  in  Series  in  the  Battery 

(0.0  <_  FRCELL  <_  1.0) 

as:  FRCELL  -*0.0  the  minimum  number  of  cells 
tend  to  be  selected 

as:  FRCELL  ■*  1.0  the  maximum  number  of  cells 
tend  to  be  selected 


FRCELL  = 


BRDSTD  = Standard  Normalized  Battery  Discharge 
Current  - Hours”! 

TBDSTD  * Standard  Battery  Discharge  Temperature  - °F 


CBAVAL(JB) 


Table  of  Available  Storage  Cell  Capacities,  from  a 
Given  Manufacturer,  in  increasing  order  of 
size  - Amp-Hours/Cell 

(JB  * 1,30  maximum) 

Maximum  desired  capacity  of  each  battery  - 
Amp-Hours 

Maximum  Heat  Dissipation  of  a Single  Zener 
Diode  - Watts 

Heat  Dissipation  Derating  Factor  for  a Single 
Zener  Diode 


ACELL  ■ Single  Solar  Cell  Area  - cmc 


CBMAX 


HDZMX 
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Step  la  Compare  Raw  Power  Bus  Voltage  with  Allowable  Minimum 

I_f:  VBUS  < VBUSMN, 

Then:  VBUS  = VBUSMN  + 1*.0,  and. 

Then : Print  out  the  following  statement: 

"VBUS  ADJUSTED  TO  ( . . . VBUS  . . . ) VOLTS" 


Step  lb  Obtain  Yearly  Temperature  Extremes 

If:  ITAPE  = 0 

Then:  GO  TO  STEP  2 

If:  ITAPE  t 0 

Then:  Obtain  TTABMX  and  TTABMN  from  'MERGE'  File,  and. 
Then:  GO  TO  STEP  5 

Where:  TTABMX  = Maximum  Value  of  the  Ambient  Temperature  - 
°F 


TTABMN  * Minimum  Value  of  the  Ambient  Temperature  - 
°F 


Step  2 Calculate  Daily  Temperature  Increment  Extremes 
DTTA  ■ DTTA1  {DATE} 

DTTAMX  ■ Maximum  Value  of  DTTA  over  the  range: 

DATE  ■ 1;365  days 

DTTAMN  ■ Minimum  Value  of  DTTA  over  the  range: 

DATE  *il;365  days 

Where:  DTTVl  ■ Average  Daily  Temperature  Increment  - °f 

DATE  ■ Days  from  start  of  the  year  - days  (1-365) 
DTTA1  ■ Input  table  of  DTTA  as  a function  of  DATE 

Step  3 Calculate  Hourly  Temperature  Increment  Extremes 
DTTAMB  « DTAMB1  (TIMEH) 

DTABMX  ■ Maximum  Value  of  DDTAMB  over  the  range: 

• TIMEH  » 0;2fc  hours 
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Step  3 (contd) 

DTABMN  = Minimum  Value  of  DTTAMB  over  the  range: 

TIMEH  = 0;2l  hours 

Where:  DTTAMB  = Average  Hourly  Temperature  Increment  - °F 
TIMEH  = Daily  Time  - Hours  after  Midnight  (Q-2L) 
DTAMB1  = Input  Table  of  DTTAMB  as  a function  of  TIMEH 


Step  - Calculate  Ambient  Temperature  Extremes 
TTABMX  = TTAVE  + DTTAMX  + DTABMX 
TTABMN  = TTAVE  + DTTAMN  + DTABMK 


Step  5 Obtain  Power  Conditioning  and  Distribution  Group  (PCDG) 
Characteristics  at  the  Raw  Power  Bus 

(Based  on  ambient  temperature  extreme  which  will  yield  the 
largest  values  of  PCD  Group  Current) 

XX(J,K)  = Power  Conditioning  and  Distribution  Group  Voltage  - 
VDC 

XI ( J,K)  = Power  Conditioning  and  Distribution  Group  Current  at 
XX(J,K)  - Amperes 

1 * Lamp  Off 

J » 1,51  (Number  of  Data  Points) 

K = PCD  Group  Load  Selector  2 = Effective  Load, 

Lamp  Flashing 

3 = Lanp  On 


Step  6 Calculate  PCD  Group  Loads 

IF0FF  ■ F{XI ( J,l) ,XX(J,K)}  at  XX(J,K)  ■ VBUS 
IP0N  - F{XI(J,2),XX( J,K)}  at  XX(J,K)  ■ VBUS 
PP0N  - VBUS  • IF0N 
PF0FF  ■ VBUS  • IF0FF 
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Step  6 (contd) 

Where;  IF0FF  = PCD  Group  Lamp-Off  Operating  Current  - Amperes 

IF0N  = PCD  Group  Lamp-Flashing  Operating  Current  - 
Amperes 

PF0FF  = PCD  Group  Lamp-Off  Load  - Watts 
PF0N  * PCD  Group  Lamp-Flashing  Load  - Watts 

Step  7 Obtain  Free  Format  Data  on  Week  Number  and  Compare  with 
Reference 

LWEEK  = Weeks  after  Start  of  the  year  - (1,52) 

If:  LWEEK  < 0,  OR, 

If:  LWEEK  > 52, 

Then:  GO  TO  STEP  53 

NWEEK  = LWEEK 

Step  8 Calculate  Date  After  Start  of  the  Year 
DATE  ■ (7.0  • NWEEK)  - 6.0 

Step  9 Obtain  Terminator  Characteristics 

SRT  ■ Sunrise  Time  - Hours  after  Midnight 
SST  ■ Sunset  Time  - Hours  after  Midnight 

THETLA  ■ Buoy  Latitude  - Radians  J * south 

H0URT  ■ Terminator  Hour  Angle  - Radians 
ET  ..quation  of  Time  Difference  - Hours 
DECL  ■ Solar  Declination  Angle  - Radians 
ALPHA EQ  ■ Solar  Vector  Location  - Radians 
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Step  11  Calculate  Daily  Tine  Increment 
DTIMEH  = (SST  - SRT)/10.0 

Where:  DTIMEH  = Daily  Time  Increment  for  Clear  Day  Solar 
isolation  Calculations  - Hours 

Step  12  Initialize  Daily  Time  and  Time  Increment  Counter 
LTIME  = 1 
TIMEH  = SRT 

Where:  LTIME  = Time  Increment  Counter 

Step  13  Compare  Time  Increment  Counter  With  Reference 

If:  LTIME  = 1, 

IQDTC  =0.0) 

SALT  = O.oj  , AND,  GO  TO  STEP  l6 
QS0L(1)  = 0.0) 

Where:  QDTC  = Clear  Day  Solar  Insolation  Incident  on 

Solar  Array  - Watts/Meter^ 

SALT  = Solar  Altitude  - Radians 


Step  lL  Compare  Time  Increment  Counter  With  Reference 

LTIME  ■ 11, 

QDTC  = 0.0) 

SALT  « 0.0)  , AND  GO  TO  STEP  l6 
QSOL(ll)  ■ 0.0) 


Step  l^a  Compare  Environmental  Tape  Index  with  Reference 

.If:  ITAPE  ■ 0 

Then:  GO  TO  STEP  15 


Step  lU b Obtain  Solar  Insolation  from  "MERGE"  File 

QDT  ■ Solar  Insolation  Incident  on  Solar  Array  - 
Watts/Meter^ 
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Step  lUc  Calculate  Solar  Insolation  Incident  on  Solar  Array 
QSOL(LTIME)  = QDT 
GO  TO  STEP  16 

Where:  QSOL  * Solar  Insolation  Incident  on  Solar  Array  at 
TIMEC  - Watts /Meter2 


Step  15  Obtain  Clear  Day  Solar  Insolation  and  Solar  Altitude 

QDTC  “ Clear  Day  Solar  Insolation  Incident  on  Solar 
Array  - Watts/Meter^ 

SALT  * Solar  Altitude  - Radians 


Step  l6  Calculate  Clear  Day  Solar  Insolation  Array  and  Solar  Altitude 
Array 

QSOLC( LTIME)  « QDTC 
SALTA (LTINE)  « SALT 
TIMEC(LTIME)  ■ TIMER 

Where:  QSOLC  ■ Clear  Day  Solar  Insolation  Incident  on  Solar 
Array  at  TIMEC  - Watts/Meter^ 

SALTA  ■ Solar  Altitude  at  TIMEC  - Radians 

TIMEC  ■ Daily  Tine  - Hours  after  Midnight 

LTIME  ■ 1,11  (Number  of  data  points) 


Step  17  Increment  Daily  Time  and  Time  Increment  Counter 
LTIME  ■ LTIME  ♦ 1 


TIMER  ■ TIMER  ♦ DTIMEH 
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Step  18  Compare  Time  Increment  Counter  With  Reference 


If:  LTIME  > 11, 

Then:  GO  TO  STEP  19 


RETURN  TO  STEP  13 


Step  19 


Initialize  Day  Counter  and  Weekly  Summary  Arrays 
LDAY  = 1 
NMR(NWEEK)  = 0 
TJT(NWEEK)  ■ 0.0 
TKT(NWEEK)  =0.0 
TLT(NWEEK)  =0.0 
QQSOLT(NWEEK)  =0.0 
QSOLMX(NWEEK)  =0.0 

Where:  LDAY  = Day  of  the  Week  Indicator  - (Range  1,7) 


NMR  = Weekly  Number  of  Battery  Operating  Mode 
Reversals  - (Charging  or  Discharging) 


TJT  = Weekly  Duration  of  Solar  Occulatations  - Hours 


TKT  * Weekly  Duration  of  Share-Mode  Operations  - 
Hours 


TLT  » Weekly  Duration  of  Battery-Charging  Periods  - 
Hours 


Q3S0LT  * Weekly  Total  of  Solar  Insolation  Incident  on 
Solar  Array  - Watt-Hours /Mete r^ 


QSOLMX  ■ Maximum  Solar  Insolation,  Incident  on  Solar 
Array,  Encountered  during  week  - Watts/Meter" 


Step  19a  Compare  Environmental  Tape  Index  with  Reference 


If:  I7APE  ■ 0 

Then:  GO  TO  STEP  20 


GO  TO  STEP  23 
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Step  20  Obtain  Cloud  Cover  Conditions 
CT(LDAY)  = Cloud  type 

0.0  = Cirrus  or  Cirrostratus  Clouds 

1.0  - Stratus  Clouds 

2.0  - Other  Cloud  Types 
TC(LDAY)  = Total  Cloud  Cover 

1.0  * 1/10  of  sky  covered 

2.0  ■ 2/10  of  sky  covered 

9-0  = 9/10  of  sky  covered 

10.0  * 10/10  of  sky  covered 
ICT  - 1 + IFIX(CT) 

Where:  ICT  = Cloud  Type  indicator 

1 ■ Cirrus  or  Cirrostratus  Clouds 

2 ■ Stratus  Clouds 

3 * Other  Cloud  Types 

Step  21  Calculate  Cloud  Cover  Modifier 
If:  TC  ■ 0.0, 

Then:  CCfW(LTIME)  ■ 1.0  AND  GO  TO  STEP  22 

If:  SALTA(LTIME)  < w/4.0. 

Then:  I SALT  ■ 1 

If:  SALT(LTIME)  > ir/l».0. 

Then:  ISALT  » 2 

CCMM(LTIME)  » PC(ICT.ISALT)  ♦ Pl(lCT.ISALT)  • TC  ♦ ... 

♦ P2(ICT, ISALT)  • (TC*»2.0)  ♦ ... 

♦ P3(ICT, ISALT)  • (TC*»3.0) 

For  LTIME  • 1,11 
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' Step  21  (contd) 

Where:  ISALT  = Solar  Altitude  Indicator 

CCMM(LTIME)  = Cloud  Cover  Modifier 

P0,  PI,  P2,  P3  * Polynomial  Coefficients  obtained  from 

input  data  tables  "Cloud  Cover  Modifier 
Polynomial  Coefficients". 

Step  22  Calculate  Solar  Insolation  Incident  on  Solar  Array 
QSOL(LTIME)  = CCMM(LTIME)  * QSOLC(LTIME) 

For  LTIME  =1,11 


Step  2 j Calculate  Total  Daily  Solar  Radiation  Incident  on  Solar  Array 


J^LTIME«11 

f QS0L{ LTIME)  • DTIMEH 

LTIME=1 


(Using  Simpsons'  Rule 
( of  Integration 


Where: 


QQSOL  * Total  Daily  Solar  Insolation  Incident  on  Solar 
Array  - Watt-Hours /Meter? 


Step  2L  Calculate  Maximum  Solar  Radiation 

QSOLM(LDAY)  ■ AMAX  (QSOL(LTIME) ) over  the  range:  LTIME  « 1,11 

Where:  QSOLM  ■ Maximum  Solar  Insolation,  Incident  on  Solar 
Array,  Encountered  during  the  Day  - 
Watts /Meter* 

Step  25  Initialise  Time  Increment  and  Flasher  Load  Counters 

LTIME  » 2 
JTOFF  ■ 1 
JTON  ■ 1 

JOFFMX  ■ 0 
JONMX  ■ 0 
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JTOFF  = Flasher  Load  Turn-off  Period  Counter 
JTON  * Flasher  Load  Turn-on  Period  Counter 
JOFFMX  = Flasher  Turn-off  Counter 
JONMX  = Flasher  Turn-on  Counter 

Step  26  Compare  Solar  Insolation  Incident  on  Solar  Array  With  Reference 

If:  QSOL(LTIME)  > QOFF 

Then:  GO  TO  STEP  27 

GO  TO  STEP  31 

Step  27  Compare  Lamp  Flasher  Condition  Indicator  With  Reference 

if:  INDFLS  ■ 1 

Then:  GO  TO  STEP  28 

GO  TO  STEP  31 

Step  28  Calculate  Tine  of  Lamp  Flasher  Turn-Off 

TOFF (JTOFF)  - F (TIMEC(LTIME),  QSOL(LTIME)} 

at:  QSOL(LTIME)  ■ QOFF  and  during  the  tine 

interval:  TIMEC(LTIME  - 1)  to  TIMEC(LTIME) 

Where:  TOFF  ■ Tine  at  which  Lamp  Flasher  Turns  Off  - Hours 
after  Midnight 

Step  29  Reset  Lamp  Flasher  Condition  Indicator 
INDFLS  - 0 

Step  30  Increment  Load  Counters 

JTOFF  ■ JTOFF  ♦ 1 
JOFFMX  ■ JTOFF  - 1 
GO  TO  STEP  36 
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Step  31  Compare  Solar  Insolation  Incident  on  Solar  Array  With  Reference 

If:  QSOL(LTIME)  <_  QOII 

Then:  GO  TO  STEP  32 

GO  TO  STEP  36 

Step  32  Compare  Lamp  Flasher  Condition  Indicator  With  Reference 

I_f : INDFLS  = 0 

Then:  GO  TO  STEP  33 

GO  TO  STEP  36 

Step  33  Calculate  Tine  of  Lamp  Flasher  Turn-On 

TON(JTON)  = F{TIMEC(LTIME) , QSOL(LTIME)} 

at:  QSOL(LTIME)  = QOII,  and  during  the  time 

interval:  TIMES (LT I HE  - 1)  to  TIMES (LTIME) 

Where:  TON  « Tine  at  which  lamp  flasher  turns-on  - Hours  after 
Midnight 

Step  3I4  Reset  Lamp  Flasher  Conditions  Indicator 
INDFLS  = 1 

stcE-JS.  Increment  Load  Counters 

JTON  = JTON  ♦ 1 
JONMX  = JTON  - 1 

Step  36  Increment  Daily  Time  Counter 
LTIME  » LTIME  ♦ 1 

Step  37  Compare  Daily  Time  Counter  With  Reference 

If:  LTIME  >11, 

Then:  GO  TO  STEP  30 

RETURN  TO  STEP  26 


Step  38 


^eP  .39 


Step  1*0 


Step  1*1 


Step  1*2 


Ste£j£ 
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Calculate  Daily  Mode  Reversals  of  Battery 
DNMR  = JOFFMX  + JONMX 

Where:  DNMR  = Daily  Operational  Mode  Reversals  of  Battery 

Calculate  Daily  Solar  Occultation  Periods 
TJ  « SRT  + (21*. 0 - SST) 

Where:  TJ  * Duration  of  Daily  Solar  Occultation  Periods  - Hours 

Initialise  Daily  Share-Mode  and  Battery  Charging  Periods 
TL  « 0 

TK  * TOFF  (1)  - SRT 

Where:  TL  ■ Duration  of  Daily  Battery  Charging  Periods  - Hours 
TK  « Duration  of  Daily  Share-Mode  Periods  - Hours 

Initialise  Load  Counters 
JTQN  ■ 1 
JWFF  • 1 

Compare  Load  Counters 

If:  JTOFF  ■ JTON, 

Then:  GO  TO  STEP  1*3 

GO  TO  STEP  1*6 

Calculate  Daily  Battery  Charging  Periods 
DTL  ■ TON ( JTON ) - TOFF( JTOFF) 

TL  - TL  ♦ DTL 


Where:  DTL  ■ Battery  Charge  Period  Increment  - Hours 
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Ste£  1*1*  Increment  Turn-Off  Load  Counter 
JTOFF  = JTOFF  + 1 

Step  1*5  Compare  Turn-Off  Load  Counter  With  Reference 

I£:  JTOFF  > JOFFMX, 

Then:  00  TO  STEP  1*8 

RETURN  TO  STEP  1*2 

Step  1*6  Calculate  Daily  Share-Mode  Periods 
DTK  = T0FF( JTOFF)  - TOll(JTON) 

TK  = TK  ♦ DTK 

Where:  DTK  * Share-Mode  Period  Increment  - Hours 

Step  1*7  Increment  Turn-On  Load  Counter 
JTON  = JTON  1 
RETURN  TO  STEP  1*2 

Step  1*8  Calculate  Weekly  Summary  of  Operational  Mode  Periods 
TJT(IIWKEK)  ■ TJT(NWEEK)  ♦ TJ 
TKT(.TWEEK)  - TKT(HWEEK)  ♦ TK 
TLT(NVEEK)  ■ TLT(HWEEK)  ♦ TK 

Step  1*9  Calculate  Weekly  Summary  of  Battery  Mode  Reversals 
NMR(HWEEK)  ■ NMR(NVEEK)  ♦ DNMR 

Step  50  Calculate  Weekly  Summary  of  Solar  Insolation 
QQSOLT(NVEEK)  ■ QQSOLT(NVEEK)  ♦ QftSGL 
QSQLMX(NWEEK)  * AMAX  [qSOLMX(NWEEK) , QSOLM(LDAY)] 
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Step  51  Increment  Day  Counter 
LDAY  = LDAY  + 1 


Step  52  Compare  Day  Counter  With  Reference 

If:  LDAY  > 7, 

Then:  RETURN  TO  STEP  7 

RETURN  TO  STEP  20 


Step  53  Calculate  Yearly  Solar  Occultation  Load  Energy 


EPCDJ  = PFON  » 


NWEEK 

£ 


TJT(I) 


Where:  EPCDJ  = Yearly  Solar  Occultation  Load  Energy  - 
Watts-Hours 


Step  5L  Calculate  Yearly  Share-Mode  Load  Energy 


NWEEK 

EPCDK  * PFON  *. 

1*1 


TKT(I) 


Where:  EPCDK  * Yearly  Share-Mode  Load  Energy  - Watt-Hours 


Step  55  Calculate  Yearly  Battery  Charging  Mode  Load  Energy 


EPCDL  « PFOFK  * 


TLT(I) 


Where:  EPCDL  * Yearly  Battery  Charging  Mode  Load 
Energy  - Watt-Hours 
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Step  56  Calculate  Maximum  Energy  Storage  Group  Temperature 
TTESMX  = TTABMX  + DTTESG 

Where:  TTESMX  = Energy  Storage  Group  'iaximum  Temperature  - 
°F 

Step  37  Obtain  Normalized  Battery  Current  Rates 

BRR(J,K,L)  = BCQT  {VCC(J,K,L),  QBB(K),  TBB(L)} 

Where:  BRR  = Normalized  Battery  Current  Rates  (Expressed 

as  the  Ratio  of  Battery  Current  to  Battery 
Capacity)  - Hours“l 

VCC  = Cell  Voltage  - VDC 

QBB  = Battery  State-of-Charge 

TBB  = Battery  Temperature  - °F 

BCQT  = Input  Table  of  BRR  as  a function  of  VCC,  QBB 
and  TBB 

J * 1,  9 (Data  Points) 

K = 1,  NQBB  (Data  Points) 

L = 1,  NTBB  (Data  Points) 

NQBB  = Number  of  QBB  entries  in  BCQT 
NTBB  = Number  of  TBB  entries  in  BCQT 

Step  56  Calculate  Minimum  Storage  Cell  Discharge  Voltage 
VCDMN  = F {BRR(J,K,L),QBB(K),TBB(L)} 
at:  BRR(J,K,L)  ■ - |BRDEST| 

QBB(K)  ■ QBRES 
TBB(L)  ■ TTESMX 

Where:  VCDMN  * Minimum  Storage  Cell  Discharge  Voltage  - 
VDC 
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Step  59  Calculate  Minimum  Storage  Cell  Charge  Voltage 
VCCMN  = F {BRR(J,K,L),QBB(K),TBB(L)} 
at:  BRR( J,K,L)  * BRCEST 
QBB(K)  = QBRES 
TBB(L)  * TTESMX 

Where:  VCCMN  * Minimum  Storage  Cell  Charge  Voltage  - VDC 

Step  60  Calculate  Minimum  Battery  Potential  Efficiency 
ETABVN  = VCDMN /VCCMN 

Where:  ETABVN  ■ Minimum  Battery  Potential  Efficiency 

Step  6l  Calculate  Battery  State-of-Charge  Increment 
DQBB1  « (1.0  - QBRES)/10.0 

Where:  DQBB1  ■ Battery  State-of-Charge  Increment 

Step  62  Initialize  Battery  SOC  Counter  and  SOC  Values 
1QBB  » 1 
QBB1  - QBRES 

Where:  IQBB  ■ Battery  State-of-Charge  Counter 
QBB1  ■ Battery  3 1 at e-o f- Ch ar ge 
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Step 


Step  66 


Calculate  Instantaneous  Battery  Coulombic  (Charge) 
Efficiency 


ETA(IQBB)  = A (BRR,  QBB,  TBB} 
at:  BRR  = BRCEST 
TBB  = TTESMX 
QBB  = QBB1 


Where:  ETA(IQBB)  = Instantaneous  Battery  Coulombic 

Efficiency 


A = A series  of  Input  Data  Tables  (Al,  A2,  A3, 
Al* , A5)  giving  ETA  as  a function  of  BRR, 
QBB  and  TBB 


Step  61*  Increment  Battery  SOC  Counter  and  SOC  Values 
IQBB  = I QBB  + 1 
QBB1  = QBB1  + DQBB1 


Step  65  Compare  Battery  SOC  Counter  With  Reference 


If:  IQBB  > 11, 

Then:  GO  TO  STEP  67 


RETURN  TO  STEP  63 


Calculate  Average  Battery  Coulombic  (Charge)  Efficiency 

.IQBB-11 


1 ETA(IQBB)  • 00BB1 

IQBB-1  ( 


ETABQ 

Where:  ETABQ  ■ Average  Battery  Coulombic  Efficiency 


Using  Simpsons' 
Rule  of 
Integration 


Step  67  Calculate  Battery  Energy  Charge/ Discharge  Ratio 
RATBAT  ■ 1.0/(ETABVN  • ETABQ) 

Where:  RATBAT  » Battery  Energy  Charge/Discharge  Ratio 
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Step  68  Compare  Battery  Charger  Type  With  Reference 


If:  ICHRT  = 0, 

Then:  GO  TO  STEP  75 


Step  69  Calculate  Battery  Discharge  Line  Efficiency 
ETAD  = VBUS/(VBUS  + 1.0) 

Where:  ETAD  = Battery  Discharge  Line  Efficiency 


Step  70  Compare  Raw  Power  Bus  Operating  Level  With  Reference 

VCHIO  = VCHIOT  {TTESG}  at:  TTESG  = TTESMX 

If:  VBUS  < VCHIO, 

Then:  GO  TO_STEP  71* 

Where:  VCHIO  * Battery  Charger  Input  Voltage  at  Turn-on 

(Minimum  Voltage  Drop  at  zero  current  level)  - 
VDC 

VCHIOT  - Input  Table  of  VCHIO  as  a function  of  TTESG 
TTESG  ■ Energy  Storage  Group  Temperature  - °F 


[ 
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Step  71  Compare  Raw  Power  Bus  Operating  Level  With  Reference 

VCHISA  ■ VCHIST  {TTESG}  at:  TTESG  ■ TTESMX 

If:  VBUS  < VCHISA, 

Then:  GO  T0~STEP  72 

GO  TO  STEP  73 

Where:  VCHISA  ■ Battery  Charger  Input  Voltage  wherein  Charger 
Changes  From  "Saturated"  operation  to  "Active" 
Operation  - VDC 


VCHIST  ■ Input  Table  of  VCHISA  as  a function  of  TTESG 
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Steg  72  Calculate  Battery  Charger  Efficiency 
ETACHG  = 1.0  - (VCHIO/VBUS ) 

GO  TO  STEP  7 6 

Where:  ETACHG  = Battery  Charger  Efficiency 

Step  73  Calculate  Battery  Charger  Efficiency 
ETACHG  = ( VCHISA  - VCHIO)/VBUS 
GO  TO  STEP  76 

Step  '7u  Calculate  Battery  Charger  Efficiency 
ETACHG  ■ 0.0 
GO  TO  STEP  76 


Step  7?  Calculate  Battery  Discharge  Line  Efficiency  and  Charger 
Efficiency 

ETAD  « 1.0 

ETACHG  -1.0 


Step  76  Calculate  Duration  of  Yearly  Share-Mode  Loads 

EK 

TK7T  » 2-j 

>1 


Where:  TKTT  « Duration  of  Yearly  Share-Mode  Loads  - Hours 
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Step  77  Calculate  Duration  of  Yearly  Battery  Charge  Mode  Loads 


NWEEK 

TLTT  = 

1=1 


TLT(I) 


Where:  TLTT  = Duration  of  Yearly  Battery  Charge  Mode 
Loads  - Hours 


Step  78  Calculate  Power  Source  Group  (PSG)  Average  Power  Level 


PPSGAV  = 


RATBAT  * ( EPCDJ  ♦ EPCDK)  + (FT AD  » ETACHG  » EPCDL) 
(ETAD  « ETACHG  * TLTT ) * (RATBAT  « TKTT) 


Where:  PPSGAV  = PSG  Average  Power  Level  During  Year  - Watts 


Step  79  Calculate  PSG  Energy  Requirement 
EPSG  = PPSGAV  * (TKTT  ♦ TLTT) 

Where:  EPSG  = PSG  Yearly  Energy  Requirement  - Watt-Hours 

Step  80  Calculate  Solar  Array  Energy  Requirement 
ESA  ■ EPSG  j (TOifB;-STo) 

Where:  ESA  ■ Solar  Array  Yearly  Energy  Requirement  - 
Watt-Hours 

Step  8l  Calculate  Maximum  Solar  Array  Temperature 

TSAFMX  - TTABMX  ♦ DTTPSO 
TSARMX  - TSAFMX  ♦ 1*59.67 
TSAKMX  - (5.0/9. 0)  * TSARMX 
TSACMX  ■ TSAKMX  - 273.15 

Where:  TSAFMX  ■ Maximum  Solar  Array  Temperature  - °F 
TSARMX  ■ Maximum  Solar  Array  Temperature  - °R 
TSAKMX  ■ Maximum  Solar  Array  Temperature  - °K 
TSACMX  ■ Maximum  Solar  Array  Temperature  - °C 
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Step  82  Calculate  Maximum  Instantaneous  Solar  Radiation 

QDTMX  = AMAX  (QSOLMX(l))  over  the  range:  1=1,  NWEEK 

Where:  QDTMX  = Maximum  Instantaneous  Total  Solar  Radiation 
Incident  on  Solar  Array  - Watts/Meter2 

Step  83  Obtain  Solar  Array  Electrical  Characteristics 
Under  the  conditions: 

(a)  TSAC  = TSACMX 

(b)  QDT  = QDTMX 

(c)  DATEM  = DURAM  * 365.21*2 

(d)  NS1  = 50 

(e)  NP  = 50 

(f)  NESP  » 1 

Where:  TSAC  = Solar  Array  Temperature  - °C 

QDT  ■ Instantaneous  Total  Solar  Radiation  Incident 
on  Solar  Array  - Watts/Meter2 

DATEM  * Elapsed  Time  from  start  of  mission  - Days 

NS1  = No.  of  Solar  Cells  in  Series  in  each  circuit 

NP  0 No.  of  Solar  Cells  in  Parallel  in  each  circuit 

NESP  « Number  of  Parallel  Electrical  Circuits  in 
Solar  Array 

Then:  V2(L)  ■ Solar  Array  Voltage  - VDC 

I2{L)  ■ Solar  Array  Current  at  V2(L)  - Amperes 
L ■ 1,  MFINAL 

MSAPWR  ■ Solar  Array  Maximum  Power  - Watts 

MAXV  ■ Solar  Array  Voltage  at  Maximum  Power  Point  - 
VDC 

MAXI  ■ Solar  Array  Current  at  Maximum  Power  Point  - 
Amperes 
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Step  83a 


Step  81* 


Step  85 


Step  86 


Step  87 


Step  88 


Calculate  Single  Cell  Maximum  Power 
PWRMX1  = MSAPWR/(NS1  * NP  * NESP) 

Where:  PWRMX1  * Single  Cell  Maximum  Power  Output  - Watts 

Calculate  Solar  Cell  Efficiency  for  End  of  Mission 

-1* 

ETAEOM  = (PWRMXl/(QDT  * ACELL  * 10  ) 

Where:  ETAEOM  = End-of-Mission  Solar  Cell  Efficiency 


Calculate  Total  Yearly  Solar  Radiation 
NWEEK 

QQSTOT  = QQSOLT(I) 

in 

Where:  QQSTOT  •*  Total  Yearly  Solar  Radiation  Incident  on  Solar 
Array  - Watt-Hours/Meter2 


Calculate  Nominal  Estimate  of  Total  Solar  Cell  Area 

PCA 

ASCTNM  | 

QQSTOT  • ETAEOM  • (1.0  - SARES) 

Where:  ASCTNM  » Nominal  estimate  of  total  solar  cell  area  - 
Meter2 


Calculate  Nominal  Estimate  of  Total  Solar  Cells 

XNSCNM  ■ ASCTTW _ 

ACELL  • (1.0  * 10  ) 

Where:  XNSCNM  » Nominal  Estimate  of  Total  Solar  Cells 


Calculate  Solar  Array  Open  Circuit  Voltage 
VSC0C  » F (V2(L),I2(L))  at  I2(L)  - 0.0 
Where:  VSC0C  ■ Array  Open  Circuit  Voltage  - VDC 
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Step  8g  Calculate  Number  of  Solar  Cells  in  Series  Required 

ymq  - WUS  *1.0 
' / V3C0C  + MAXV  \ 

\ 2.0  * N31  f 

NS  = IFIX  ( XNS ) + 1.0 

Where:  XNS  = Estimate  of  Cells  in  Series 

NS  = Number  of  Solar  Cells  in  Series  Required 


Step  90  Calculate  Number  of  Solar  Array  Electrical  Sections  in  Parallel 


XNESP  = 


XNSCNM 

FLO AT (NS)  » FLOAT (NPREQ) 


NESP  = IFIX(XNESP)  + 1.0 


Where:  XNESP  = Estimate  of  Electrical  Sections  in  parallel 


Step  91  Calculate  Total  Solar  Cells  Required 
NSCT0T  * NESP  * NS  * NPREQ 

Where:  NSCT0T  * Total  Number  of  Solar  Cells  Required  for  the 
Solar  Array 


Step  92  Calculate  Total  Solar  Cell  Area 

ASCT0T  * NSCT0T  * ACELL  ( 1.0  * 10"S 

Where:  ASCT0T  * Total  Solar  Cell  Area  Required  for  Solar 
Cell  Array  - Meter*-’ 

Step  21  Calculate  Total  Solar  Array  Area 

ASA  * _ 

CELPAC  • (9.290  * 10“?) 

p 

Where:  ASA  * Total  Solar  Array  Area  Required  - Foot  ‘ 
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Step  91*  Calculate  Solar  Array  Weight 
DWDA  * DWDAT  * ASA 
WSA  = ASA  * DWDA 

p 

Where:  DWDA  = Solar  Array  Specific  Weight  - Lbs/Ft 

WSA  * Solar  Array  Weight  - Lbs 
DWDAT  * Input  Table  of  DWDA  as  a function  of  ASA 


Step  95  Calculate  Solar  Array  Cost 
ASATP  = ASA  * FLOAT  (NSAP) 

DCDA  = PCDAT  {ASATP} 

CSA  = ASA  * DCDA 

2 

Where:  ASATP  * Total  Area  of  Solar  Arrays  to  be  procured  - Ft 
CSA  * Solar  Array  Cost  - $ 

DCDA  ■ Solar  Array  Specific  Cost  - $/Ft^ 

DCDAT  s Input  Table  of  DCDA  as  a Function  of  ASATP 

Step  96  Calculate  Minimum  Solar  Array  Temperature 
TSAFMN  ■ TTABMN  ♦ DTTPSG 
TSARMN  TSAFMN  ♦ 1*59.67 
TSAKMN  ■ (5. 0/9.0)  • TSARMN 
TSACMN  » TSA not  - 273.1) 

Where:  TSAFMN  » Minimum  Solar  Array  Temperature  - °F 

TSARMN  ■ Minimum  Solar  Array  Temperature  - °R 

TSAKMN  ■ Minimum  Solar  Array  Temperature  - °K 

TSACMf  ■ Minimum  Solar  Array  Temperature  - °C 
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Step  97 


Obtain  Solar  Array  Current  Voltage  Characteristics  at  the 
Beginning  of  Life 

Under  the  conditions: 


(a)  TSAC  = TSACMN 

(b)  QDT  = QDTMX 

(c)  DATEM  = 0.0 

(d)  NP  = NPREQ 


(e) 

(f) 


Then: 


NS  « MS  ( „ 

I from  above 

NESP  = NESP  {Calculations 

V2(L)  = Solar  Array  Voltage  - VDC 
I2(L)  = Solar  Array  Current  at  V2(L)  - Amperes 
L « 1,  MFINAL 

MSAPWR  ■ Solar  Array  Maximum  Power  - Watts 

MAXV  * Solar  Array  Voltage  at  Maximum  Power  Point  - VDC 
MAXI  * Solar  Array  Current  at  Maximum  Power  Point  - Amperes 


Step  98  Calculate  Solar  Array  Open  Circuit  Voltage  (at  BOL) 

VSAOC  « F (V2(L) , I2(L) } at  I2(L)  * 0.0 

Where:  VSAOC  * Solar  Array  Open  Circuit  Voltage  - VDC 

Step  99  Calculate  Maximum  Blocking  Diode  Power  Loss 

KSABDP  ' (vBli'°  l.o)‘  MSAPW 

Where:  MSABDP  * Maximum  Power  Loss  in  all  Solar  Array  Blocking 
Diodes  - Watts 
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Step  100 


Step  101 


Step  102 


Step  103 


i 

i 


■» 


Calculate  Solar  Array  Electrical  Section  Blocking  Diode  Rating 
PESBD  = MS ABDP/ FLOAT  (NESP) 

Where:  PESBD  = Electrical  Section  Blocking  Diode  Power 
Rating  - Watts 


Initialize  Weekly  Counter  and  Battery  State 
NWEEK  =1.0 
BSTATE  =0.0 
EBTHMX  =0.0 
EBTHMN  = 0.0 

Where:  BSTATE  = Relative  Energy  State  of  the  Battery  - 
Watt-Hours 

EBTHMX  * Maximum  Battery  Energy  State  - Watt-Hours 
EBTHMN  * Minimum  Battery  Energy  State  - Watt-Hours 

Compare  Weekly  Counter  With  Reference 

If:  NWEEK  < 0 ; OR 

If:  NWEEK  > 52  ; OR 

If:  NWEEK  > LWEEK  : 

Then:  GO  TO  STEP  113 


Calculate  Average  Weekly  Solar  Insolation 
QQSOLA  ■ QQSTOT/ FLOAT  (NWEEK) 

Where:  QQSOLA  ■ Average  Weekly  Solar  Radiation  Incident  on  Solar 
Array  - Watt-Hcurs/Meter^ 
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Step  10^  Calculate  Modified  PSG  Power  Output 


PPSGL 


PPSGAV  » QQSOLT(  iiWKEK ) 
QQSOLA 


PPSGK 


PPSGL  » QOFF 
QSOLMX(NWEEK) 


Where:  PPSGL  = Weekly  Average  PSG  power  output  during  battery 
charging  periods  - Watts 

PPSGK  = Weekly  Average  PSG  power  output  during  Share-Mode 
Operational  periods  - Watts 


Step  10^)  Calculate  Battery  Load  Profile  Power  Levels 
PBJ  * -PFON/ETAD 
PGK  = ~(PF0N  - PPSGK )/ETAD 
PBL  « ( ETACHG/RATBAT ) * (PPSGL  - PFOFF) 

Where:  PBJ  * Battery  Discharge  Power  during  Solar  Occultation 
periods  - Watts 

PBK  * Battery  Discharge  Power  during  Share-Mode  Operation 
periods  - Watts 

PBL  * Battery  Charge  Power,  adjusted  to  rtf  fleet  Net 

gain  in  discharge  energy,  during  Battery  Charging 
periods  - Watts 


Step  106  Compare  Share-Mode  Battery  Discharge  Power  with  Reference 

I£:  PBK  > 0.0, 

Then:  PBK  * 0.0 


3-J0 


Step  107  Calculate  Battery  Energy  Profile  Levels 
EBX(l)  = PBJ  * TJL(NWEEK) 

EBX ( 2 ) = PBK  * TKT(NWEEK) 

EBX( 3)  = PBL  * TLT(NWEEK) 

EBDX(l)  = EBX(l)/7.0 
EBDX(2)  = EBX(2)/7.0 
EBDX( 3)  « EBX(3)/7.0 

Where:  EBX(I)  = Battery  Discharge  Energy  during  a particular 

operational  period  - Watt-Hours 

EBDX(I)  * Average  Daily  Battery  Discarge  Energy 

1=1=  During  Solar  Occultation  periods 

* 2 = During  Share-Mode  periods 

- 3 = During  Battery  charging  periods 

Step  108  Initialize  Period  Counter 

1=1 

Step  1C9  Compare  Period  Counter  With  Reference 

If:  I > 3, 

Then:  GO  TO  STEP  112 

Step  110  Calculate  Battery  Energy  State 

BSTATE  * BSTATE  + EBDX(I) 

EBTHMX  = AMAX  (EBTHMX,  BSTATE) 

EBTHMN  = AMIN  ( EBTHMN,  BSTATE) 


501*0-27 


M 


Step  111  Increment  Period  Counter 

1 = 1 + 1 

RETUB N TO  STEP  109 

Step  112  Increment  Weekly  Counter 

NWEEK  = NWEEK  + 1 
RETURN  TO  STEP  102 

Step  113  Calculate  Theoretical  Battery  Discharge  Energy  Requirement 
EBDTH  = ABS  ( EBTHMX  - EBTHMN ) 

Where:  EBDTH  = Theoretical  Battery  Discharge  Energy  -Watt-Hours 

Step  111*  Calculate  First  Criterion  for  Energy  Storage  Group  Discharge  Energy 
EESGl  = EBDTH/ (1.0  - QBRES) 

Where:  EESGl  = First  Criterion  for  ESG  Discharge 
Energy  - Watt-Hours 


Step  115  Calculate  Yearly  Battery  Mode  Reversals 
NWEEK 

NMRT  = NMR(I) 

1=1 


Where:  NMRT  = Total  Yearly  Reversals  of  Battery  Operating  Mode 

Step  life  Calculate  Total  Mission  Battery  Cycle  Requirements 

. FLOAT  [(^BffiAH)  . (Jjggjggy]  . 1 

Where:  NCYCLE  ■ Battery  Charge/Discharge  Cycle  Requirement 
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Step  117  Calculate  Theoretical  Battery  Depth-of-Discharge 
DOD  = DODT  {In  NCYCLE} 

Where:  DOD  = Theoretical  Battery  Depth  of  Discharge 

DODT  = Input  Table  of  DOD  as  a function  of  the  natural 
logarithm  of  NCYCLE 

Step  118  Compare  Theoretical  Depth  of  Discharge  With  Reference 

If:  DOD  < 0.0, 

Then:  DOD  * 0.001 

If:  DOD  > 1.0, 

Then:  DOD  * 1.0 

Step  119  Calculate  Second  Criterion  for  Battery  Discharge  Energy  Requirement 
EESG2  = EBDTH/DOD 

Where:  EESG2  * Second  Criterion  for  ESG  Discharge 
Energy  - Watt-Hours 

Step  120  Calculate  PSG  Maximum  Power  Output 
PPSGMX  « MSAPWR  - MSAPDP 

Where:  PPSGMX  ■ PSG  Maximum  Power  Output  - Watts 

Step  121  Calculate  Maximum  Battery  Charge  Power 

PBCHMX  ■ ETACHG  * (PPSGMX  - PF0FF) 

Where:  PBCHMX  ■ Maximum  Power  Available  for  Battery 
Charging  - Watts 
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Step  122 


Step  123 


Step  12U 


Step  125 


Step  126 


Calculate  ESG  Minimum  Temperature 
TTESMN  = TTABMN  + DTTESG 

Where:  TTESMN  - Minimum  Temperature  of  ESG  - °F 


Calculate  Maximum  Storage  Cell  Discharge  Voltage 
VCDMX  = F { BRR ( J , K ,L ) , QBB ( K ) , TBB ( L ) } 
at:  BRR(J,K,L)  = -|BRDEST| 

QBB(K)  = 1.0 
TBB(L)  = TTESMN 

Where:  VCDMX  = Maximum  Storage  Cell  Discharge  Voltage  - VDC 

Calculate  Maximum  Storage  Cell  Charge  Voltage 
VCCMX  = F {BRR( J,K,L),QBB{ K) ,TBB(L) } 
at:  BRR(J,K,L)  = BRCEST 
QBB(K)  = 1.0 
TBB(L)  * TTESMN 

Where:  VCCMX  * Maximum  Storage  Cell  Charge  Voltage  - VDC 

Calculate  .Maximum  Battery  Potential  Efficiency 
ETABVX  = VCDMX /VCCMX 

Where:  ETABVX  * Maximum  Battery  Potential  Efficiency 

Calculate  Third  Criterion  for  Battery  Discharge  Ener^  Requirement 

rrcr,  , PBCHMX  » ETABVX 
EESG3  dqouiiv 


Where:  EESG3  * Third  criterion  for  ESG  Discharge 
Energy  - Watt-Hours 
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Step  127  Calculate  Actual  Battery  Energy  Discharge  Capability 
EESG  = AMAX(EESG1,  EESG2,  EESG3) 

Where:  EESG  = ESG  Discharge  Energy  Capability  - Watt-Hours 


Compare  Battery  Charger  Type  With  Reference 

If:  ICHRT  > 0, 

Then:  GO  TO  STEP  133 


Step  130  Compare  Shunt  Limiter  Type  With  Reference 

I£:  ISH  > 0, 

Then:  GO  TO  STEP  132 


Calculate  Maximum  Battery  Charge  Voltage 


VBCHMX  - VSAOC  - 1.0 


GO  TO  STEP  139 

Where:  VBCHMX  = Maximum  Battery  Charge  Voltage  - VDC 


Step  132  Calculate  Maximum  Battery  Charge  Voltage 

VBCHMX  ■ VBUS  ♦ (0.75)  * {(VSAOC  - 1.0)  - VBUS) 
GO  TO  STEP  139 


Step  133  Compare  Shunt  Limiter  Type  with  Reference 

If:  ISH  > 0, 

Then:  ISH  • 0 


Step  13l*  Calculate  Battery  Charger  Reference  Voltages 


VCHIO  - VCHIOT  (TTESG) 
VCHISA  ■ VCHIST  (TTESG) 


| at: 


TTESG  ■ TTESMH 
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Step  13^  (contd) 


Where:  VCHIO  = Battery  Charger  Input  Voltage  at  Turn-On  - VDC 

VCHISA  = Battery  Charger  Input  Voltage  at  which  charger 
changes  from  "Saturated"  to  "Active" 
operation  - VDC 

VCHIOT  = Input  Table  of  VCHIO  as  a function  of  TTESG 
VCHIST  = Input  Table  of  VCHISA  as  a function  of  TTESG 


Step  135  Compare  Solar  Array  Open  Circuit  Voltage  With  Reference 


If:  ( VSAOC  - 1.0)  > VCHIO, 

Then:  GO  TO  STEP  136 

Print  Out  Error  Message: 

"VBUS  too  low  to  turn-on  battery  charger" 

Return  to  Main  Program  and  Stop  Computer  Run 


Compare  Solar  Array  Open  Circuit  Voltage  With  Reference 

If:  (VSAOC  - 1.0)  > VCHISA, 

Then:  GO  TO  STEP  138 


Step  137  Calculate  Maximum  Battery  Charge  Voltage 
VBCHMX  « (VSAOC  - 1.0)  - VCHIO 


GO  TO  STEP  139 


Step  138  Calculate  Maximum  Battery  Charge  Voltage 
VBCHMX  ■ F (VESA,  TTESG) 


at:  VESA  ■ VSAOC  - 1.0 


TTESG  ■ TTESMH 
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Step  138  (contd) 

Where:  VESA  = VCHIT  {VCHOOA,  TTESG} 

VCHOOA  = VBCHMX 

VESA  = Estimate  of  Battery  Charger  Input  Voltage  in 
"active"  condition  - VDC 

VCHOOA  = Battery  Charger  Output  Voltage,  in  "active" 
condition  at  zero  current  - VDC 

VCHIT  * Input  Table  of  VESA  as  a function  of  VCHOOA  and 
TTESG 

Step  139  Calculate  Storage  Cell  Minimum  Charge  Voltage 
VCCHMN  = F {BRR(J,K,L),QBB(K),TBB(L)) 
at:  BRR(J.K.L)  = 0.0 
QBB(K)  = 1.0 
TBB(L)  - TTESMX 

Where:  VCCHMN  * Storage  Cell  Minimum  Charge  Voltage  - VDC 

Step  lUO  Calculate  Maximum  Number  of  Storage  Cells  in  Series 
XNCBMX  - VBCHMX/ VCCHMN 

Where:  XNCBMX  * Estimated  Maximum  Number  of  Storage  Cells  in 
Series 

Step  lhl  Calculate  Storage  Cell  Maximum  Charge  Voltage 
VCCHMX  ■ F {BRR( J,K,L),QBB(K),TBB(L)} 
at:  BRR(J,K,L)  - 0.0 
QBB(K)  ■ 1.0 
TBB(L)  » TTESMN 

Where:  VCCHMX  ■ Storage  Cell  Maximum  Charge  Voltage  - VDC 


! 
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Step  1 UZ*  Calculate  Minimum  Number  of  Storage  Cells  in  Series 

XNCBMN  = VBUSMN/VCCHMX 

Where:  XNCBMN  = Estimated  Minimum  Number  of  Storage  Cells  in 
Series 

Step  lj*3  Calculate  Estimated  Storage  Cells  in  Series 

XNCELL  = XNCBMN  + (FRCELL)  * (XNCBMX  - XNCBMN) 

Where:  XNCELL  = Estimated  Storage  Cells  in  Series 

Step  ILL  Calculate  Actual  Number  of  Storage  Cells  in  Series 

HCELL  = 1F1X( XNCELL) 

YNCELL  = FLOAT (NCELL) 

DNCELL  * XNCELL  - YNCELL 

If:  DNCELL  >0.5, 

Then : NCELL  = NCELL  + 1 

Where:  NCELL  = Number  of  Storage  Cells  in  Series  in  the  Battery 

Step  ll*5  Initialize  Voltage  Counter  and  Battery  State-of-Charge 
I = 1 
QDS  = 0.0 

Where:  QBS  = Dummy  Variable  used  for  Battery  State-of-Charge 

Step  1L6  Compare  Voltage  Counter  With  Reference 

If:  I > 11, 

Then:  GO  TO  STEP  1**9 
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step  lVr 


Step  lU8 


Step  1^9 


Step  150 


Step  151 


Calculate  Storage  Cell  Discharge  Voltage 
VCDSTD(I)  = F {BRR(J,K,L),QBB(K),TBB(L)} 
at:  BRR(J,K,L)  = -|BRDSTD| 

QBB(K)  * QBS 

TBB(L)  = TBDSTD 

Where:  VCDSTD  = Storage  Cell  Discharge  Voltage  under  Standard 

Conditions  of  BRDSTD  and  TBDSTD  - Volts 

Increment  Voltage  Counter  and  Battery  State-of-Charge 
QBS  * QBS  + 0.1 
1=1  + 1 

RETURN  TO  STEP  lU6 


Calculate  Average  Storage  Cell  Discharge  Voltage 


J/'I-ll 

f VCDSTD(I)  • (0.1) 
I»1 

i 


Using  Simpsons’  Rule 


Where:  VCDAVG  ■ Average  Storage  Cell  Discharge  Voltage  Under 
Standard  Conditions  - VDC 


Calculate  Total  ESC  Discharge  Capacity 
CBDSTD  * 

VCDAVG  • FLOAT (N CELL) 

Where:  CBDSTD  * Total  ESC  Discharge  Capacity  Under  Standard 
Conditions  - Amp-Hours 


Initialize  Storage  Cell  Size  Counter 
JB  « 1 


Where:  JB  • Tabular  Location  Size  of  Available  Storage  Cell 
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Compare  Storage  Cell  Size  Counter  With  Reference 


1£:  JB  > 30 

Then:  GO  TO  STEP  151* 


Step  153  Compare  Storage  Cell  Size  With  Reference 


If:  CBAVAL(JB)  > 0.0, 

Then:  JB  = JB  + 1,  AND 
Then:  RETURN  TO  STEF  152 


Step  15l*  Calculate  Maximum  Number  of  Storage  Cell  Table  Entries 


JBTOT  = JB 


Where:  JBTOT  * Max.  No.  of  Entries  in  Storage  Cell  Table 


Initialize  Storage  Cell  Size  Counter 


JB  = 1 


Step  156  Compare  Storage  Cell  Size  Counter  With  Reference 

I£:  JB  > JBTOT, 

Then:  JBMAX  « JBTOT,  AND 
Then:  GO  TO  STEP  159  ” 

Where:  JBMAX  ■ Location  of  Max.  Available  Battery  Capacity  in 
Storage  Ceil  Table 


Step  157  Compare  Maximum  Desired  Battery  Capacity  With  Reference 

If:  CBMAX  > CBAVAL(JB), 

Then:  JB  * JB  ♦ 1,  AND 
Then:  RETURN  TO  STEP  156 


Step  1 


Calculate  Location  of  Maxinun  Desired  Battery  Capacity  in 
Available  Storage  Cell  Table 


JBMAX  ■ JB  - 1 
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Calculate  Estimate  of  Batteries  Required 


XNBATT 


CBDSTD 

CBAVAL(JBMAX) 


Where:  XNBATT  = Estimate  of  Number  of  Batteries  in  Parallel 


Step  l60  Compare  Estimate  of  Batteries  Required  With  Reference 

If:  XNBATT  >1.0, 

Then  JBB  = JBMAX,  AND 
Then:  GO  TO  STEP  l^T" 

Where:  JBB  * Location  of  Selected  Design  Capacity  in  Storage 
Cell  Table 


Step  161  Initialize  Counter 


J * 1 


Step  162  Compare  Counter  With  Reference  And  Select  Battery  Attribute 

I£:  J ■ JBMAX, 

Then:  NBATT  ■ 1 , AND. 

Then:  CBDT  ■ CBDSTD,  AND. 

Then:  CBD  ■ CBDT,  AND, 

Then:  GO  TO  STEP  ITT" 

Where:  NBATT  ■ Number  of  Batteries  in  Parallel 

CBDT  ■ Total  Discharge  Capacity  of  all 
batteries  - Amp-Hours 

CBD  * Discharge  Capacity  of  a single  battery  - Amp-Hours 


Step  163  Calculate  Location  of  Selected  Battery  Design  Capacity  in  Storage 
Cell  Table 

JBB  > JMAX  - J 


Step  164  Calculate  Estimate  of  Number  of  Batteries  in  Parallel 


XNBATT  * 


CBAVAL(JBB) 
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Step  165 

Compare  Estimate  of  Batteries  Required  With  Reference 

If:  XNBATT  < 1.0 

Then:  J = J + 1,  AND, 

Then:  RETURN  TO  STEP  162 

Step  166 

Compare  Estimate  of  Batteries  Required  With  Reference 

If:  XNBATT  <10.0 

Then:  GO  TO  STEP  172 

Step  167 

Increment  Maximum  Desired  Battery  Capacity  Location  in  Storage 

Cell  Table 

JBMAX  * JBMAX  + 1 

Step  168 

Compare  Maximum  Desired  Battery  Capacity  Location  With 
and  Select  Battery  Attributes 

If:  JBMAX  > JBTOT, 

Then:  NBATT  « 10,  AND, 

Then:  CBDT  « CBDSTD,  AND, 

Then:  CBD  * CBDT/FLOATTNBATT ) ; AND. 

Then:  GO  TO  STEP  171* 

Reference 

Step  169 

Calculate  Location  of  Selected  Battery  Design  Capacity 
Cell  Tabl^I.f- 

JBB  « JBMAX 

in  Storage 

Step  170 

Calculate  Estimate  of  Number  of  Batteries  in  Parallel 

XNBATT  ■ CBDSTD/ CBAVAL(JBB) 

Step  171 

Compare  Estimate  of  Batteries  Required  With  Reference 

If:  XNBATT  > 10.0, 

Then:  RETURN  TO  STEP  167 
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Step  172 


Step  173 


Calculate  Actual  Number  of  Batteries  in  Parallel 

NBATT  = IFIX(XNBATT) 

YNBATT  = FLOAT (NBATT) 

DNBATT  = XNBATT  - YNBATT 

If:  DNEATT  >.0.5, 

Then:  NBATT  = NBATT  + 1 


Calculate  Battery  and  ESC  Storage  Capacity 
CBD  = CBAVAL(JBB) 

CBDT  * CBD  * FLOAT( NBATT) 


Step  17*<  Calculate  Maximum  Allowable  Battery  Charging  Current 
EBDA  * CBDT  • VCDAVG  • FLOAT (NCELL) 

DO DA  » EBDTH/EBDA 
YICHMX  ■ CBDT  • BRCHMX 
ZICHMX  ■ YICHMX /FLOAT (NBATT) 

Where:  EBDA  ■ Total  Battery  Energy  Watt-Hours 

DO DA  ■ Maximum  Battery  Depth-of-Discharge 

YICHMX  ■ Maximum  Allowable  Battery  Charging  Current  for 
the  ESC  - Amperes 

ZICHMX  * Maximum  Allowable  Charging  Current  for  Single 
Battery  - Amperes 


Step  175  Calculate  Battery  Weight 


DWDE  = DWDET  { CBD } 

EBTOT  = CBDT  * VCDAVE  * FLOAT  (NCELL) 

WBATT  = DWDE  * EBTOT 

Where:  DWDE  = Battery  Specific  Weight  - Ibs/Watt-Hour 

EBTOT  = Total  Battery  Energy  in  ESG  - Watt-Hours 
WBATT  = Battery  Weight  - Lbs 

DWDET  * Input  Table  of  DWDE  as  a function  of  CBD 

Step  176  Calculate  Battery  Cost 

DCDE  * DCDET  {CBD,  NBATP} 

CBATT  * DCDE  * EBTOT 

Where:  DCDE  * Battery  Specific  Cost  - $/Watt-Hour 

NBATP  * Total  Number  of  Eatteries  to  be  procured 

CBATT  * BatVry  Cost  - $ 
f 

DCDET  ■ Input  table  of  DCDE  as  a function  of  CBD  and  NBATP 

f 

Step  177  Compare  Battery  Charger  Type  With  Reference 

If:  ICHRT  ■ 0, 

Then:  PCHO  ■ 0.0,  AND, 

Then:  WCHG  =*  0.0,  AND. 

Then:  CCHG  ■ 0.0,  AND. 

Then:  GO  TO  STEP  iir 

Where:  PCHO  ■ Maximum  Load  for  a Single  Battery  Charger  - Watts 
WCHG  ■ Weight  of  all  chargers  - Lbs 
CCHG  ■ Cost  of  all  chargers  - $ 
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Step  178  Calculate  Single  Charger  Maximum  Load 
PCHO  ~ P BCHMX / FLO AT ( N BATT ) 


Step  179  Calculate  Battery  Charger  Weight 
DWDCH  = DWDCHT  {PCHO} 

WCHG  = DWDCH  * PCHO  * FLOAT(NBATT) 

Where:  DWDCH  = Battery  Charger  Specific  Weight  - Lb/Watt 

DWDCHT  * Input  Table  of  DWDCH  as  a function  of  PCHO 

Step  l80  Calculate  Battery  Charger  Cost 
DCDCH  = DCDCHT  {PCHO,  NCHGP} 
at:  NCHGP  = NBATP 

Where:  NCAGP  « Number  of  Battery  Chargers  Procured 

DCDCH  * Battery  Charger  Specific  Cost  - $/Watt 

DCDCHT  * Input  Table  of  XDCH  as  a Function  of  PCHO 
and  NCHGP 

Step  l8l  Compare  Shunt  Limiter  Type  With  Reference 

If:  ISH  » 0, 

Then:  WSL  ■ 0.0,  AND. 

Then:  CSL  ■ 0.0,  AND. 

Then:  GO  TO  STEP  215 

Where:  WSL  ■ Total  Weight  of  Shunt  Limiters  - Lbs 
CSL  ■ Total  Cost  of  Shunt  Limiters  - $ 
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Step  182  Calculate  Shunt  Limiter  Operating  Point 
VSLOP  = VBCHMX 
XISLOP  = F { I2(L ) , V2( L ) } 

at:  V2(L)  = VSLOP 

Where:  VSLOP  = Shunt  Limiter  Operating  Point  Voltage  - VDC 

Xl'SLOP  = Shunt  Limiter  Operating  Point  Current  - Amperes 

Step  183  Compare  Shunt  Limiter  Type  With  Reference 

if.:  ISH  > 2, 

Then:  GO  TO  STEP  208 

Step  18^  Calculate  Power  Source  Group  Type 
IPSG  = 1 

Where:  IPSG  = Power  Source  Group  Type 

0 = One  Shunt  Limiter  for  the  Solar  Array 

1 = One  Shunt  Limiter  for  each  electrical  section  of 

the  solar  array 

Step  183  Calculate  Zener  Diode  Reference  Power  Level 
PZRF25  * HDER  * HDZMX 

Where:  PZRF25  * Single  Zener  Diode  Reference  Power  Level  at 
25°C  - Watts 

Step  166  Calculate  Estimate  of  Zener  Diodes  in  a Single  String 
XNZS  « PPSGMX/ ( PZRF25  * FLOAT(NESP)) 

Where  XNZS  ■ Estimate  of  Zener  Diodes  in  a Single  String 
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Step  187 


Step  188 


Calculate  Actual  No.  of  Zener  Diodes  in  a Single  String 

NZS  = IFIX(XNZS) 

YNZS  = FLOAT (NZS) 

DNZS  = XNZS  - YNZS 

If:  DNZS  1 0.5, 

Then:  NZS  * NZS  + 1 

If:  NZS  < 1, 

Then:  NZS  * 1 

Calculate  Zener  Diode  Operating  Temperature 
rcZ  = TSACMN 

Where:  TCZ  = Zener  Diode  Operating  Temperature  - °C 

Calculate  Single  Zener  Diode  Operating  Point 

VZOP  = vslop/float(nzs ) 

XIZOP  = XISLOP/FLOAT(NESP) 

Where:  VZOP  * Zener  Diode  Operating  Point  Voltage  - VDC 

XIZOP  * Zener  Diode  Operating  Point  Current  - Amperes 


Step  190  Compare  Shunt  Limiter  Type  With  Reference 

K:  ISH  > 1, 

Then:  GO  TO  STEP  201 


Step  191  Initialize  Zener  Diode  Voltage  Counter 
LZ  ■ 1 


Step  189 


Step  192  Initialize  Zener  Diode  Breakdown  Reference  Voltage 
VZB30  - VZOP 

Where:  VZB30  ■ Zener  Diode  Breakdown  Voltage  at  30°C  - VDC 


i 

I 
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Step  193  Compare  Zener  Diode  Voltage  Counter  With  Reference 

Ijf:  LZ  > 25 

Then:  VZBR  * VZB,  AND. 

Then:  GO  TO  STEP  205 

Where:  VZBR  = Single  Zener  Diode  Breakdown  Voltage  -VDC  (at 
Operating  Temperatures) 

VZB  = Estimate  of  Zener  Diode  Breakdown  Voltage  - VDC 


Step  19l*  Calculate  Zener  Diode  Temperature  Coefficient 
TCO  = ZTCOEF  {VZB 30) 

Where:  TCO  = Zener  Diode  Temperature  Coefficient  - (%/°C) 

ZTCOEF  » Input  Table  of  TCO  as  a function  of  VZB30 

Step  195  Calculate  Zener  Diode  Dynamic  Impedance 
ZZ  « ZDIMP  {TCZ,  VZB30} 

Where:  ZZ  * Zener  Diode  Impedance  - Ohms 

ZDIMP  * Input  Table  of  ZZ  as  function  of  TCZ  and  VZB30 

Step  196  Calculate  Estimate  of  Zener  Diode  Breakdown  Voltage 
VZB  « VZOP  - (XIZOP  • ZZ) 


Step  197  Calculate  Estimate  of  Reference  Temperature  Zener  Breakdown 
Voltage 

VZB 301  - VZB  • £l.O  - — *1^q  * 3°.‘ QJ j 

Where:  VZB 301  ■ Estimate  of  Zener  Diode  Breakdown  Voltage 
at  30°C  - VDC 
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Calculate  Breakdown  Voltage  Residual 


Step  199 


Step  200 


Step  201 


Step  20? 


DVZB . abs  (Brngjaa ) 

Where:  DVZB  = Residual  in  Estimate  of  Zener  Diode  Breakdown 
Voltage  at  30°C  - VDC 


Compare  Breakdown  Voltage  Residual  With  Reference 

If:  DVZB  _<  0.01, 

Then:  VZBR  = VZB,  AND. 

Then:  GO  TO  STEP  205 


Increment  Zener  Diode  Voltage  Counter  and  Adjust  Reference 
Breakdown  Voltage 

LZ  * LZ  ♦ 1 

VZB30  = VZB301 

RETURN  TO  STEP  193 


Calculate  Temperature  Compensated  (TC)  Zener  Reference  Current 
1ZRF25  » CURZ  (HDZMX) 

Where:  IZRF25  * Zener  Diode  Reference  Current  at  25°C  - Amperes 
CURZ  ■ Input  Table  of  IZRF25  as  a function  of  HDZMX 


Calculate  TC  Zener  Reference  Voltage 


VZRF25  ■ PZRF25/IZRF25 

Where:  VZRF25  ■ TC  Zener  Reference  Voltage  at  25°C  - VDC 
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Step  203 


Step  2Ql4 


Step  205 


Step  206 


Calculate  TC  Zener  Breakdown  Voltage  Ratio 
RATVB  = F {TCZ,  RATI} 

at:  RATI  = 0.0 
RATI  = TCZIV  {RATV,  TCZ} 

Where:  RATV  = Zener  Diode  Voltage  Ratio 

RATI  = Zener  Diode  Current  Ratio 
TCZIV  = Input  Table  of  RATI  as  a function  of  RATV  and  TCZ 
RATVB  = Zener  Breakdown  Voltage  Ratio 

Calculate  TC  Zener  Breakdown  Voltage 
VZBR  = RATVB  * VRF25 

Calculate  Zener  Diode  Temperature  at  Breakdown 
TZBR  = TCZ 

Where:  TZBR  = Zener  Diode  Temperature  at  Breakdown  Voltage  - °C 

Calculate  Total  Weight  of  Shunt  Limiter 
DWDNZ  = DWDNZT  {PRF25,  ISH} 

UZTOT  = NZS  * NKSP 
WSL  = DWDNZ  * FLOAT (HZT0T) 

Where:  DWDNZ  * Zener  Diode  Specific  Weight  - Lbs/Zener 

NZTOT  * Total  Number  of  Zener  Diodes 

DWDNZT  * Input  Table  of  DWDNZ  as  a function  of  PRF25 
and  ISH 
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Step  207  Calculate  Total  Cost  of  Shunt  Limiter 

NZTP  = NZTOT  » NSAP 

DCDNZ  = DCDNZT  {PRF25 , NZTP,  ISH} 

CSL  = DCDNZ  • FLOAT (NZTOT) 

Where:  NZTP  = Total  No.  of  Zener  Diodes  Procured 

DCDNZ  = Zener  Diode  Specific  Cost  - $/Zener 

DCDNZT  = Input  Table  of  DCDNZ  as  a function  of  PRF25, 
NZTP  and  ISH 

GO  TO  STEP  215 

Step  208  Calculate  Power  Source  Group  Type 
IPSG  = 0 


Step  209  Calculate  Shunt  Limiter  Reference  Temperature 
TSHREF  « TSACMN 

Where:  TSHREF  « Shunt  Limiter  Reference  Temperature 

Step  210  Calculate  Active  Shunt  Limiter  Impedance 
ZSH  « ZSHTAB  {TSH} 

at:  TSH  ■ Tf  HREF 

Where:  ZSH  ■ Active  Shunt  Limiter  Dynamic  Impedance  - Ohms 

ZSHTAB  ■ Input  Table  of  ZSH  as  a function  of  TSH 

Step  211  Calculate  Active  Shunt  Limiter  Turn-On  Voltage 
VSHTOR  « VSLOP  - XISLOP  • ZSH 

Where:  VSHTOR  ■ Shunt  Limiter  Tum-on  Voltage  - VDC  at  TSHREF 
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Step  212  Calculate  Active  Shunt  Limiter  Load 
PSL  = PPSGMX 

Where:  PSL  = Active  Shunt  Limiter  Load  - Watts 

Step  213  Calculate  Shunt  Limiter  Weight 
DWDPS  = DWDPST  {PSL} 

WSL  = DWDPS  * PSL 

Where:  DWDPS  = Shunt  Limiter  Specific  Weight  - Lbs/Watt 

Step  21^4  Calculate  Total  Cost  of  Shunt  Limiter 
NSLP  = NSAP 

DCDPS  = DCDPST  (PSL, NSLP} 

CSL  - DEDPS  * PSL 

Where:  NSLP  * No.  of  Active  Shunt  Limiters  Procured 

DCDPS  * Active  Shunt  Limiter  Specific  Cost  - $/Watt 

DCDPST  * Input  Table  of  DCDPS  as  a function  of  PSL 
and  NSLP 

Step  215  Calculate  Total  Power  System  Weight 
WPWR  » WSA  ♦ WBATT  + WCHG  ♦ WSL 
Where:  WPWR  = Total  Power  System  Weight  - lbs 

Step  Pl6  Calculate  Total  Power  System  Cost 
CPWR  ■ CSA  ♦ CBATT  ♦ CCHG  ♦ CSL 
Where:  CPWR  » Total  Power  System  Cost  - $ 
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Step  217 

Step  218 


Print  Design  Synthesis  Output  Information 


Set  Performance  Analysis  Battery  Parameters 
CB  = CBD 

XN  = FLOAT (NCELL) 

XICHMX  = ZICHMX 


RETURN  TO  MAIN  PROGRAM 
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3. 1 Terminator  Characteristics 


The  Terminator  Characteristics  routine  is  used  to  calculate  the  sunrise 
and  sunset  times  for  a specific  day  of  the  year.  The  equations  use  solar 
vector  location  angles,  time  zone  number,  and  buoy  latitude  and  longitude 
to  determine  terminator  hour  angles  and  times. 


PROGRAM  ALGORITHMS 


Step  1 

Obtain 

Exact  Date 

DATE  = 

Date;  Days  from  the  start  of  the  year  (1,365) 

Step  2 

Calculate  Solar  Vector  Location  in  Equatorial  Plane 

ALPHEQ 

= OMEGA  » DATE 

Where: 

ALPHEQ  * Solar  Vector  Location  - Radians 

OMEGA  = (2  • it)/ 365. 2*42 

TT  ■ 3.1*4159 

Note: 

There  are  365.2*42  days  per  tropical  year  as  measured 
from  Vernal  Equinox  to  Vernal  Equinox 

Step  3 

Calculate  Solar  Radiation  Variables 

VAR(I) 

« FAO(I)  ♦ FA1(I)  • COS (ALPHEQ)  ♦ ... 

♦ FA2(I)  • C0S(2.0  • ALPHEQ)  ♦ ... 

♦ FA3(I)  • COS ( 3.0  • ALPHEQ)  ♦ ... 

♦ FB1(I)  • SIN(ALPHEQ)  ♦ ... 

♦ FB2(I ) • SZN(2.0  • ALPHEQ)  ♦ ... 

♦ FB3(I)  • SZN( 3.0  • ALPHEQ) 

DECL 

■ VAR(l)  • rt/180.0 

ET 

■ VAR(2) 

Where: 

DECL  a Solar  Declination  Angle  - Radians 

ET  a Equation  of  Tine  Difference  - Hours 

FA,  FB  ■ Fourier  Coefficients  obtained  from  input 
data  tables  "Solar  Radiation  Fourier 

Coefficients" 
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Step  3a  Calculate  Solar  Radiation  Variables 
IF:  ITAPE  * 0, 

THEN:  APPSC  = VAR(3)  * 3.152UQ0Q,  and. 

THEN:  ATMEXC  = VAR(U),  and. 

THEN:  SDF  = VAR(5) 

p 

Where:  APPSC  = Apparent  Solar  Constant  - Watts/Meter 
(at  AMO) 

ATMEXC  = Atmospheric  Extinction  Coefficient  - Air  Mass-^ 
SDF  = Sky  Diffuse  Factor 

Step  U Obtain  Buoy  Latitude 

THELAD  * Buoy  latitude  - degrees  South} 

Step  5 Convert  Buoy  Latitude 

THETLA  - THELAD  • tt/180.0 

Where:  THETLA  ■ Buoy  Latitude  - Radians 

Step  6 Calculate  Terminator  Hour  Angle 

IF:  [THETLA  _>  (n/2.0)  - DECL], 

THEN:  HOURT  ■ 7,  AND.  Go  to  Step  7 

HOURT  ■ ARCCOS  (-1.0  • TAN  (THETLA)  • TAN  (DECL)) 

Where:  HOURT  ■ Terminator  Hour  Angle  - Radians 

StegJ.  Convert  Terminator  Hour  Angle 

HOURA  « HOURT  • 12.0/w 

Where:  HOURA  ■ Terminator  Hour  Angle  - Hours 

Step  6 Obtain  Buoy  Location  Time  Zone  Number 

T£N  ■ Time  Zone  Number  (Hours  behind  Greenwich  Mean  Time) 
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Step  9 Obtain  Buoy  Longitude 

THELOD  = Buoy  Longitude  - degrees 


+ West 
East 


Step  10  Calculate  Time  of  Sunrise  and  Sunset  at  Buoy  Location 
SRT  = 12.0  - HOURA  - ET  - TZN  + (THEL0D/15.0 ) 

SST  = 2U.0  - SRT 

Where:  SRT  = Sunrise  Time  - Hours 
SST  * Sunset  Time  - Hours 
RETURN  TO  DESIGN  SYNTHESIS  DRIVER  ROUTINE 
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3.2  Clear  Day  Solar  Insolation 

The  Clear  Day  Solar  Insolation  routine  is  used  to  compute  the  intensity 
of  solar  radiation  incident  on  the  buoy  solar  array  for  a particular  day 
and  time.  The  insolation  is  calculated  for  a clear  day  (i.e.,  no  cloud 
cover)  using  solar  array  tilt  angles,  buoy  location  hour  angles,  sky 
diffuse  factor,  clearness  numbers,  and  surface  reflectivity. 

PROGRAM  ALGORITHMS 

Step  1 Obtain  Geometrical  and  Temporal  Information 

TIMEH  = Daily  Time  - Hours  after  Midnight  - (0,21*) 

TZN  = Buoy  Location  Time  Zone  Number  - (Hours  behind 
Grennwich  Mean  Time) 

THELOD  * Buoy  Longitude  - Degrees  |+ 

HOURT  * Terminator  Hour  Angle  - Radians 
DECL  * Solar  Declination  Angle  - Radians 

THETLA  * Buoy  Latitude  - Radians 

Step  2 Calculate  Buoy  Location  Hour  Angle 

BHOURD  ■ 15.0  • (TIMEH  - 12.0  ♦ TZN  ♦ ET)  - THELOD 
BHOUR  ■ BHOURD  • */l80.0 

Where:  BHOURD  ■ Buoy  Location  Hour  Angle  - Degrees 
BHOUR  ■ Buoy  Location  Hour  Angle  - Radians 

Step  3 Compare  Buoy  Location  Hour  Angle  With  Terminator  Hour  Angle 
(Test  for  Solar  Occulation) 

IF:  ASB( BHOUR)  > ABS( HOURT) 

THEN:  GO  TO  STEP  23 
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Step  U Calculate  Direction  Cosines  of  Direct  Solar  Radiation 

COS(THETZS)  = COS ( BHOUR ) * COS(DECL)  * COS(THETLA)  + 

....+  SIN(DECL)  * SIN(THETLA) 

COS(THW)  = COS(DECL)  * SIN (BHOUR) 

IF:  COS (BHOUR)  >{^jSA)). 

THEN;  KS  “ 1.0 

IF:  COS  (BHOUR)  < {^S*)  } ■ 

THEN:  KS  ■ -1.0 

/ p pi  0. 5 

COS(THS)  ■ KS*  Ul-COS(THETZS)]  - [C0S(THV.-}]^ 

Where:  THETZS  = Angle  between  the  local  zenith  and  the  solar 
Vector  - Radians 


THW,  THS  = Additional  Direction  Angles  - Radians 


Step  5 Calculate  Solar  Altitude 
S ALT  « AP  CS I N ( COS { THETZS ) ) 

Where:  SALT  ■ Solar  Altitude  (Angle  between  the  solar  vector  and 
the  Horizontal,  i.e.,  Earth's  surface)  - Radians 


Step  6 Calculate  Solar  Azimuth 
IF:  COS (THS)  > 0, 

THEN:  SAZM  » ARCSIN(COS(THW)/COS(SALT) ) , AND  GO  TO  STEP  > 
IF:  COS (THS)  < 0, 

THEN:  SAZM  ■ n - ARCSIN(COS(THW)/COS(SALT) ) 

Where:  SAZM  ■ Solar  Azimuth  (Angle  between  the  Solar  Vector 
Projected  onto  the  Horizontal  Surface  and  the 
South-Pointing  Vector  on  the  Horizontal 
Surface)  • Radians 


Step  7 Obtain  Clearness  Number 
CN  ■ Clearness  Number 

■ 0. 7-9.9  for  an  industrial  atmosphere 

■ 0.85-1.10  for  non-industrial  atmospheres 
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Step  8 Obtain  Solar  Radiation  Variables 

r 

APPSC  = Apparent  Solar  Constant  at  AMO  - Watts /Meter^ 
ATMKXC  = Atmosphere  Extinction  Coefficient  - Air  Mass'1 
SDF  = Sky  Diffuse  Factor 


Step  9 Calculate  Intensity  of  Direct  Normal  Solar  Radiation 
QDN  = APPSC  * CN  « EXP(-ATMEXC/COS(THETZS) ) 

Where:  QDN  = Direct  Normal  Solar  Radiation  Intensity  - Watts/ 
Me  ter ^ 


Step  30  Obtain  Solar  Array  Pointing  Angles 

PHIAID  = Surface  Tilt  Angle  from  Horizontal  - Degrees 

(Angle  between  local  Zenith  and  Solar  Array  Normal) 

PHIAAD  = Surface  Azimuth  Angle  from  South  - Degrees 

(Angle  between  South  pointing  vector  and  projection  of 
array  normal  on  horizontal  surface) 


(+  if 
\-  if 


West  of  South 
East  of  South 


Step  11  Convert  Solar  Array  Pointing  Angles 

PHIAI  * PHIAID  * tt/180.0 
PHIAA  = PHIAAD  * tr/l80.C 

Wher  . PHIAI  ■ Surface  Tilt  Angle  ••  Radians 

PHIAA  * Surface  Azimuth  Angle  - Radians 


Step  12  Calculate  Direction  Cosines  of  Array  Normal 

(Reference  Axis:  Vertical,  Horizontal  to  West,  Horizontal  to 
South) 

ETAA  ■ COS (PHIAI) 

ETAB  - SIN (PHIAA)  * SIN(PHIAl) 

ETAC  ■ COS (PHIAA)  * SIN (PHIAI) 

Where:  ETAA,  ETAB,  ETAC  are  Array  Normal  Direction  Cosines 
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Step  13  Calculate  Solar  Array  Tilt  Angle 

COS(TILT)  = ETAA  * COS(THETZS)  + ... 

+ ETAB  * COS(THW)  + ... 

+ ETAC  * COS(THS) 

Where:  TILT  = Solar  Array  Tilt  Angle  - Radians 

(Angle  between  Solar  Vector  and  Solar  Array 
Normal ) 

Step  lL  Calculate  Intensity  of  Direct  Solar  Radiation  Incident  on  the 
Solar  Array 

IF:  COS  (TILT)  >0.0, 

THEN:  QD  = QDN  * COS(TILT) 

IF:  COS  (TILT)  <_  0.0, 

THEN:  QD  = 0.0 

Where:  QD  = Direct  Solar  Radiation  Incident  on  Solar  Array  - 
Watts/Meter2 

Step  15  Calculate  Sky  Brightness 

BS  = SDF  * QDN/ (CN  **  2.0) 

Where:  BS  = Sky  Brightness  - Watts/Meter2 

Step  16  Obtain  Horizontal  Surface  (Ground/Ocean)  Reflectivity 

REFLH  * Horizontal  Surface  Reflectivity  for  Solar  Radiation 

Step  17  Calculate  Horizontal  Surface  Brightness 

BG  ■ REFLH  * (BS  ♦ QDN  • COS(THETZS)) 

p 

Where:  BG  ■ Horizontal  Surface  Brightness  - Watts/Meter 

Step  18  Calculate  Intensity  of  Horizontal  Surface  Diffuse  Radiation 
Incident  on  Solar  Array 

QDG  ■ BG  • ((1  - ETAA)/2.0) 

Where:  QDG  ■ Horizontal  Surface  Diffuse  Radiation  Incident  on 
Solar  Array  - Watts /Meter^ 
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Calculate  Intensity  of  Sky  Diffuse  Radiation  Incident  on  a 
Horizontal  Solar  Array 

QDSH  = QDN  * SDF 

Where:  QDSH  = Sky  Diffuse  Radiation  Incident  on  a Horizontal 
Solar  Array  - Watts/Meter2 


Step  20  Calculate  Intensity  of  Sky  Diffuse  Radiation  Incident  on  a 
Vertical  Solar  Array 

YV  = 0.1*5 

IF:  COS (TILT)  > (-0.20), 

THEN:  YV  = 0.55  + 0.1*37  * COS(TILT)  + 0.313  * ( COS ( TILT ) **2. 0 ) 

QDSV  = QDN  * (SDF  * YV  + ( REFLH  * (SDF  + COS(THETZS) ) )/2.0) 

Where:  QDSV  = Sky  Diffuse  Radiation  Incident  on  a Vertical 
Solar  Array  - Watts/Meter2 

Step  21  Calculate  Intensity  of  Sky  Diffuse  Radiation  Incident  on  Solar 
Array 

QDS  = QDSV  + (QDSH  - QDSV)  * COS(SALT) 

Where:  QDS  * Sky  Diffuse  Radiation  Incident  on  Solar  Array 
- Watts/Meter2 


Step  22  Calculate  Intensity  of  Total  Clear  Day  Solar  Insolation  Incident 
on  Solar  Array 

QDTC  = QD  + QDG  + QDS 

Where:  QDTC  ■ Total  Clear  Day  Solar  Radiation  Incident  on 
Solar  Array  - Watts/Meter2 

Return  to  Design  Synthesis  Driver  Routine 


Step  23  Calculate  Occultation  Conditions  for  Solar  Insolation 

QDN  ■ 0.0 
QD  ■ 0.0 
QDG  ■ 0.0 
QDS  ■ 0.0 
QDTC  ■ 0.0 


Return  to  Design  Synthesis  Driver  Routine 
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3. 3 Solar  Array  Electrical  Characteristics 

The  Solar  Array  Electrical  Section  is  made  up  of  the  solar  array,  the 
solar  array  isolation  diodes,  and  the  power  source  series  resistance. 

The  characteristics  of  these  elements  are  calculated  for  the  environ- 
mental conditions  in  which  the  subsystem  will  operate  and  are  then 
combined  into  a single  solar  array  current-v oltage  curve.  Performance 
data  are  stored  for  a single  solar  cell,  for  an  isolation  diode,  and  for 
the  series  resistance  that  is  typical  of  these  in  the  buoy  solar  array. 
The  data  are  projected  from  the  component  level  into  the  electrical 
configuration  determined  by  the  Design  Synthesis  driver  program.  Equa- 
tions are  also  included  to  estimate  performance  when  the  array  is  mis- 
oriented  from  the  sun  vector  and  to  estimate  performance  degradation  due 
to  cloud  cover,  temperature  extremes,  and  environmental  effects. 

PROGRAM  ALGORITHMS 


Step  1 Obtain  Elapsed  Time  From  Start  of  Mission 

DATEM  = Elapsed  time  from  start  of  mission  - days 

Step  2 Obtain  Current  Degradation  Factors  for  Solar  Array 

CDEGA  = Solar  Array  Current  Degradation  Factor  Due  to 
Fabrication  Losses  - Percent  (from  zero) 

CDEGB  = Solar  Array  Current  Degradation  Factor  Due  to  Terres- 
trial Performance  Extrapolation  Uncertainty  - 
Percent  (from  zero) 


Step  3 Calculate  Current  Degradation  Factor  Due  to  Environmental  Effects 
ODEGC  « SADEGC( DATEM) 

Where:  CDEGC  = Solar  Array  Current  Degradation  Factor  Due  to 
Environmental  Effects  - Percent  (from  zero) 

SADEGC  * Table  of  Solar  Array  Input  Current  Degradation 
Due  to  the  Environment  (in  Percent  from  zero) 
as  a Function  of  DATEM 
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Step  1*  Calculate  Solar  Array  Current  Degradation  Factor 

1.0  * 106  - (100.0  - CDEGA ) » ( 100.0  - CDEGB)  » ( 100.0  - CDEGC ) 

CDEU  = " ^ n " ” ■ i r r ■'  i i . 

1.0  * 10 

Where:  CDEG  = Solar  Array  Current  Degradation  Factor  - 
Dimensionless 

Step  5 Obtain  Voltage  Degradation  Factor  for  Solar  Array 

VDEGA  = Solar  Array  open  circuit  voltage  degradation  due  to 
temperature  uncertainty  - Percent  (from  zero) 

Step  6 Calculate  Voltage  Degradation  Factor  Due  to  Environmental  Effects 
VDEGB  = SADEGV(DATEM) 

Where:  VDEGB  = Solar  Array  Open  Circuit  Voltage  Degradation 
Factor  due  to  Environmental  Effects  - 
Percent  (from  zero) 

SADEGV  = Table  of  Solar  Array  Open  Circuit  Voltage 

Degradation  due  to  the  Environment  (in  percent 
from  zero)  as  a function  of  DATEM 


Ste£_I 


Step  8 


Calculate  Solar  Array  Voltage  Degradation  Factor 


VDEG 


1.0  » 10  - (100.0  - VDEGA)  * ( 100.0  - VDEGB) 

i.o  • io1* 


Where:  VDEG  * Solar  Array  Voltage  Degradation  Factor  - 
Dimensionless 


Obtain  Solar  Cell  Spectral  Correction  Factor 

SPEC0R  = Solar  Cell  Spectral  Correction  Factor  - Deminsionless 
(Corrects  for  differences  between  Spectrum  of  Solar 
Radiation  Incident  on  Solar  Cell  and' Spectral  Response 
of  Solar  Cell) 


Step  9 Obtain  Total  Solar  Radiation 

QDT  ■ Total  Solar  Radiation  Incident  on  Solar  Array  - 
Watts/Meter^ 
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r 

"top  10  Calculate  Effective  Solar  Insolation 
X = SPECOR  * QDT/10.0 

Where:  X = Effective  Solar  Insolation  Incident  on  Solar  Cell  - 
Milliwatts/cm^ 

Step  11  Calculate  Modified  Solar  Insolation 
XX  = X * (1.0  - CDEG) 

p 

Where:  XX  = Modified  Solar  Insolation  - mw/em 

Step  12  Obtain  Single  Solar  Cell  Area 

ACELL  = Single  Solar  Cell  Area  - cm^ 

Step  13  Obtain  Solar  Array  Temperature 

TSAC  = Solar  Array  Temperature  - °C 

Step  ll*  Calculate  Short  Circuit  Current  Temperature  Coefficient  for  a 
Single  Solar  Cell 

ALPHA C * (7.1*28  * 10"7  - (1.83  * 10"9)  * TSAC  * (XX)  * ACELL/L.O 

Where:  ALPHAC  = Short  Circuit  Current  Temperature  Coefficient  - 
Amperes/°C-cell 

Step  15  Calculate  Solar  Cell  Series  Resistance 

RCELLC  * F[RSCELL,  TEMTAB]  at  TEMTAB  = TSAC 

Where:  RCELLC  ■ Solar  Cell  Series  Resistance  - Ohms 
(at  Temperature  TSAC) 

RSCELL  ■ Internal  Table  of  Soiar  Cell  Series  Resistance 
as  a function  of  Cell  Temperature 

TEMTAB  * Internal  Table  of  Temperature  Range  Associated 
with  RSCELL 


3-64 


Step  16  Calculate  Solar  Cell  I-V  Curve  Correction  Factor 
ROCELL  = F[ROE,  SUNLIT]  at  SUNLIT  = XX 


Where:  ROCELL  = Solnr  Cell  I-V  Curve  Correction  factor  at  Solar 
Insolation  Level:  XX 

ROE  = Internal  Table  of  So’ar  Cell  I-V  Curve  Correc- 
tion Factor  as  a Function  of  Solar  Insolation 

SUNLIT  = Internal  Table  of  Solar  Insolation  Range 
Associated  with  ROE 


Step  IT  Calculate  Open  Circuit  Voltage  Temperature  Coefficient  for  a 
Single  Solar  Cell 

BETAA  = F [BETAB(or  BETAC  or  BETAD)]  at  XX  and  TSAC 
BBETA  = BETAA/ 1000.0 

Where:  BBETA  * Open  Circuit  Voltage  Temperature 

Coefficient  - (Volts/°C)  at  XX 
and  TSAC 

BETAA  ■ Open  Circuit  Voltage  Temperature 
Coefficient  - (mv/°C)  at  XX 
and  TSAC 

BETAB,  BETAC,  BETAD  * Internal  Tables  of  Solar  Cell  Open 

Circuit  Voltage  As  a Function  of 
Solar  Insolation  and  Cell 
Temperatures 

SUNMW,  SONMW,  SENMW  * Internal  Tables  of  Solar  Insola- 
tion Ranges  Associated  with  (BETA) 
Tables 

BTEMP,  CTEMP,  DTEMP  ■ Internal  Tables  of  Solar  Cell 

Temperature  Ranges  Associated  with 
(BETA)  Tables 

Internal  Tables  BTEMP,  SUNMW  AND  BETAB  used  when: 

(100  <_  XX  <_  5^0  mw/cm2)  and  (-60  < TSAC  <_  l60°C) 

Internal  Tables  CTEMP,  SONMW,  BETAC  used  when: 

(5  1 XX  < 253  mw/cm2)  and  (-1*0  <.  TSAC  <_  60°C) 

Internal  Tables  DTEMP,  SENMW,  BETAD  used  when: 

(5  < XX  < 100  mw/cm2)  and  (-lUO  < TSAC  < -1*0°C) 
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Step  18  Obtain  Single  Cell  ISC,  VOC  Data 

IISC  = Solar  Cell  Short  Circuit  Current  - Amperes/cell 
(at  lU 5 mv/cm^  Solar  Insolation  and  60°C) 

WOC  ■ Solar  Cell  Open  Circuit  Voltage  - Volts/cell 
(at  ll*5  mv/cm2  Solar  Insolation  and  60°C) 

Step  19  Calculate  ISC,  VOC  Shift  Due  to  Degradation 
Cl  = CDEG  * IISC 

C2  = VDEG  * WOC 

Where:  Cl  = Solar  Cell  Short  Circuit  Current  Shift  - Amps/cell 
C2  = Solar  Cell  Open  Circuit  Voltage  Shift  - Volts/cell 

Step  20  Obtain  Single  Circuit  (of  Solar  Cells)  Arrangement 
NS  = No  of  Solar  Cells  in  Series  in  Each  Circuit 

NP  = No  of  Solar  Cells  in  Parallel  in  Each  Circuit 

Step  21  Calculate  Cell  Electrical  Circuit  Parameters 

ALPHA  = ALPHAC  * NP 
BETA  » BBETA  * NS 
RCELL  >»  (O.lll*  ♦ RCELLC ) * NS/NP 
RHO  « ROCELL  * NS/NP 

Where:  ALPHA  ■ Short  Circuit  Current  Temperature 
Coefficient  for  a Single  Circuit  - 
Amperes/°C-circuit 

BETA  ■ Open  Circuit  Voltage  Temperature 
Coefficient  for  a Single  Circuit  - 
Volts/#C 

RCELL  * Single  Circuit  Series  Resistance  - Ohms 

RHO  ■ Series  Resistance  Temperature  Correction 
Factor 
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Step  22  Calculate  Modified  Electrical  Circuit  Short  Circuit  Current 
ISC  = IISC  * NP  * (1.0  - CDEG) 

Where:  ISC  = Modified  Electrical  Circuit  Short  Circuit 
Current  - Amperes/circuit 


Step  23  Calculate  Short  Circuit  Current  Difference  (for  an  Electrical 
Circuit ) 

DISC  = ISC  * ( (X/1U5.0)  - 1.0)  + ALPHA  * (TSAC  - 60.0) 

Where:  DISC  = Short  Circuit  Current  Difference  due  to  current 
degradation,  solar  insolation  changes  and  temp- 
ature  changes  - Amperes/circuit 


Step  2L  Calculate  Electrical  Circuit  Voltage  and  Series  Resistance 
Correction  Factors 

C3  ■ BETA  * (TSAC  - 60.0)  ♦ DISC  * RCELL 
CU  = RHO  * (TSAC  - 60.0) 

Where:  C3  * Electrical  Circuit  Voltage  Correction  Factor  - 
Volt8/circuit 

CU  * Electrical  Circuit  Series  Resistance  Correction 
Factor  - Ohms 


Step  25  Obtain  Reference  Solar  Cell  Current-Voltage  Characteristics 


IJ(J)  ■ Reference  Solar  Cell  Current  Data  point  - 
Amperes 

VV(J)  ■ Reference  Solar  Cell  Voltage  Data  point  - 
Volts 


Internal 

Tables 


Where:  J ■ 1,30 


Step  26  Calculate  Solar  Cell  Electrical  Circuit  Current-Voltage 
Characteristics 


I(J)  ■ NP  • (II(J)  - Cl)  ♦ DISC 

V(J)  ■ NS  • (W(J)  - C 2)  - C3  - (CU  • I(J)) 


3-67 


50^0-27 


d 

J = 1,30 

I(J)  = Electrical  Circuit  Current  - Amperes  at  the  given 
level  of  V(J) 

V ( J ) = Electrical  Circuit  Voltage  - Volts 

Step  27  Obtain  Solar  Array  Voltage  Increment 

VSAINC  = Solar  Array  Voltage  Increment  - volts 

Step  28  Redefine  Electrical  Circuit  Current  - Voltage  Array  in  Selected 
Voltage  Increments  as  follows: 

a)  Set:  Counter  L = 1 and  voltage  V2(L)  - 0.0 

b)  Establish:  Current  Il(L)  at  V2(L) 

I1(L)  * F U(J),  V(J)}  at  V(J)  = V2(L) 

c)  Increment:  Counter  L * L ♦ 1 and  voltage  V2(L  •*■!)=  V2(L)  ♦ 
VSAINC  and  Establish:  Current  Il(L  ♦ l)  at 
V2(L  ♦ 1), 

Until:  I1(L  1)  <,  0.0 

d)  Redefine:  Last  V2(L)  at  Il(L)  * 0.0 
V2(L)  * F {I(J),  V(J)}  at  I(J)  = 0.0 

e)  Set:  Current-Voltage  Matrix  Dimension  to  last  counter 
value:  MFINAL  ■ L 

Step  29  Obtain  Number  of  Solar  Cell  Electrical  Circuits  in  Solar  Array 

NESP  ■ Number  of  Electrical  Circuits  in  Solar  Array 
(assumed  in  parallel) 
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Step  30  Calculate  Solar  Array  Current-Voltage  Characteristics 
I2(L)  = ( II ( L ) • NESP)  at  V2(L) 

L = 1,  MFINAL 

Where:  Il(L)  = Electrical  Circuit  Current  - Amperes  at  V2(L) 
12 ( L ) = Solar  Array  Current  - Amperes  at  V2(L) 

V2(L)  = Circuit  or  Array  Voltage  - Volts 


Obtain  Voltage  Data  for  Calculation  of  Solar  Array  Maximum  Power 
Point 

XV  = Initial  Voltage  for  Max  Power  Point  Calculations  - Volts 

DXN  s Voltage  Increment  for  Maximum  Power  Point  Calculation  - 
Volts 


Initialize  Calculation  Value  of  Solar  Array  Maximum  Power 
MSAPWR  = 0.0 

Where:  MSAPWR  ■ Solar  Array  Maximum  Power  - Watts 


Step  33  Calculate  Solar  Array  Power  and  Current 
XI  « F { I 2 { L ) , V2( L ) } at  V2(L)  ■ XV 
SAPWR  « XI  • XV 

Where:  XV  « Solar  Array  Voltage  - Volts 

XI  ■ Solar  Array  Current  - Amperes 
SAPWR  ■ Solar  Array  Power  - Watts 


Step  3^  Compare  Solar  Array  Power  With  Maximum  Power 

IF:  SAPWR  > MSAPWR 
THEN:  MSAPWR  ■ SAPWR 


XV  ■ XV  ♦ DXV 


Repeat  Step  33  until:  SAPWR  < MSAPWR 
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Recalculate  Solar  Array  Current  and  Power 
MSAPWR  =0.0 


Step  j5 

XV  = XV  - DXV 
REPEAT  STEP  33  ONLY 

Step  36  Compare  Solar  Array  Power  With  Maximum  Power 

IF:  ■ SAPWR  _>  MSAPWR, 

THEN : MSAPWR  = SAPWR 

DXV  = DXV/ 10.0 

XV  = XV  ♦ DXV 

REPEAT  STEP  33  ONLY  UNTIL:  SAPWR  < MSAPWR 

Step  37  Calculate  Solar  Array  Maximum  Power  Point  Characteristics 
MAXV  = XV  - DXV 
MAXI  = MSAPWR /MAXV 

Where:  MAXV  * Solar  Array  Voltage  at  Max  Power  Point  - Volts 

MAXI  = Solar  Array  Current  at  Maximum  Power  Point  - 
Amperes 

RETURN  TO  DESIGN  SYNTHESIS  DRIVER  ROUTINE 
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3.1*  Power  Conditioning  and  Distribution  Group 

The  Power  Conditioning  and  Distribution  Group  is  made  up  of  the  lamp 
flasher  and  the  housekeeping  regulator.  The  characteristics,  of  these  two 
subassemblies  are  computed  as  a function  of  the  selected  lamp  flasher 
pattern  and  the  lamp  flasher  condition  (on,  off,  or  flashing).  These 
characteristics  are  then  shifted  for  the  combined  effects  of  wiring  and 
connector  series  resistance  to  give  a single  set  of  current-voltage  curves 
at  the  unregulated  bus. 


PROGRAM  ALGORITHMS 


Step  1 Obtain  Flasher  Pattern  Type 

IF:  IFTYPE  « 0. 

THEN:  GO  TO  STEP  3 

Where:  IFTYPE  is  the  type  of  flasher  pattern 
= 0:  Non-Standard  Pattern 
> 0:  Standard  Pattern 


Step  2 


Calculate  Standard  Flasher  Pattern 


TL1(J)  ■ TL0  (IFTYPE,  J} 

(1  <_  IFTYPE  <_  15) 

(1  < J < l6) 


GO  TO  STF.P  4 


115  standard  pattern  types 
Up  to  16  steps  per  pattern 
Alternate  On/Off  steps 


Where:  TL0  is  an  input  table  containing  TL1  as  a function  of 
IFTYPE  and  J 

TL1  B Selected  Lamp  Flasher  Pattern 


Step  3 Calculate  Non-Standard  Flasher  Pattern 
Tll(J)  * TLL1(J)  (1  < J < 16) 

Where:  TLLl(J)  is  the  input  data  containing  up  to  1 6 alternate 
on-off  steps  for  the  Non-Standard  Flasher  Pattern 
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Step  U 


st*E  > 


Step  6 


Step  7 


Calculate  Total  Duration  of  Lanp  Illumination  and  Lamp  Shut-Off 
15 

TLON  = [tLI(J)] 

J-1,3,5... 


TLOFF  = 


E 


Ctli(j)] 


Where:  TLON  3 Total  duration  of  lamp  illumination 

TLOFF  * Total  duration  of  lamp  shut-off 


in  a single 
flasher  period 


IF:  TLON  < 0,  and;  TLOFF  < 0 

THEN:  Stop  program  and 

Print:  "No  flasher  pattern  entries" 


Calculate  Lamp  Duty  Cycle 
DL  - TLON /(TLON  ♦ TLOFF) 
Where:  DL  ■ Lamp  Duty  Cycle 


Obtain  Lamp  Characteristics 

VLR  ■ Lamp  Voltage  Rating  - VDC 

ILR  ■ Lamp  Current  Rating  - Amperes 

CLS  ■ Cola-Filament  Lamp  Surge  Coefficient 


Calculate  Actual  Lamp  Current 
IL  » CLS  • ILR 

Where:  IL  ■ Actual  Lamp  Current  - Amperes 


Note:  If  DL  » 1.0  then  CLS  ■ 1.0 
If  DL  < 1.0  then  CLS  > 1.0 
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Step  8 Calculate  Actual  Lamp  Resistance 
VLR 


RL  = 


IL 


Where:  RL  * Actual  Lamp  Resistance  - Ohms 


Step  9 Calculate  Average  Lamp  Current 
IL  = IL  * DL 

Where:  IL  = Average  Lamp  Current  - Amperes 

Step  10  Calculate  Effective  Lamp  Resistance 
RL  = VLR/<IL) 

Where:  RL  = Effective  Lamp  Resistance  - Ohms 

Step  11  Obtain  Raw  Power  Bus  Voltage  Limits  and  User  Load  Cable  Resistance 
VMINIV  * Minimum  Raw  Power  Bus  Voltage  - VDC 
VMAXIV  = Maximum  Raw  Power  Bus  Voltage  - VDC 
RLL  * User  Load  Cable  Resistance  - Ohms 

Step  12  Calculate  PCD  Group  Voltage  Increment 
VINCIV  = {VMAXIV  - VMINIV)/50.0 
Where:  VINCIV  * PCD  Group  Voltage  Increment  - VDC 

Step  13  Obtain  PCD  Equipment  Temperature  Characteristics 
TTAMB  * Ambient  Temperature  - °F 
DTTPCD  * PCD  Equipment  Temperature  Rise  - °F 
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Step  lk  Calculate  PCD  Equipment  Temperature 
TTPCD  = TTAMB  + DTTPCD 

Where:  TTPCD  = PCD  Equipment  Temperature  - °F 

Step  13  Compare  Raw  Power  Bus  Minimum  Voltage  With  Reference 

VRIO  = VRIOT  {TTPCD} 

VRISA  = VRISAT  {TTPCD} 

IF:  VMINIV  < VRIO, 

THEN:  GO  TO  STEP  16 

IF:  (VMINIV  > VRIO),  And:  (VMINIV _<  VRISA) , 

THEN:  Go  to  Step  2k 

IF:  VMINIV  > VRISA, 

THEN:  GO  TO  STEP  29 

Where:  VRIO  = Minimum  (No  Current)  Voltage  Drop  - VDC 
Across  Lamp  Regulator  in  "Saturated" 

Condition 

VRISA  = Voltage  level  at  which  lamp  regulator  - VDC  changes 
from  "Saturated"  condition  operation  to  "Active" 
operation 

VRIOT  = Input  Table  of  VRIO  as  a function  of  TTPCD 
VRISAT  = Input  Table  of  VRISA  as  a function  of  TTPCD 

Step  16  Initialize  Counter  and  Lamp  Regulator  Voltage 
J = 1 

VRI(J)  = VMINIV 

Step  17  Calculate  Lamp  Regulator  Current 
IRI(J,1)  = 0.0 
IRI(J,2)  = 0.0 
IRI(J,3)  * 0.0 
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Step  IT  (Contd) 

Where:  VRI(J)  = Lamp  Regulator  Input  Voltage  - VDC 

IRI(J,K)  = Lamp  Regulator  Input  Current  - Volts 

When:  K=1  - Lamp  Off 

K=2  - Lamp  Flashing  - Effective 
K=3  - Lamp  On 


Step  l8  Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare 
With  Reference 

J = J + 1 

VRI(J)  = VRI(J-I)  + VINCIV 

(VRI(J)  > VRIO)  And: 

(VRI(J)  < VMAXIVI 
GO  TO  STEP  20 
VRI(J)  > VMAXIV 
GO  TO  STEP  32 

Step  19  Calculate  Lamp  Regulator  Currents 
IRI( J,l)  = 0.0 
IRI(J,2)  = 0.0 
IRI( J,3)  = 0.0 
REPEAT  STEPS  18  AMD  19 

Step  20  Calculate  Lamp  Regulator  Current 
IRl(J,l)  = 0.0 

IRI(J,2)  = (VRI(J)-VRI0)/(RL  ♦ ZRS) 

IRI(J,3)  = ( VRI ( J ) -VRIO ) / ( RL  ♦ ZRS) 

Where:  ZRS  ■ Regulator  Impedance  in  ’’Saturated”  Condition  - Ohms 


IF: 

IF: 

THEN: 

IF: 

THEN: 


ZRS  * ZRST  {TTPCD} 

ZRST  * Input  Table  of  ZRS  as  a function  of  TTPCD 


Step  21  Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare  with 
Reference 

J = J + 1 

VRI(J)  = VRI( J - 1)  + VINCIV 

IF:  (VRI(J)  > VRISA)  AND: 

IF:  (VRI(J)  < VMAXIV ) 

THEN:  Go  to  Step  22 
IF:  VRI(J)  > VMAXIV 

THEN:  GO  TO  STEP  32 

REPEAT  STEPS  20  AND  21 


Step  22  Calculate  Lamp  Regulator  Currents 
VLB  = VLBT  {VRI,  TTPCD} 

ZRA  = ZRAT  {TTPCD} 

IRI(J,1)  = 0.0 
IRI(J,2)  = Vlb/(MI  ♦ ZRA) 

IRI( J,3 ) = Vlb/(RL  ♦ ZRA) 

Where:  V = Regulator  Output  Voltage  at  Zero  Current  - Volts 
Ld 

ZRA  = Regulator  impedance  in  "Active"  region  - Ohms 
VLBT  = Input  Table  of  VLB  as  a function  of  VRI  and  TTPCD 
ZRAT  = Input  Table  of  ZRA  as  a function  of  TTPCD 


Step  23  Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare 
With  Reference 

J ■ J ♦ 1 

VRI(J)  = VRI( J - 1)  ♦ VINCIV 

IF:  VRI(J)  > VMAXIV 

THEN:  GO  TO  STEP  32 


REPEAT  STEPS  22  AND  23 


Step  2l* 


Initialize  Counter  and  Lamp  Regulator  Voltage 


Step  23 


Step  26 


Step  27 


J = 1 

VRI(J)  = VMINIV 


Calculate  Lamp  Regulator  Currents 
IRI( J,l)  = 0.0 

IRI(J,2)  = (VRI(J)  - VRIO ) /( RL  + ZRS) 
IRI(J,3)  = (VRI(J)  - VRIO )/ ( RL  + ZRS) 


Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare 
with  Reference 

J = J + 1 

VRI(J)  = VRI( J - l)  + VINCIV 

IF:  (VRI(J)  > VRISA),  AND: 

_IF:  (VRI(J)  < VMAXIV) 

THEN : c;o  TO  STEP  27 

IF:  VRI(J)  > VMAXIV 

THEN : GO  TO  STEP  32 

REPEAT  STEPS  25  AND  26 

Calculate  Lamp  Regulator  Currents 

IRI ( J ,1)  = O.G 

IRI(J,2)  = VLB/ ( RL  + ZRA) 

IRI(J,3)  * VLB/ ( RL  + ZRA) 
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Step  28  Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare  with 
Reference 

J = J + 1 

VRI(J)  = VRI( J - 1)  + VINCIV 

IF:  VRI(J)  > VMAXIV 

THEN:  GO  TO  STEP  32 

REPEAT  STEPS  27  AND  23 

Step  29  Initialize  Counter  and  Lamp  Regulator  Voltage 
J = 1 

VRI(J)  = VMINIV 

Step  30  Calculate  Lamp  Regulator  Currents 
IRI ( J,l)  * 0.0 
IRl(J ,2)  = VLB/(RL  + ZRA ) 

IRI(J,3)  = VLB/ ( RL  ♦ ZRA) 

Step  31  Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare  with 
Reference 

<J  * J + 1 

VRI(J)  - VRI( J - l)  ♦ VINCIV 

IF:  VRI(J)  > VMAXIV 

THEN:  Go  to  Step  32 

REPEAT  STEPS  30  AND  31 

Step  32  Calculate  PCD  Group  Current 

IHI(J)  ■ IHIT  {VRI(J),  TTPCD) 

XI( J,K)  « IHI(J)  ♦ IRI( J,K) 
for:  J ■ 1,  50;  K ■ 1,3 


Where:  Xl(J,K)  = PCD  Group  Current  - Amperes 

IHl(J)  = Housekeeping  Load  Regulator  Input  Current  - 
Amperes 

IHIT  = Input  Table  of  IHI(J)  as  a function  of  VRI(J) 
and  TTPCD 

Step  33  Calculate  PCD  Group  Voltage 

XX(J,K)  = VRI(J)  + Xl(J,K)  * RLL 

for:  J = 1,50;  K ■ 1,3 

Where:  XX(J,K)  = PCD  Group  Voltage  - VDC 
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4.  PERFORMANCE  ANALYSIS 


The  Performance  Analysis  segment  of  the  DSPA  program  uses  known  power 
system  arrangements,  electrical  sizes,  and  physical  characteristics  to 
determine  the  response  (operational  characteristics)  of  the  equipment  to  a 
given  stimulus  (load  and  environmental  profiles).  A block  diagram  of  the 
Performance  Analysis  driver  program  is  shown  in  Figure  4.1.  As  shown, 
the  object  of  the  performance  analysis  methodology  is  to  obtain  the  raw 
power  bus  operating  point  for  each  time  increment  during  a mission 
period.  Once  this  operating  point  is  obtained,  all  power  system  opera- 
tional characteristics  are  determinable. 

To  obtain  the  operating  point,  the  power  system  is  divided  into  groups. 
These  groups  and  the  equipment  comprising  them  are: 

a.  Power  Source  Group: 

1)  Solar  Array. 

2)  Shunt  Limiter. 

3)  Solar  Array  Isolation  Diodes. 

4)  Solar  Array  Cable. 

b.  Energy  Storage  Group: 

1 ) Batteries . 

2)  Battery  Chargers. 

3)  Battery  Isolation  Diodes. 

4)  Battery  Cables. 

c.  Power  Conditioning  and  Distribution  Group: 

1)  User  Loads. 

2)  Load  Cablo. 

Current-voltage  (I-V)  characteristics  are  determined  for  each  equipment 
group.  The  Power  Conditioning  and  Distribution  Group  characteristics 
are  then  deducted  from  those  of  the  Power  Source  Group  to  obtain  a 
Difference  Curve.  The  intersection  of  the  Difference  Cu.-ve  with  the 


4-1 


FIGURE  1»-1.  POWER  SYSTEM  PERFORMANCE  ANALYSIS 
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Energy  Storage  Group  characteristics  is  an  estimate  of  the  voltage  and 
current  on  the  raw  power  bus  (operating  point).  During  the  time  incre- 
ment under  examination,  the  battery  is  in  a state-of-charge,  discharge, 
or  open-circuit  depending  on  the  value  and  sign  of  the  operating  point 
current. 

After  the  operating  point  is  obtained,  the  operational  characteristics 
of  the  power  system  equipment  are  determined.  The  computer  program  is 
then  ready  for  examination  of  power  system  response  for  the  next  time 
increment.  The  process  is  repeated  until  the  end  of  the  mission  period* 
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PROGRAM  ALGORITHMS 


Step  1 

Obtain 

Pertinent  Mission  and  Equipment  Information 

QON 

= Solar  Insolation  Level  for  Lamp  Flasher  Turn-On  - 
Watts/Meter^ 

QOFF 

= Solar  Insolation  Level  for  Lamp  Flasher  Turn-Off  - 
Watts /Met er^ 

ISH 

= Shunt  Limiter  Type 

Where: 

0 = None 

1 = Ordinary  Zener  Diode 

2 = Temperature  Compensated  Zener  Diode 

3 = Active  Shunt  Limiter 

VBUSI 

* Initial  Estimate  of  Power  System  Operating  Point  Voltage  at 

Raw  Power  Bus  - VDC 

Step  2 

Obtain 

Free  Format  Data  Head-In  Card  Type 

NCTYPE 

= Free  Format  Card  Type 

Where : 

0 = Start-Up  Data 

1 * Time-Variant  Data 

<0  = Termination  Data 

IF: 

NCTYPE  » 0 

THEN: 

GO  TO  STEP  3 

IF: 

NCTYPE  ■ 1 

THEN 

GO  TO  STEP  6 

IF: 

NCTYPE  < 0 

THEM: 

Terminate  Performance  Analysis  Calculations  and  return  to 
main  program 

Step  3: 

Obtain 

Free  Format  Start-Up  Data 

START 

■ Starting  "Time"  of  Computer  Run  - YYDDDHHMM 

Where: 

YY  « Year  {00-99} 

DDD  ■ Date  { 001—365 } days  from  beginning  of  the  year 

HH  ■ Hours  {00-24}  Hours  after  Midnight 

MM  * Minutes  {00-60}  minutes  after  the  hour 
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Step  3 (Contd) 

HLLA  = Nominal  Time  Increment  for  Performance  Analysis 
Calculations  - hours  (Default  value  =1.0  hours) 

ACCQB  = Battery  State-of-Charge  Accuracy  Requirement  for  the 
Predictor-Corrector  Routine  - (Default  value  = 0.01) 

CT  = Cloud  Type  (at  Start-Up) 

0.0  = Cirrus  or  Cirrostratus  Clouds 

1.0  = Stratus  Clouds 

2.0  = Other  Cloud  Types 

TC  = Cloud  Cover  (at  Start-Up) 

0.0  “ No  Cloud  Cover 

(0.0  < TC  £ 10.0)  * Tenths  of  Sky  Covered  by  Clouds 
INDFLS  = Lamp  Flasher  Condition  Indicator 

0 = Lamp  Flasher  Off  (Lamp  Not  Flashing) 

1 = Lamp  Flasher  On  (Lamp  Flashing) 

Step  1*  Calculate  Starting  "Time"  Information  for  Computer  Run 

YEAR  = YY  {00-99}  - Year 

DATE  ■ DDD  {001-365)  - Days  from  beginning  of  year 

TIMEH  « HH  ♦ - Hours  after  Midnight 

DvJ«  U 

DATEM  * 0.0  - Elapsed  Time  Since  "Start"  of  Computer  run  - days 
Step  5 Initialize  "Starting"  Time  Reference  Data 
YEAR1  « YEAR 
DATE1  ■ DATE 


DATEM1  « DATEM 


Go  To  Step  10 


Step  6 Obtain  Free  Format  Time-Varying  Data 


DURA  * Duration  of  a Significant  Time  Interval  (=  DDDHHMM) 

Where:  DDD  = Days  (>_  000) 

HH  = Hours  (>  00) 

MM  = Minutes  (>  00) 

NTS  a Number  of  time  steps  during  Significant  Time  Interval  (DURA) 
for  printing  out  Performance  Analysis  Data 

CT  = Cloud  Type  (During  Significant  Time  Interval) 

TC  » Cloud  Cover  (During  Significant  Time  Interval) 

Step  7 Initialize  Time  Step  Counter  and  Durations  of  Various  Intervals 


LNTS  « 1 


DURAH  = (24. C).  * DDD  ♦ HH  + 


HINT  ■ DURAH/NTS 


Where:  DURAH  » Duration  of  & Significant  time  interval  - Hours 

HINT  = Duration  of  a print-out  time  step  - Hours 
during  a significant  time  interval 

Step  7a  Initialize  Reference  Value  of  Print-Out  Time  Step  Duration 
H2  ■ 0.0 

Where:  H2  » Reference  Value  of  ?rint-0ut  Time  Step  Duration  - Hours 
Step  8 Calculate  Performance  Analysis  Calculation  Time  Interval 
H - AMIN  ( HLLA-H2 , HINT) 


Where:  H ■ Performance  Analysis  Calculation  Time  Interval  - Hours 
Step  9 Calculate  "Time"  of  Performance  Analysis  Calculation 
DATEM  - DAT Q41  ♦ (H/24.0) 

DATE?  • DATE1  ♦ (H/24.0) 

TIMEH  ■ (DATE2  - IFIX  (DATE?))  • 24.0 
DATE  - IFIX  (DATE2) 
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Step  9 (Contd) 

IF:  DATE  > 365. 0 

THEN:  DATE  = DATE  - 365. 0,  AND, 

THEN:  YEAR  = YEAR1  +1.0 

Step  9a  Initialize  Low  Insolation  Load  Selector 

KLL  3 3 

Where:  KLL  * Load  Selector  at  Low  Values  of  Solar  Insolation 


1 * Lamp  Off 
3 * Lamp  On 

Step  10  Obtain  Power  Sources  Group  Characteristics  at  the  Raw  Power  Bus 

NESP  * Number  of  Electrical  Circuits  in  Solar  Array 

RSA  * Solar  Array  Electrical  Circuit  Cable  Resistance  - Ohms 

MSAPWR  * Solar  Array  Maximum  Power  - Watts 

QDT  3 Total  Solar  Insolation  Incident  on  Solar  Array  - 
Watts/Meter2 

SX( I ) 3 Power  Source  Group  Voltage  - VDC 
SY(I)  3 Power  Source  Group  Current  at  SX(I)  - Amperes 
I ■ 1,  NAPSG 

Where : IF:  ISH  ■ 0;  THEN:  NAPSG  3 MFINAL 

IF:  ISH  > 0;  THEN:  NAPSG  ■ NFINAL 

IF;  QDT  < 0.0;  THEN:  NAPSG  3 KFINAL 

MFINAL  OR  NFINAL  are  calculated  as  part  of  the  Power  Source  Group 
current-voltage  estimate 

KFINAL  is  the  maximum  extent  of  the  SX(l)  and  SY(I)  arrays. 

Step  11  Determine  Lamp  Flasher  Condition 

IF:  QDT  < 0.0 

THEN:  INDFLS  • 1 AND:  KL  ■ 2 

IF:  (0.0  < QDT  < QON) 

THEN:  INDFLS  • 1 AND:  FL  ■ 3 

IF:  (QON  < QDT  < QOFF)  AND:  INDFLS  • 0 

THEN:  INDFLS  ■ 0 AND:  KL  - 1 
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Step  11  (Contd) 


IF: 

THEN: 

(QON  < QDT 
INDFLS  = 1 

< QOFF)  AND:  INDFLS  = 1 
AND:  KL  = 3 

IF: 

THEN: 

QDT  > QOFF 
INDFLS  = 0 

AND: 

KL  = 1 

IF: 

THEN: 

NCTYPE  = 0 
KL  = 2 

AND: 

KL  = 3 

Where: 

KL  = Load 

Selector  Indicator 

1 = Lamp 

Off 

(During  Daylight) 

2 - Lamp  Flashing  (At  Night) 

3 = Lamp  Off/On  (At  Low  Solar  Insolation  Levels) 

Step  12  Determine  Power  Conditioning  and  Distribution  Group  Level  Selector 

IF:  KL  = 1 

THEN:  K * 1 

IF:  KL  = 2 

THEN:  K = 2 

IF:  KL  = 3 AN£:  KLfj  * 1 

THEN:  K = 3 

Where:  K = Power  Conditioning/Distribution  Group  Load  Selector 

1 = Lamp  Off 

2 * Elfective  Load  - Lamp  Flashing 

3 * Lamp  On 

Step  13  Obtain  Power  Conditioning  and  Distribution  Group  Characteristics  at 
the  Raw  Power  Bus 

XX(J)  ■ Power  Conditioning  and  Distribution  Group  Voltage  - VDC 

XI (J ) ■ Power  Conditioning  and  Distribution  Group  Current  at  XX(J)  - 

Amperes 

J ■ 1,50  (Number  of  Data  Points) 

K ■ 1,  2,  or  3 (See  Above) 

DL  ■ Lamp  Duty  Cycle 


w 

\ 

' 


50140-27 


Step  1 U Rearrange  Voltage  Data  into  One  Array  in  Ascending  Order 

DIFIV(L.l)  * {F  SX(I),  and  XX(J)} 

(in  ascending  order) 

1=1,  NAPSG 

J - 1,  50 

L * 1,  (50  + NAPSG) 

Where:  DIFIV(L,l)  « Difference  Curve  Voltages  - VDC 
Step  15  Calculate  Difference  Curve  Current  Values 
DIFIV(L,2)  = SY(I)  - XI(J) 

Where:  DIFIV(L,2)  * Difference  Curve  Current  at  Voltage  Level 

DIFIV(L.l)  - Amperes 

SY(I)  and  XI(J)  are  selected  at  voltage  DIFIV(L,1) 

Note:  If  DIFIV(L,1)  are  beyond  the  limits  of  one  of  the  current 
arrays,  then  use  extrapolation  methods  based  on  existing 
current  data. 

Step  l6  Obtain  Energy  Storage  Group  Characteristics  at  the  Raw  Power  Bus 
ICHRT  = Battery  Charger  Type 

0 ■ No  Charger 

1 ■ Constant  Voltage  Charger  with  Current  Limit 
NBATT  ■ Number  of  Batteries  in  Parallel 

TRESLT(LY,l)  ■ Energy  Storage  Group  Voltage  - VDC 

TRESLT(LY,2)  ■ Energy  Storage  Group  Current  - Amperes 
at  TRESLT(L.l) 


IF: 

ICHRT  » 

0;  THEN; 

LY  ■ 

1,9 

IF: 

ICHRT  ■ 

1;  THEN: 

LY  ■ 

1,10 

CB  ■ Maximum  Discharge  Capacity  - Ampere-Hours 

(for  each  Battery)  at  a standard  discharge  rate  to  a 
standard  terminal)  voltage  per  cell  at  a standard 
temperature. 


4 


50U0-27 


Step  16  (Contd) 

Standard  Discharge  Rate 

Standard  Min  Cell 
Terminal  Voltage 

Standard  Temperature 


NiCd  Batteries 

Pb-Acid  Batteries 

c/2 

c/20 

0.5 

1.0 

?0°F 

70°F 

TTESG  = Energy  Storage  Group  Temperature  - °F 

WB  = State  of  Charge  of  each  battery  - dimensionless 

Step  17  Initialize  Power  System  Operating  Point 

IF:  NCTYPE  = 0, 

THEN:  VBUS  * VBUSI 

Where:  VBUS  = Power  System  Voltage  at  Operating  Point  - VDC 
(at  Raw  Power  Bus) 

Step  18  Calculate  Current  Difference 

DELCUR  = DIFIV(M,2)  - TRESLT(M,2) 

Where:  DIFIV(M,2)  *=  F{DIFIV(L,2) , DIFIV(l.l)}  at  DIFIV(L,1)  - VBUS 

TRESLT(M,2)  - F {TRESLT(LY,2) , TRESLT(LY.l)} 
at  TRESLT(LY.l)  ■ VBUS 

DELCUR  ■ Difference  in  current  level  between  difference 
curve  and  ener^  storage  group  - Amperes 

Step  19  Compare  Current  Difference  with  Reference  and  Calculate  New  Operating 
Point  Voltage 

IF:  DELCUR  < 0.0005  (Amperes) 

THEN:  Go  to  Step  20 

VBUS2  « VBUS 

VBUS  ■ VBUS  ♦ DVBUS 

VBUS 3 ■ VBUS 

Where:  DVBUS  * Increment  added  to  operating  point  voltage  estimate 
by  Nevton-Raphaon  closure  routine  - VDC 

VBUS2,  VBUS 3 » Estimates  of  VBUS 

REPEAT  STEPS  1 8 AND  19 
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Step  20  Determine  Operating  Point  Stability 

IF:  SDIF  < SESG 

THEN:  GO  TO  STEP  21 

Where:  SESG  = Slope  of  the  Energy  Storage  Group  Current-Voltage 
Characteristics  (TRESLT)  at  VBUS  - Amperes/VDC 

SDIF  = Slope  of  the  Difference  curve  (DIFIV)  Current-Voltage 
Characteristics  at  VBUS  - Amperes/VDC 

VBSAVE  * VBUS 

Compare  VBUS3  to  VBUS2  and  continue  in  the  same  general  "direction" 
(from  VBUS2  to  VBUS3J  with  respect  to  voltage  (by  repeating  steps  18, 
19,  20)  until  a stable  point  is  found  or_  until  the  appropriate  lower 
or  upper  voltage  limits  (based  on  Difference  Curve  Voltages)  are 
reached.  If  no  stable  point  is  obtained,  set  VBUS  * VBSAVE  and 
repeat  steps  18,  19,  20  in  the  opposite  direction  (relative  to 
voltage)  until  a stable  point  is  obtained  or  until  the  appropriate 
lower  or  upper  voltage  is  reached. 

If  no  stable  operating  point  is  estimated,  then  terminate  run  and 
print  cut  the  following  information: 

DIAGNOSTIC: 

"RUN  TERMINATED  - HO  STABLE  OPERATING  POINT" 

"TRESLT(LY,1)  ■" 

"TRESLT(LY,2)  ■" 

LY  ■ 1,  10  to  (10  * NBATT) 

"DIFIV(L,1 ) »" 

"DIFIV(L,2)  «" 

L ■ 1,  (50  ♦ NAPSG ) 

Step  21  Calculate  Energy  Storage  Group  Current 

XITT  = DIFIV(M,2)  at  DIFIV(M,1)  ■ VBUS 
Where:  XITT  ■ Energy  Storage  Group  Current  - Amperes  at  VBUS 
Step  22  Obtain  Energy  Storage  Unit  Current-Voltage  Characteristics 

TRESV(LY)  ■ Energy  Storage  Unit  Voltage  - VDC 

TRESI(LY)  ■ Energy  Storage  Unit  Current  - Amperes 
at  TRESV(LY) 

IF:  ICHRT  ■ 0;  THEN:  LY  ■ 1,  9 
JF:  ICHRT  ■ 1;  THEN:  LY  ■ 1,  10 
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Step  23  Calculate  Battery  Current 

BCUR  = F (TRESI(LY),  TRESV(LY)}  at  TRESV(LY)  = VBUS 
Where:  BCUR  = Current  for  each  Battery  - Amperes 
♦ = Battery  Charging 
0 = Battery  Open  Circuit 
- = Battery  Discharging 

Step  2l*  Obtain  Battery  Current-Voltage  Characteristics 

VG(LJ)  * Battery  Potential  - VDC 

XIB(LJ)  * Battery  Current  - Amperes 
at  VB(LJ) 

LJ  = 1,  9 (Number  of  data  points) 

Step  25  Calculate  Battery  Potential 

VBAT  « F{VB(LJ),  XIB(LJ)}  at  XIB(LJ)  = BCUR 
Where:  VBAT  * Potential  of  each  battery  - VDC 
Step  26  Calculate  Power  Conditioning  and  Distribution  Croup  Current 
XIPCD  - F {XI(J),  XX(J))  at  XX(J)  • VBUS 
J ■ 1,  50  (data  points) 

K » 1,  2 or  3 (from  Step  12) 

Where:  XIPCD  ■ Input  Current  to  Power  Conditioning  and  Distribution 
Group  - Amperes 

Step  27  Calculate  Power  Sources  Group  Current 

XIPSG  • F{SY(I),  SX(I)}  at  SX(I)  ■ VBUS 
I ■ 1,  NAPSO 

Where:  XIPSG  ■ Power  Sources  Group  Output  Current  - Amperes 
at  the  Operating  Point 
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Step  28 


Step  29 


Step  30 


Step  31 


Step  32 


StefcJ * 


Obtain  Shunt-Limiter  Current-Voltage  Characteristics  at  the  Raw  Power 
Bus 

ZV(IZ)  = Shunt  Limiter  Voltage  - VDC 

ZI(IZ)  = Shunt  Limiter  Current  - Amperes 
at  ZV(IZ) 

IZ  = 1,  20  (data  points) 

Calculate  Shunt  Limiter  Current 

XIZ  = F{ZI( IZ ) , ZV(IZ)}  at  ZV(IZ)  = VBUS 

Where:  XIZ  = Shunt  Limiter  Current  at  the  operating  point  - 
Ampe  res 

Calculate  Solar  Array  Current 
XISA  * XIZ  ♦ XIPSG 

Where:  XISA  = Solar  Array  Current  - Amperes 
Calculate  Solar  Array  Electrical  Circuit  Current 
XI EC  • XISA/NESP 

Where:  XIEC  * Solar  Array  Electrical  Circuit  Current  - Amperes 

Calculate  Solar  Array  Potential 

VDIODE  • ADI  (XIEC) 

VSA  - VBUS  ♦ VDIODE  ♦ (XIEC  • RSA) 

Where:  VDIODE  ■ Electrical  Section  Blocking  Diode  Voltage  Drop  at 
XIEC  - VDC 

VSA  ■ Solar  Array  Potential  - VDC 
ADI  ■ Input  Table  of  VDIODE  as  a function  of  XIEC 
Calculate  Equipment  Power  Levels 
PESO  » | VBUS  • XITT| 

PBATT  ■ |VBAT  • BCUR| 

PPCD  ■ VBUS  • XIPCD 
PPSG  • VBUS  • XIPSG 
PSL  • VBUS  • XIZ 
PSA  ■ VSA  • XISA 
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Step  33  (Contd) 

Where:  PESO  = Energy  Storage  Group  Power  - Watts 

PBATT  = Power  Level  of  each  Battery  - Watts 

PPCD  = Power  Conditioning  and  Distribution  Group  Power  - 
Watts 

PPSG  = Power  Source  Group  Power  - Watts 
PSL  = Shunt  Limiter  Power  - Watts 
PSA  * Solar  Array  Power  - Watts 
Step  jk  Calculate  Equipment  Power  Margins 
MARSA  = MSAPWR  - PSA 

Wiiere:  MARSA  = Solar  Array  Power  Margin  - Watts 
Step  36  Initialize  Battery  Charge  Efficiency 
IF:  NCTYPE  * 0;  THEN:  ETA  * 1.0 

Where:  ETA  * Instantaneous  Charge  Efficiency  for  each  Battery  - 
dimensionless 


Step  j'f  Compare  Battery  Charge  Efficiency  With  Reference 

IF:  ETA  <_  0.0;  THEN:  ETA  ■ 0.00001;  AND:  Go  To  Step  k2 

Step  38  Compare  Battery  Current  With  Reference 

IF:  BCUR  < 0.0;  THEN:  ETA  ■ 1.0;  AND:  Go  To  Step  k2 
Step  39  Calculate  Normalized  Battery  Charge  Rate 
CHRN  ■ BCUR/CB 

Where:  CHRN  ■ Normalised  Battery  Charge  Rate  - Hour"1 
Step  to  Calculate  Instantaneous  Battery  Charge  Efficiency 
ETA  ■ A{TTESG,  CHRN,  QB) 

Where:  A ■ A serien  of  Input  Data  Tables  (Al,A2,A3,Ab,A5,A6)  giving 
instantaneous  charge  efficiency  as  a function  of  battery 
temperature,  normal ised  battery  charge  rate  and  battery 
state-of-charge 


Step  4l  Compare  Battery  Charge  Efficiency  With  Reference 
IF:  ETA  <0.0,  THEN:  ETA  = 0.00001 
Step  42  Calculate  Rate  of  Charge  of  Battery  State-of-Charge  With  Time 
WBDT  = ( BCUR  * ETA  )/CB 

Wliere : D<iBDT  * Rate  of  Charge  of  each  Battery  State-of-Charge  with 
Time  - Hours”! 

Step  44  Compare  Card  Type  with  Reference 

IF:  NCTYPE  * 0; 

THEN : H = 0;  AND, 

THEN:  GO  TO  STEP  49 


Step  43  Compare  Load  Selection  Indicator  With  Reference  and  Determine 
Reference  Calculation  Interval 


IF:  KL  « 1,  OR, 

IF:  KL  * 2, 

THEN:  GO  TO  STEP  49 

IF:  KLL  - 3, 

THEN:  HI  » H 

Where:  HI  » Performance  Analysis  Reference  Calculation  Interval  - 
Hours 

Step  46  Compare  Low  Insolation  Load  Selector  With  Reference  and  Calculate 
Performance  Analysis  Calculation  Interval 

IF:  KLL  » 1 

THEN:  GO  TO  STEP  47 

H ■ DL  • HI 

GO  TO  STEP  49 

Step  47  Calculate  Performance  Analysis  Calculation  Interval 
H ■ (1.0  - DL)  • HI 

Step  49  Calculate  Battery  State-of-Charge  Increment 
DEL4B  - H • DviBDT 

Where:  DEL4B  ■ State-of-Charge  Increment  for  each  Battery 
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Step  51  Compare  Card  Type  and  Low  Insolation  Load  Selector  With  References 

IF:  1ICTYPE  = 1,  AND. 

IF:  KLL  * 3, 

THEN:  GO  TO  STEP  57 

Step  53  Calculate  New  Battery  State-of-Charge 
•iB  = ♦ DELQB 

Step  Calculate  Battery  Electrolyte  Freezing  Temperature 

SPGH  * SPGR1  {QB} 

TBFRZ  * TBFRZ1  {SPGR} 

Where:  SPGR  * Electrolyte  Specific  Gravity  for  each  Battery 

TBFRZ  * Electrolyte  Freezing  Temperature  for  each  Battery  - ®F 

SPGH1  « Input  Table  of  SPGR  as  a function  of  QB 

TBFRZ1  * Input  Table  of  TBFRZ  as  a function  of  SPGR 

Step  Compare  Card  Type  With  Reference 

IF:  iJCTYPE  « 0 

THEN:  GO  TO  STEP  2 

GO  TO  STEP  63 

Step  57  Calculate  Battery  State-of-Charge  Relative  Increment 
D*B  ■ ABS(DEHlBAiB) 

Where:  DQB  ■ State-of-Charge  Relative  Increment  for  each  Battery 

Step  59  Compare  Battery  50C  Relative  Increment  With  Reference 

IF:  (H  • ACCvjB/OQB)  > (HIST  - h2  - 0.01) 

THEM:  Go  To  Step  60 

IF:  (OQB)  < (0.7  • ACCvJB),  OR: 

IF:  (DOB)  > (ACC4M) 

THEN:  H ■ (H  • ACC0B/DQB),  AND: 

THEN:  RETURN  TO  STEF  9 
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Step  bJ 


Stei  bl 


Step  62 


Step  62a 


Step  63 


Step  6k 


Step 


Step  fab 


Calculate  New  Mattery  State-of-Charge 
4M  = ul B ♦ DELQB 

Calculate  Mattery  Electrolyte  Freezing  Temperature 
SPGR  * SPGRIWm} 

TBFRZ  = TMFRZltSPGR} 

Compare  Load  Selection  Indicator  and  Low  Insolation  Load  Selection 
With  Reference  and  Calculate  Low  Insolation  Load  Selector 

IF:  KL  - 3,  AND. 

IF:  KLL  - 3, 

THEN:  KLL  ■ 1,  AND 
THEN:  CO  TO  STEP  12 

IF:  QB  - 0 

THEN:  Print  Error  Message,  AND 
THEN:  Return  to  Main  Program 

Calculate  Time  Reference  Data 

DATEM1  - DATUM 

DATE1  » DATE 

TIMEH  ■ 2h.O  • (DATE1  - IFIX(DATE)) 

YEAR1  ■ YEAR 

Calculate  Time  Elapsed  Since  Last  Print-Out 
H2  ■ H2  ♦ H 

Compare  Elapsed  Time  With  Duration  of  Time  Until  Next  Print-Out 
IF:  H2  > HINT 

THEN:  Print  Performance  Analysis  Output  Information 
THEN : CO  TO  STEP  66 

RETURN  TO  STEP  8 

Increment  Time  Step  Counter  and  Compare  With  Reference 

LNT3  ■ LOTS  ♦ 1 

IF:  LNTS  > NTS 

THK5:  00  TO  STEP  2 


RETURN  TO  STEP  7a 
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I4. 1 Power  Sources  Group 

The  Power  Source  Group  is  made  up  of  the  solar  array,  the  shunt  limiter, 
the  solar  array  isolation  diodes,  and  the  power  source  series  resistance. 
The  characteristics  of  these  elements  are  calculated  for  the  environ- 
mental conditions  in  which  the  subsystem  is  operated  and  are  then  com- 
bined into  a single  power  source  current-voltage  curve  at  the  unregulated 
bus.  Performance  data  are  stored  for  a single  solar  cell,  for  an  isola- 
tion diode,  for  the  specific  shunting  device  to  be  used,  and  for  the 
series  resistance  that  are  typical  of  those  in  the  buoy  solar  array. 

The  data  are  projected  from  the  component  level  into  the  electrical 
configuration  of  the  buoy  solar  array,  and  the  program  solves  sets  of 
equations  that  are  designed  to  >redict  solar  array/shunt  limiter  I-V 
characteristics  at  the  unregulated  bus.  The  Power  Source  Group  program 
also  includes  equations  to  estimate  the  array  performance  when  the  array 
is  misoriented  from  the  sun  vector  and  to  estimate  performance  degrada- 
tion due  to  cloud  cover,  temperature,  and  environmental  effects. 


PROGRAM  ALGORITHMS 

Step  1 

Obtain  the  exact  date  and  time  of  year 

TIKEH  ■ Daily  Time  - Hours  after  Midnight 

Range  of  values:  0-2L 

DATE  ■ Date  - Days  from  start  of  the  year 

Range  of  Values:  1-365 

Step  la 

IF:  ITAPE  i 0 

THEN:  Obtain  TTAMB  from  'MERGE'  file 

THEM:  GO  TO  STfcP  5 

Step  2 

Obtain  the  Average  Yearly  Temperature 

TTAVE  ■ Average  Yearly  temperature  in  selected  location  - *F 

W8 
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Step  j Calculate  Average  Daily  Temperature 

TTA  = TTAVE  + DTTA 

Where:  TTA  = Average  Daily  Temperature  - °F 

DTTA  ■ Average  Daily  Temperature  Increment  - °F 

DTTA  is  obtained  from  input  data  in  a table  of  average  daily 
tenperature  increment  as  a function  of  the  date,  ie., 

DTTA  ■ DTTA1{DATE) 

Step  1*  Calculate  Ambient  Temperature 

ITAMB  * TTA  ♦ DTTAM3 

Where:  TTAMB  ■ Ambient  Temperature  at  selected  location  - °F 

DTTAMB  ■ Average  hourly  temperature  increment  at 
selected  location  - °F 

DTTAMB  is  obtained  from  a table  of  Average  hourly  temperature 
increment  as  a function  of  the  time,  ie. , 

DTTAMB  ■ DTAMBl(TIMEH) 

Step  5 Obtain  Hover  Source  Croup  Equipment  Temperature  Characteristics 
DTTPSG  ■ Hover  Source  Group  Equipment  Temperature  Rise  - ®F 

Step  6 Calculate  Solar  Array  Temperature 

TSAF  ■ TTAMB  ♦ DTTPSG 

Where:  TSAF  ■ Solar  Array  Temperature  - °F 
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Step  7 Convert  Solar  Array  Temperature 

TSAR  = TSAF  + U59.67 
TSAK  = (5. 0/9.0)  # TSAR 
TSAC  = TSAK  - 273.15 

Where:  TSAR  = Solar  Array  Temperature  - °R 
TSAK  = Solar  Array  Temperature  - °K 
TSAC  = Solar  Array  Temperature  - °C 

Step  8 Calculate  Solar  Vector  Location  in  Equatorial  Plane 

AliPHEQ  * OMEGA  • DATE 

Where:  ALPHEQ  * Solar  Vector  Location  - Radians 
OMEGA  = (2  • 7T)/365.2U2 
tt  * 3. 11*159 

Note:  There  are  365.21*2  days  per  tropical  year  as  measured  from 
Vernal  Equinox  to  Vernal  Equinox 

Step  9 Calculate  Solar  Radiation  Variables 

VAR(I)  - FAO(I)  ♦ FA1(I)  • COS(ALPHEQ)  ♦ ... 

♦ FA2(I)  * COS(2.0  • ALPHEQ)  ♦ ... 

♦ FA3(I)  # C0S(3.0  • ALPHEQ)  ♦ ... 

♦ FBI ( I ) * SIN(ALPHEQ)  ♦ ... 

♦ FB2( I ) • SIN(2.0  • ALPHEQ)  ♦ ... 

♦ FB3(I)  • SIN( 3.0  • ALPHEQ) 
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Step  9 (contd) 

DECL  = VAR  (1)  • tt/180. 0 

ET  = VAR  (2) 

APPSC  = VAR  (3)  * 3. 152U808 

ATMEXC  = VAR  (1*) 

SDF  = VAR  (5) 

Where:  DECL  = Solar  Declination  Angle  - Radians 


ET  = Equation  of  Time  Difference  - Hours 

APPSC  * Apparent  Solar  Constant  - Watts/Meter2 
(at  AMO) 

ATMEXC  = Atomsphere  Extinction  Coefficient  - Air  Mass"1 
SDF  = Sky  Diffuse  Factor 

FA,  FB  = Fourier  Coefficients  obtained  from  input  data 
tables  "Solar  Radiation  Fourier  Coefficients" 


Step  10  Obtain  Buoy  Latitude 

THELAD  « Buoy  latitude  - degrees  ( + North 

I - South 


Step  11  Convert  Buoy  Latitude 

THETLA  - THELAD  • n/l80.0 

Where:  THETLA  ■ Buoy  Latitude  - Radians 

Step  12  Calculate  Terminator  Hour  Angle 

IF:  THETLA  > [(n/2.0)  - DECL] 

THEM:  HCURT  ■”* 

Go  to  Step  13 

HOURT  » ARCCOS  (-1,0  • TAM  (THETLA)  • TAM  (DECL)) 


Where:  HOURT  ■ Terminator  Hour  Angle  - Radians 


Convert  Terminator  Hour  Angle 
HOURA  = HOURT  * 12.0/tt 


Where:  HOURA  = Terminator  Hour  Angle  - Hours 

Step  ll*  Obtain  Buoy  Location  Time  Zone  Number 

TZN  = Time  Zone  Number  (Hours  behind  Greenwich  Mean  Time) 


Obtain  Buoy  Longitude 


THELQD  = Buoy  Longitude  - degrees 


♦ West 
- East 


Step  16  Calculate  Time  of  Sunrise  and  Sunset  at  Buoy  Location 
SRT  = 12.0  - HOURA  - ET  - TZN  ♦ (THEL0D/15.0) 

SST  = 21*. 0 - SRT 

Where:  SRT  * Sunrise  time  - Hours 
SST  * Sunset  time  - Hours 

Step  IT  Calculate  Buoy  Location  Hour  Angle 

BHOURD  « 15.0  * (TIMEH  - 12.0  ♦ TZN  ♦ ET)  - THELOD 
BHOUR  « BHOURD  • ir/180.0 

Where:  BHOURD  ■ Buoy  Location  Hour  Angle  - Degrees 
BHOUR  * Buoy  Location  Hour  Angle  - Radians 

Step  18  Test  for  Solar  Occultation 

IF:  ABS  (BHOUR)  _>  ABS  (HOURT) 

Go  to  Step  118 

Step  19  Calculate  Direction  Cosines  of  Direct  Solar  Radiation 

COS  (THETZS ) ■ COS  (BHOUR)  • COS  (DECL)  • COS  (THETLA)  ♦ ... 
♦ SIR  (DECL)  • SIR  (THETLA) 
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Step  19  (contd) 

Where:  THETZS  = Angle  between  the  local  zenith  and  the  solar 
vector  - radians 

COS  (THW)  = COS  (DECL)  * SIN  (BHOUR) 

IF:  COS  (BHOUR)  > 

THEN:  KS  = 1.0 

cos  < |Ty (ffli  i 

THEN:  KS  = -1.0 

I 9 9 ) 0 • 5 

COS  (THS)  = KS  * [l  - COS  (THETZS)]^  - [COS  (THW)]  ) 

Where:  THW,  THS  = Additional  Direction  Angles  - Radians 


Step  ?0  Calculate  Solar  Altitude 

SALT  = ARCSIN  (COS  (THETZS)) 

Where:  SALT  = Solar  Altitude  (Angle  between  the  solar  vector  and 
the  Horizontal,  i.e..  Earth's  surface)  - Radians 


Step  21  Calculate  Solar  Azimuth 
IF:  COS  (THS)  > 0 

THEN:  SAZM  - ARCSIN  [COS  (THW) /COS  (SALT)] 

CO  TO  STEP  22 

IF:  COS  (THS)  <0 

THEN:  SAZM  ■ tr  - ARCSIN  [COS  (THW)/COS  (SALT)] 

Where:  SAZM  * Solar  Azimuth  (Angle  between  the  Solar  Vector. 

Projected  onto  the  Horizontal  Surface  and  the 
South-Pointing  Vector  on  the  Horizontal 
Surface)  - Radians 
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Step  2?  Obtain  Cloud  Cover  Conditions 
CT  = Cloud  Type 

0.0  = Cirrus  or  Cirrostratus  Clouds 

1.0  = Stratus  Clouds 

2.0  = Other  Cloud  Types 
TC  * Total  cloud  Cover 

1.0  = 1/10  of  sky  covered 

2.0  = 2/10  of  sky  covered 
• 

• 

9.0  * 9/10  of  sky  covered 

10.0  * 10/10  of  sky  covered 
ICT  * 1 ♦ IFIX  (CT) 

Where:  ICT  * Cloud  Type  indicator 

1 * Cirrus  of  Cirrostratus  Clouds 

2 » Stratus  Clouds 

3 ■ Other  Cloud  Types 

Step  23  Calculate  Cloud  Cover  Modifier 
IF:  TC  « 0.0, 

THEW:  CCM  » 1.0  and  go  to  Step  2U 

IF:  SALT  < wA.O; 

THEW:  I SALT- ■ 1 

IF:  SALT  > "A.0; 

THEW:  I SALT  ■ 2 

Where:  ISALT  ■ Solar  Altitude  Indicator 

CCM  - PO  (ICT, ISALT)  ♦ PI  (ICT, ISALT)  • TC  ♦ ... 

♦ P2  ( ICT, ISALT)  • (TC*»2.0)  ♦ ... 

♦ P3  (ICT, ISALT)  • (TC»»3.0) 
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Step  23  (contd) 

Where:  CCM  = Cloud  Cover  Modifier 

PO,  PI,  P2,  P3  = Polynomial  Coefficients  obtained  from 

input  data  tables  "Cloud  Cover  Modifier 
Polynomial  Coefficients" 

Step  2k  Obtain  Clearness  Number 
CN  * Clearness  Number 

« 0. 7*0.9  for  an  industrial  atmosphere 
= 0.85-»..0  for  non-industrial  atmospheres 

Step  23  Calculate  Intensity  of  Direct  Normal  Solar  Radiation 
QDN  * APPSC  • CN  • CCM  • EXP  (-ATMEXC/COS  (THETZS)) 

Where:  QDN  = Direct  Normal  Solar  Radiation  Intensity  - Watts/Meter2 

Step  26  Obtain  Solar  Array  Pointing  Angles 

PHIAID  « Surface  Tilt  Angle  from  Horizontal  - Degrees 

(Angle  between  local  Zenith  and  Solar  Array  Normal) 

PHIAAD  ■ Surface  Azimuth  Angle  from  South  - Degrees 

(Angle  between  South  pointing  vector  and  projection 
of  array  normal  on  horizontal  surface) 

■ + if  West  of  South 

■ - if  East  of  South 

Step  27  Convert  Solar  Array  Pointing  Angles 
PHIAI  ■ PHIAID  • s/180.0 
PHIAA  ■ PHIAAD  • ir/100.0 
Where:  PHIAI  ■ Surface  Tilt  Angle  - Radians 

PHIAA  ■ Surface  Azimuth  Angle  - Radians 
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Step  28 


Step  ?9 


Step  30 


Step  jl 


Step  32 


Calculate  Direction  Cosines  of  Array  Normal 

(Reference  Axis:  Vertical,  Horizontal  to  West,  Horizontal  to  South) 
ETAA  = COS  (PHIA1 ) 

ETAB  = SIN  (PHIAA)  * SIN  (PHIAI) 

ETAC  = COS  (PHIAA)  * SIN  (PHIAI) 

Where:  ETAA,  ETAB,  ETAC  = Array  Normal  Direction  Cosines 


Calculate  Solar  Array  Tilt  Angle 

COO  (TILT)  = ETAA  * COS  (THETZS)  + ... 

+ ETAB  * COS  (THW)  + ... 

+ ETAC  * COS  (THS ) 

Where:  TILT  = Solar  Array  Tilt  Angle  - Radians 

(Angle  between  Solar  Vector  and  Solar  Array  Normal) 


Calculate  Intensity  of  Direct  Solar  Radiation  Incident  on  the 
Solar  Array 

IF:  COS  (TILT)  > 0.0 

THEN:  QD  = QDN  • COS  (TILT) 

IF:  COS  (TILT)  < 0.0 

THEN:  QD  = 0.0 

Where:  QD  * Direct  Solar  Radiation  Incident  on  Solar 
Array  - Watts/Meter2 


Calculate  Sky  Brightness 
BS  * SUF  • QDN/ (CN**2.0 ) 

p 

Where:  BS  * Sky  Brightness  - Watts/Meter 

Obtain  Horizontal  Surface  (Ground/Ocean)  Reflectivity 
REFLH  » Horizontal  Surface  Reflectivity  for  Solar  Radiation 
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Step  33  Calculate  Horizontal  Surface  Brightness 
BG  = REFLH  * (BS  + QDN  * COS  (THETZS)) 

2 

Where:  BG  = Horizontal  Surface  Brightness  - Watts/Meter 


Step  3j*  Calculate  Intensity  of  Horizontal  Surface  Diffuse  Radiation 
Incident  on  Solar  Array 

QDG  = BG  • ((1  - ETAA )/2.0) 

Where:  QDG  - Horizontal  Surface  Diffuse  Radiation  Incident 
on  Solar  Array  - Watts/Meter2 


Step  35  Calculate  Intensity  of  Sky  Diffuse  Radiation  Incident  on  a 
Horizontal  Solar  Array 

QIJSH  = QDN  * SDF 

Where:  QDSH  * Sky  Diffuse  Radiation  Incident  on  a Horizontal 
Solar  Array  - Watts/Meter2 


Step  jC  Calculate  Intensity  of  Sky  Diffuse  Radiation  Incident  on  a 
Vertical  Solar  Array 

YV  = O.I45 

IF:  COS  (TILT)  > (-0.2); 

THEM:  YV  « 0.55  ♦ O.J4 37  * COS  (TILT)  ♦ 0.313  * ((COS  (TILT)»«2.0) 

QDSV  « QDN  * (SDF  » YV  ♦ (REFLH  • (SDF  ♦ COS  (THETZS) ) )/2.0) 

Where:  QDSV  * Sky  Diffuse  Radiation  Incident  on  a Vertical 
Solar  Array  - Watts/Meter2 


Step  37  Calculate  Intensity  of  Sky  Diffuse  Radiation  Incident  on  Solar 
Array 

QDG  « QDGV  ♦ (QDSH  - QDSV)  • COS  (SALT) 

Where:  QDS  ■ Sky  Diffuse  Radiation  Incident  on  Solar 
Array  - Watts/Meter2 
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Step  38  Calculate  Intensity  of  Total  Solar  Insolation  Incident  on  Solar 
Array 

QDT  = QD  + QDG  + QDS 

Where:  QDT  = Total  Solar  Radiation  Incident  on  Solar 
Array  - Watts/Meter*5 


Step  39  Obtain  Elapsed  Time  From  Start  of  Mission 

DATUM  = Elapsed  time  from  start  of  mission  - days 


Step  ^0  Obtain  Current  Degradation  Factors  for  Solar  Array 

CDFGA  = Solar  Array  Current  Degradation  Factor  Due  to 
Fabrication  Losses  - Percent  (from  zero) 

'DEC®  = Solar  Array  Current  Degradation  Factor  Due  to 
Terrestrial  Performance  Extrapolation 
Uncertainty  - Percent  (from  zero) 


Step  Li  Calculate  Current  Degradation  Factor  Due  to  Environmental  Effects 
ESC  = SADEGC  (DATEM) 

Where:  CDEGC  * Solar  Array  Current  Degradation  Factor  Due  to 

Environmental  Effects  - Percent  (from  zero) 

SADEGC  * Table  of  Solar  Array  Input  Current  Degradation 
Due  to  the  Environment  (in  Percent  from  Zero) 
as  a function  of  DATEM 


Step  h<?  Calculate  Solar  Array  Current  Degradation  Factor 

CDEC  . 1.0  * I0b  - (100.0  - CDEGA)  » (100.0  - CDEGB)  • (100. U - CDHCC) 

1.0  • KT 

Where:  CDEG  ■ Solar  Array  Current  Degradation 
Factor  - Dimensionless 
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Step  I43  Obtain  Voltage  Degradation  Factor  for  Solar  Array 

VDEGA  = Solar  Array  open  circuit  voltage  degradation  due  to 
temperature  uncertainty  - Percent  (from  zero) 


Step  I4I4  Calculate  Voltage  Degradation  Factor  Due  to  Environmental  Effects 
VDEGB  = SADEGV  (DATEM) 

Where:  VDEGB  = Solar  Array  0pe:n  Circuit  Voltage  Degradation 

Factor  due  to  Environmental  Effects  - 
Percent  (from  zero) 

SADEGV  = Table  of  Solar  Array  Open  Circuit  Voltage 

Degradation  due  to  the  Environment  (in  percent 
from  zero)  as  a function  of  DATEM 


Step  I45  Calculate  Solar  Array  Voltage  Degradation  Factor 

VDEG  * 1-0  * - (100.0  - VDEGA)  » (100.0  - VDEGB) 

1.0  * V 

Where:  VDEG  * Solar  Array  Voltage  Degradation 
Factor  - Dimensionless 


Step  I46  Obtain  Solar  Cell  Spectral  Correction  Factor 

SPEC0R  ■ Solar  Cell  Spectral  Correction  Factor  - Dimensionless 
(Corrects  for  differences  between  Spectrum  of  Solar 
Radiation  Incident  on  Solar  Cell  and  Spectral  Response 
of  Solar  Cell) 


Step  h7  Calculate  Effective  Solar  Insolation 
X ■ 5PEC0R  • QDT/10.0 

Where:  X ■ Effective  Solar  Insolation  Incident  on  Solar 
Cell  - Millivatts/Csr 
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Step  ^8  Calculate  Modified  Solar  Insolation 
XX  = X * (1.0  - CDEG) 

Where:  XX  - Modified  Solar  Insolation  - Mw/Cnf 

Step  1*9  Obtain  Single  Solar  Cell  Area 

ACELL  = Single  Solar  Cell  Area  - Cm' 

Step  V)  Calculate  Short  Circuit  Current  Temperature  Coefficient  for  a 
Single  Solar  Cell 

AI.PHAC  = ((7.1*28  • 10‘7)  - (1.83  * 10'9)  * TSAc)  • (XX)  * ACELLA.O 

Where:  ALPHAC  * Short  Circuit  Current  Temperature 
Coefficient  - Amperes/°C-cell 

Step  51  Calculate  Solar  Cell  Series  Resistance 

RCELLC  = FtRSCELL,  TEMTAB}  at  TSAC 

Where:  RCELLC  ■ Solar  Cell  Series  Resistance  - ohms 
(at  Temperature  TSAC) 

RSCELL  » Internal  Table  of  Solar  Cell  Series  Resistance 
as  a function  of  Cell  Temperature 

TEMTAB  • Internal  Table  of  Temperature  Range  Associated 
with  RSCELL 

Step  51A  Calculate  Solar  Cell  I -V  Curve  Correction  Factor 
ROCELL  « F (ROE,  SUNLIT)  at  XX 

Where:  ROCELL  ■ Solar  Cell  I-V  Curve  Correction  Factor  at 
Solar  Iniolation  Level:  XX 

ROE  ■ Internal  Table  of  Solar  Cell  I-V  Curve  Correction 
Factor  as  a Function  of  Solar  Insolation 

SUNLIT  ■ Internal  Table  of  Solar  Insolation  Range  Associated 
with  ROE 
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Step  5."  Calculate  Open  Circuit  Voltage  Temperature  Coefficient  for  a 
Single  Solar  Cell 

BETAA  = F[BETAB(or  BETAC  or  BETAD)]  at  XX  and  TSAC 
BSETA  = BKTAA/ 1000.0 

Where:  BBETA  * Open  Circuit  Voltage  Temperature 

Coefficient  - (Volts/°C)  at  XX 
and  TSAC 

BETAA  = Open  Circuit  Voltage  Temperature 
Coefficient  - (Mv/°C)  at  XX  and 
TSAC 

BETAB,  BETAC,  BETAD  = Internal  Tables  of  Solar  Cell  Open 

Circuit  Voltage  as  a Function  of 
Solar  Insolation  and  Cell 
Temperature 

SUNMW,  S0NMW,  SENMW  «=  Internal  Tables  of  Solar  Insolation 

Ranges  Associated  with  (BETA)  Tables 

BTEMP,  CTEMP,  DTEMF  « Internal  Tables  of  Solar  Cell 

Temperature  Ranges  Associated  with 
(BETA.,)  Tables 


Internal  Tables  BTEMP,  SUNMW  AND  BETAB  used  when: 

(100  < XX  <,  51*0  Mw/Cm2)  and  (-60  < TSAC  < l60°C) 
Internal  Tables  CTEMP,  SONMW,  BETAC  used  when: 

(5  1 XX  < 253  Mw/Co2)  and  (-1*0  < TSAC  < 6o°C) 
Internal  Tables  DTEMP,  SENMW,  BETAD  used  when: 

(5  £ XX  < 100  Mw/Cb2)  and  (-11*0  < TSAC  <_  -1*0°C) 

\ 

Stet  53  Obtain  Single  Cell  ISC,  V0C  Data  l 

IISC  ■ Solar  Cell  Short  Circuit  Current  - Aaperes/Cell 
(at  ll«5  Mw/Cb^  Solar  Insolation  and  60°C) 

WOC  ■ Solar  Cell  Open  Circuit  Voltage  - Volts/Cell 
(at  1*5  Mw/c£*2  solar  Insolation  and  60*C) 
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Calculate  ICC,  VOC  Shift  Due  to  Degradation 


Step  35 


Step  56 
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Cl  = CDEG  * IISC 

C2  = VDEG  * VVOC 

Cl  = Solar  Cell  Short  Circuit  Current  Shift  - Amps/Cell 
C2  = Solar  Cell  Open  Circuit  Voltage  Shift  - Volta/Cell 

Obtain  Single  Circuit  (of  Solar  Cells)  Arrangement 
HS  = Ho.  of  Solar  Cells  in  Series  in  Each  Circuit 

NP  * Ho.  of  Solar  Cells  in  Parallel  in  Each  Circuit 

Calculate  Cell  Electrical  Circuit  Parameters 
ALPHA  = ALPHAC  * HP 
BETA  = BBETA  * HS 
RCELL  * (O.lll  ♦ BCELLC ) * NS/HP 
HHO  * ROCKLL  * HS/HP 

Where:  ALPHA  * Short  Circuit  Current  Temperature  Coefficient 
for  a Single  Circuit  - Anpcres/°C-circuit 

BETA  » Open  Circuit  Voltage  Temperature* Coefficient 
for  a Single  Circuit  - Volts/°C 

RCELL  » Single  Circuit  Series  Resistance  - Ohms 

RHO  » Series  Resistance  Temperature  Correction  Factor 

Calculate  Modified  Electrical  Circuit  Short  Circuit  Current 
ISC  - IISC  • BP  * (1.0  - CDEG) 


Where:  ISC  ■ Modified  Electrical  Circuit  Short  Circuit 
Current  - Ampere*/ Circuit 


Calculate  Short  Circuit  Current  Differences  (for  an  Electrical 
Circuit ) 

DISC  = ISC  * ( (X/1U5.0 ) - 1.0)  ♦ ALPHA  * (TSAC  - 60.0) 


Where:  DISC  = Short  Circuit  Current  Difference  due  to  current 

degradation,  solar  insolation  changes  and  tempera- 
ture changes  - Amperes/Circuit 


Step  39  Calculate  Electrical  Circuit  Voltage  and  Series  Resistance 
Correction  Factors 

C3  = BETA  * (TSAC  - 60.0)  ♦ DISC  • RCELL 
Cl*  = RHO  * (TSAC  - 60.0) 

Where:  C3  * Electrical  Circuit  Voltage  Correction 
Factor  - Volts/Circuit 

Cl*  * Electrical  Circuit  Series  Resistance  Correction 
Factor  - Ohms 

Step  60  Obtain  Reference  Solar  Cell  Current-Voltage  Characteristics 

II ( J ) * Reference  Solar  Cell  Current  Data  Point  - Amperes 
(Internal  Tables) 

VV(J)  * Reference  Solar  Cell  Voltage  Data  Point  - Volts 
(Internal  Tables) 

J « 1,  30 

Step  6l  Calculate  Solar  Cell  Electrical  Circuit  Current-Voltage 
rharact*;ristic8 

I ( J ) r HP  • (II(J)  - Cl)  ♦ DISC 

V(J)  ■ NS  * (W(J)  - C2)  - C3  - (Cl*  • I(J)) 

J ■ 1,  30 

Where:  I ( J ) ■ Electrical  Circuit  Current  - Amperes 
at  the  given  level  of  V(J) 

V(J)  ■ Electrical  Circuit  Voltage  - Volts 
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Step  62  Obtain  Solar  Array  Voltage  Increment 


Step  63 


Step  6U 


Step  65 


VSAINC  = Solar  Array  Voltage  Increment  - Volts 


Redefine  Electrical  Circuit  Current-Voltage  Array  in  Selected 
Voltage  Increments  as  follows: 

a)  Set:  Counter  L=1  and  voltage  V2(L)  = 0.0 

b)  Establish:  Current  I1(L)  at  V2(L) 

II ( L ) = F{l(J) , V(J)}  at  I ( J ) = 0.0 

c)  Increment:  Counter  L 3 L + 1 and  voltage  V2(L  ♦ 1)  = V2(L) 

+ VSAINC  and  Establish:  Current  II  (L  + l)  at 
V2(L  + 1).  Until:  I1(L  ♦ l)  < 0.0 

d)  Redefine:  Last  V2(L)  at  I1(L)  3 0.0 
V2(L)  = F{I(J),  V(J)}  at  I(J)  3 0.0 

e)  Set:  Current-Voltage  Matrix  Dimension  to  last  counter 

value  MFINAL  = L 


Obtain  Number  of  Solar  Cell  Electrical  Circuits  in  Solar  Array 

NESP  3 Number  of  Electrical  Circuits  in  Solar  Array  (assumed 
in  parallel) 


Calculate  Solar  Array  Current-Voltage  Characteristics 
I2(L)  ■ I1(L)  * NESP  at  V2(L) 

L * 1,  MFINAL 

Where:  Il(L)  ■ Electrical  Circuit  Current  - Amperes 
at  V2(L) 

I2(L)  3 Solar  Array  Current  - Amperes 
at  V?(L) 

V2(L)  ■ Circuit  or  Array  Voltage  - Volts 


Step  66  Obtain  Voltage  Data  for  Calculation  of  Solar  Array  Maximum  Power 
Plant 

XV  » Initial  Voltage  for  Max  Power  Point  Calculations  - Volts 
DXN  3 Voltage  increment  for  Max  Power  Point  Calculation  - Volts 


B mmm 


Initialize  Calculation  Value  of  Solar  Array  Maximum  Power 
MSAPWR  =0.0 


Where:  MSAPWR  = Solar  Array  Maximum  Power  - Watts 

Step  68  Calculate  Solar  Array  Power  and  Current 
XI  = F{I2(L),  V2(L)}  at  V2(L)  = XV 
SAPWR  = XI  * XV 

Where:  XV  = Solar  Array  Voltage  - Volts 

XI  * Solar  Array  Current  - Amperes 
SAPWR  = Solar  Array  Power  - Watts 

Step  A9  Compare  Solar  Array  Power  With  Maximum  Power 

IF:  SAPWR  > MSAPWR 

I MSAPWR  = SAPWR 
I XV  * XV  + DXV 

REPEAT  STEP  68  UNTIL:  SAPWR  <_  MSAPWR 

Step  70  Recalculate  Solar  Array  Current  and  Power 
MSAPWR  * 0.0 

XV  = XV  - DXV 

REPEAT  STEP  68  ONLY 

Step  71  Compare  Solar  A. ray  Power  With  Maximum  Power 

IF:  SAPWR  > MSAPWR 

MSAPWR”-  SAPWR  j 
DXV  * DXV/lO.Oj 
XV  - XV  ♦ DXV ) 

REPEAT  STEP  68  ONLY  UNTIL:  SAPWR  < MSAPWR 
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Step  72  Calculate  Solar  Array  Maximum  Power  Point  Characteristics 
MAXV  = XV  - DXV 
MAXI  = M5APWR/MAXV 

Where:  MAXV  = Solar  Array  Voltage  at  Max  Power  Point  - Volts 

MAXI  = Solar  Array  Current  at  Max  Power  Point  - Amperes 

Step  73  Obtain  Solar  Array  Electrical  Section  Cable  Resistance 

RSA  = Series  resistance  of  cable  for  an  electrical  section  of  the 
solar  array  - Ohms 

Step  7I4  Calculate  Voltage  Shift  in  Electrical  Section  Voltage  Due  to  Cable 
Resistance  and  Blocking  Diodes 

VDIODE  = ADI  { I 1 ( L ) } 

V2(L)  = V2(L)  - ( I 1 ( L ) * RSA)  - VDIODE 

L * 1,  MFIUAL 

Where:  VDIODE  » Electrical  Section  Blocking  Diode  - Volts 
Voltage  Drop  at  Current  Level  I1(L) 

ADI  * Table  (Input  Data)  of  Electrical  Section  Blocking 
Diode  as  a function  of  current 

Step  75  Obtain  Shunt  Limiter  Type 
ISH  * Shunt  Limiter  Type 

0 * Hone 

1 * Ordinary  Zener  Diode 

2 * Temperature  Compensated  Zener  Diode 

3 ■ Active  Shunt  Limiter 

Step  76  Obtain  Shunt  Limiter  Current  Voltage  Characteristics 
ZI(I)  * Shunt  Limiter  Current  at  ZV(I)  - Amperes 
ZV(I)  * Shunt  Limiter  Voltage  - Volts 
I « 1,  20 
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Step  77  Obtain  Power  Source  Group  Type 
IPSG  = Power  Source  Group  Type 

0 = One  Shunt  Limiter  for  the  Solar  Array 

1 = One  Shunt  Limiter  for  each  Electrical  Section  of  the 

Solar  Array 

Step  78  Select  Power  Source  Group  Current-Voltage  Calculation  and 
Calculate  Significant  Voltages 

VZS3  = ZV(2) 

SAOCV  = V2(MFINAL) 


IF: 

ISH  = 

0; 

(iO 

TO 

STEP 

79 

IF: 

ISH  = 

1; 

GO 

TO 

STEP 

80 

IF: 

ISH  = 

2; 

GO 

TO 

STEP 

98 

IF: 

ISH  = 

3; 

GO 

TO 

STEP 

80 

Where:  VZSB  = Shunt-Limiter  Turn-On  Voltage  - Volts 

SAOCV  = Solar  Array  Open  Circuit  Voltage  - Volts 

Step  79  Calculate  Power  Source  Group  Current-Voltage  Characteristics 
(With  No  Shunt-Limiter) 

XZI(L)  = ZI(1) 

SY(L)  - I2(L)  - XZI(L) 

SX ( L ) « V2(L) 

L ■ 1,  MFINAL 

Where:  XZI(L)  ■ Shunt  Limiter  Current  at  SX(L)  - Amperes 

SY(L)  ■ Power  Source  Group  Current  at  SY(L)  - Amperes 
SX(L)  ■ Power  Source  Group  Voltage  - Volts 
Return  to  Performance  Analysis  Routine 
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.'top  80 


Sole 

?ct  PSG 

i Current 

-V 

1 1 age 

Calculation  Based  on  PSG  Type 

_I_F: 

I PSG 

= 0; 

CO 

TO 

STEP 

81 

IF: 

I PSG 

* l; 

CO 

TO 

STEP 

89 

Step  6l  Initialize  Index  Counter  and  PSG  Voltage  and  Current 


LL  = 1 


sy(i)  = o.o 

GX(l)  = F{I2(L),  V2(L) } at  VP(L)  = SX(l) 

XZ 1(1)  = 0.0 


Step  82  Compare  Voltage  and  Current  to  Reference  Levels 

IJF:  SY(LL)  < 0.0:  and 

I_F:  CX(LL)  > SAOCV 

THEN:  CO  TO  STEP  88 

Step  83  Increment  Index  Counter  and  Calculate  PSG  Voltage  and  Current 
LL  = LL  ♦ 1 

SX(LL)  = SX(LL  - l)  ♦ VSAINC 

SY(LL)  * F{I2(L),  V2(L)}  at  V2(L)  = SX(LL) 

XZI(LL)  « 0.0 

Step  84 


Compare  PSG  Voltage  with  Shunt-Limiter  Turn-On  Voltage 

IF:  SX(LL)  < VZSB: 

THEN:  REPEAT  STEPS  82  AND  83 
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Step  85 


Step  86 


Step  87 


Step  88 


Calculate  PSG  Current  at  Shunt  Limiter  Turn-On 
SX(LL)  = V2SB 

SY(LL)  = F{I2(L),  V?(L)}  at  V2(L)  = SX(LL) 
XZI(LL)  = F{ZI (I ) , ZV(I )}  at  ZV(I)  = SX(LL) 
SY(LL)  = SY(LL)  - XZI(LL) 


Compare  PSG  Current  With  Reference 

IF:  SY(LL)  < 0.0 

THEN : GO  TO  STEP  88 

Increment  Index  Counter  and  Calculate  PSG  Voltage  and  Current 
LL  = LL  + 1 

SX(LL)  = SX(LL  - 1)  + 0.01 

3Y(LL)  = F{I2(L),  V2(L)}  at  V2(L)  « SX(LL) 

XZI(LL)  = F{ZI( I ) , ZV(I)}  at  ZV(l)  « SX(LL) 

SY(LL)  = SY(LL)  - XZI(LL) 

REPEAT  STEPS  86  AND  87 

Calculate  Maximum  PSG  Voltage 

Perform  a Straight-Line  Interpolation  between  the  last  two  sets 
of  PSG  current-voltage  vaive3  to  predict  the  voltage  at  which 
the  PSG  current  is  equal  to  zero  (VZCR) 

SX(LL)  ■ VZCR 

SY(LL)  ■ 0.0 

XZI(LL)  ■ F{ZI(I),  ZV(I))  at  ZV(I)  ■ SX(LL) 

NFINAL  - LL 

RETURN  TO  PERFORMANCE  ANALYSIS  ROUTINE 
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Step  93  Calculate  Section  Current  at  Shunt  Limiter  Turn-On 
VSECT(LL)  = VZSB 

ISECT(LL)  = F{I1(L),  V2(L)}  at  V2(L)  = VSECT(LL) 

XZIS(LL)  = F{ZI(I),  ZV( I )}  at  ZV(l)  = VSECT(LL) 

ISECT(LL)  = ISECT(LL)  - XZIS(LL) 

Step  9l*  Compare  Section  Current  With  Reference 

IF:  ISECT(LL)  <0.0 

THEN:  GO  TO  STEP~96 

Step  95  Increment  Index  Counter  and  Calculate  Section  Voltage  and  Current 
LL  * LL  ♦ 1 

VSECT(LL)  *=  VSECT(LL  - l)  ♦ 0.01 

ISECT(LL)  * F{I1(L),  V2(L)}  at  V2(L)  =*  VSECT(LL) 

XZIS(LL)  = F{ ZI(I ) , ZV(I)}  at  ZV(I)  - VSECT(LL) 

ISECT(LL)  * ISECT(LL)  - XZIS(LL) 

REPEAT  STEPS  94  AND  95 

Step  96  Calculate  Maximum  Section  Voltage 

Perform  a straight-line  interpolation  between  the  last  two  sets 
of  section  current-voltage  values  to  predict  the  voltage  at  which 
the  section  current  is  equal  to  zero  (VZCR) 

VSECT(LL)  ■ VZCR 

ISECT(LL)  ■ 0.0 

XZIS(LL)  ■ F{ZI(I)»  ZV(I)}  at  ZV(I)  ■ VSECT(LL) 
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NFINAL  - LL 


Step  97  Calculate  PSG  Current  Voltage  Characteristics 
SY(LL)  = ISECT(LL)  * NESP 
XZI(LL)  = XZIS(LL)  * NESP 
SX(LL)  = VSECT(LL) 

LL  = 1,  NFINAL 

RETURN  TO  PERFORMANCE  ANALYSIS  ROUTINE 

Step  98  Select  PSG  Current  Voltage  Type  Based  on  PSG  Type 
.IF:  I PSG  = 0;  GO  TO  STEP  99 
IF:  I PSG  = 1;  GO  TO  STEP  108 

Step  99  Initialize  Index  Counter  and  PSG  Voltage  and  Current 
LL  * 1 

SX(1)  * 0.0 

SY(l)  * F{ I?(L) , V2(L)}  at  V2(L)  « SX(l) 

XZI(l)  * 0.0 

Step  100  Compare  Voltage  and  Current  to  Reference  Levels 

IF:  SY(LL)  < 0.0;  and 

IF:  SX(LL)  > SAOCV: 

THEM:  go  TO  STEP  107 

Step  101  Increment  Index  Counter  and  Calculate  PSG  Voltage  and  Current 
LL  ■ LL  ♦ 1 

SX(LL)  » SX(LL  - 1)  ♦ VSAIHC 

SY(LL)  ■ F{l?(L),  V2(L)}  at  V?(L)  ■ SX(LL) 

XZI(LL)  - 0.0 
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Step  102  Compare  PSG  Voltage  With  Shunt-Limiter  Turn-On  Voltage 

IF:  SX(LL)  < VZSB; 

THEN:  REPEAT  STEPS  100  AND  101 

Step  lOj  Calculate  PSG  Current  at  Shunt-Limiter  Turn-On 
SX(LL)  = VZSB 

SY(LL)  = F{I2{L),  V2(L)}  at  V2(L)  = SX(LL) 

XZI(LL)  = F{ZI(I),  ZV(I)}  at  ZV(I)  = SX(LL) 

SY(LL)  = SY(LL)  - XZI<LL) 

Step  I OU  Compare  PSG  Current  with  Reference 

IF:  SY(LL)  < 0.0 

THEN:  GO  TO  STEP  107 

Step  1 Calculate  Zener  String  Voltage  Increment 

VZINC  = ZV(1»)  - ZV(3) 

VZINC  ■ Zener  String  Voltage  Increment  - Volts 
beyond  Zener  Breakdown  Voltage 

Stop  100  Increment  Index  Counter  and  Calculate  PSG  Voltage  and  Current 
LL  * LL  ♦ 1 

SX(LL)  - SX( LL  - 1)  ♦ VZiKC 

SY(LL)  - F{I2(L),  V2(L)}  at  V?(L)  ■ SX(LL) 

XZI(LL)  ■ F{ZI(I ),  ZV( I )>  at  ZV(I)  - SX(LL) 

SY(LL)  - SY(LL)  - XZI(LL) 

REPEAT  STEPS  104,  105  AND  106 
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Step  1 37 


Calculate  Maximum  PSG  Voltage 


Perform  a straight-line  interpolation  between  the  lust  two  sets 
of  PSG  Current-Voltage  values  to  predict  the  voltage  at  which  the 
PSG  current  is  equal  to  zero  (VZCR) 

SX(LL)  = VZCR 

SY(LL)  = 0.0 

XZI(LL)  = F{ZI(I),  ZV( I ) } at  ZV( I ) = SX(LL) 

NFINAL  = LL 

RETURN  TO  PERFORMANCE  ANALYSIS  ROUTINE 


Step  108  Initialize  Index  Counter  and  Section  Voltage  and  Current 


LL  = 1 


VSECT(i)  = 0.0 

ISECT(l)  * Hll(L),  VP ( L ) } at  V2(L)  = VSECT(l) 
XZIS(l)  = 0.0 


Step  10 5 


Compare  Voltage  and  Current  to  Reference  Levels 


IF:  ISECT(LL)  <_  0.0;  and 

JF:  VSECT(LL)  > SA0CV; 

THEN:  GO  TO  STEP  116 


Step  110 


Increment  Index  Counter  and  Calculate  Section  Composite  Voltage 
and  Current 

LL  * LL  ♦ 1 

VSECT(LL)  ■ VSECT( LL  - 1)  ♦ VGA INC 

ISECTUL)  ■ F{ I1(L),  V?{L) ) at  VP(L)  * VSKCT(LL) 

XZIS(LL)  - 0.0 
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Step  111  Compare  Section  Voltage  With  Shunt  Limiter  Turn-On  Voltage 

IF:  VSECT(LL)  < VZSB; 

THEN:  REPEAT  STEPS  109  AMD  110 

Step  112  Calculate  Section  Current  at  Shunt  Limter  Turn-On 
VSECT( LL)  = VZSB 

I SECT ( LL)  * F{I1(L),  V2(L)}  at  V2(L)  = VSECT(LL) 

XZIS(LL)  = F{ZI( I ) , ZV( I )>  at  ZV( I ) = VSECT(LL) 

ISECT(LL)  * ISECT(LL)  - XZIS(LL) 

Step  113  Compare  Section  Current  With  Reference 

IF:  ISECT(LL)  <.  0.0 

THEN:  GO  TO  STEP”  116 

St vC  114  Calculate  Zener  String  Voltage  Increment 
VZINC  - ZV(4 ) - ZV( 3) 

Step  115  Increment  Index  Counter  and  Calculate  Section  Voltage  and  Current 
LL  ■ LL  ♦ 1 

VSECT(LL)  «■  VSECTtLL  - l)  ♦ VZINC 
ISECT(LL)  • F(I1(L),  V2(L))  at  V2{L)  * VSKCT(LL) 

XZIS(LL)  » F{ZI( I ),  ZV(I))  at  ZV{ I ) • VSECT(LL) 

IOECT(LL)  - ISECT(LL)  - XZIS{LL) 

REPEAT  STEPS  113,  114  AND  115 


I 
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Step  I lb  Calculate  Maximum  Section  Voltage 

Perform  a straight-line  interpolation  between  the  last  two  sets 
of  section  current-voltage  values  to  predict  the  voltage  at 
which  section  current  is  equal  to  z«ro  (VZCR) 

VSKCT(LL)  = VZCR 

iCKCT(LL)  =0.0 

XZIS(Ll)  = F{ZI(I >,  ZV ( I ) > at  ZV(I)  = VSECT(LL) 

NFINAL  = LL 

Step  117  Calculate  PSG  Current  Voltage  Characteristics 
SY(LL)  = ISECT(LL)  » NESP 
XZI(LL)  = XZIS(LL)  • NESP 
SX(LL)  = VSECT(LL) 

LL  * 1,  NFINAL 

Return  to  Performance  Analysis  Routine 

Step  118  Calculate  Occultation  Conditions  for  Solar  Insolation 
QDN  - 0.0 
QD  * 0.0 
QDG  « 0.0 
QOS  « 0.0 
QDT  « 0.0 

Where:  QDN  ■ Direct,  Normal  Solar  Insolation  at  Buoy 
Location  - Watts/M2 

QD  » Direct  Solar  Insolation  Incident  on  Solar 
Array  - Watts /M2 

QDG  * Horizontal  Surface  Diffuse  Insolation  Incident 
on  Solar  Array  - Watts/M2 

QDS  ■ Sky  Diffuse  Insolation  Incident  on  Solar 
Array  - Watta/M2 

QDT  ■ Total  Solar  Insolation  Incident  on  Solar 
Array  ~ Watta/M2 
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Ster  11 8a  Calculate  Solar  Array  Maximum  Power 


MSArWR  =0.0 


Step  11)  Obtain  Solar  Array  Parameters 

VSAINC  = Solar  Array  Voltage  Increment  - Volts 

KFINAL  = Maximum  extent  of  PSG  group  Current-Voltage 
Characteristics  Matrix 

Step  120  Calculate  PSG  Current-Voltage  Characteristics 
SY( 1 ) = 0.0 
XZI(l)  = 0.0 

SX(1)  = 0.0  (INITIALIZATION) 


SY(LL)  = 0.0 
XZI(LL)  =0.0 

SX(LL)  = SX( LL  - 1)  ♦ VSAINC 
LL  « 1,  KFINAL 

RETURN  TO  PERFORMANCE  ANALYSIS  ROUTINE 
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b.l.l  Shunt  Limiters 

The  Shunt  Limiter  routine  allows  the  user  to  specify  whether  an  active 
shunt  limiter,  an  ordinary  zener  diode,  a temperature-compensated  zener 
diode,  or  no  shunting  device  will  be  used  in  conjunction  with  the  solar 
array.  The  shunt  device  selected  will  not  influence  the  solar  array 
performance  until  prevailing  conditions  require  the  limiting  of  the 
array  voltage.  The  array  voltage  is  then  clamped  at  a maximum  potential , 
altering  the  Power  Source  Group  current-voltage  characteristics.  The 
combined  solar  array/shunt  limiter  performance  curve  is  the  algebraic 
difference  between  the  solar  array  and  the  array  limiter  characteristics. 


PROGRAM  ALGORITHMS 


Step  1 Obtain  Shunt  Limiter  Type 
ISH  * Shunt  Limiter  Type 

0 = No  Shunt  Limiter 

1 * Ordinary  Zener  Diode 

2 * Temperature  Compensated  Zener  Diode 

3 ■ Active  Shunt  Limiter 


Step  2 Select  Shunt  Limiter  Current-Voltage  Characteristics 


If: 

ISH  * 

0; 

GO  TO  STEP  3 

If: 

ISH  « 

li 

GO  TO  STEP  l* 

If: 

ISH  ■ 

2; 

GO  TO  STEP  5 

If: 

ISH  ■ 

3; 

GO  TO  STEP  6 
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Step  3 Calculate  No-Shunt  Limiter  Current  Voltage  Characteristics 
FOR  I * 1,  20 
21(1 ) = 0.0 
ZV(I)  = 0.0 

Where:  ZI(I)  = Shunt  Limiter  Current  at  ZV(I)  - Amperes 
ZV(l)  » Shunt  Limiter  Voltage  - Volts 
RETURIJ  to  power  sources  group  routine 

step  Obtain  Ordinary  Zener  Diode  Current-Voltage  Characteristics 

ZI(I ) ; ZV(I) 

For  I » 1,20 

RETURN  TO  POWER  SOURCES  GROUP  ROUTINE 

S’ep  5 Obtain  Temperature  Compensated  Zener  Diode  Current-Voltage 
Characteristics 

For  I « 1,  20: 

ZI(I);  ZV(I) 

RETURN  TO  POWER  SOURCES  GROUP  ROUTINE 

Step  6 Obtain  Active  Shunt  Limiter  Current-Voltage  Characteristics 
For  I » 1,  20: 

ZI(l)i  ZV(I) 

RETURN  TO  POWER  SOURCES  GROUP  ROUTINE 
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Step  1 


Step  2 


For  the  Initial  Load  Line  Analysis  Calculation 
go  to  Step  1 

For  Subsequent  Load  Line  Analysis  Calculations  go  to  Step  3 


Obtain  Zener  Diode  Operational  Requirements 

VZRR  = Breakdown  Voltage  of  a Single  Zener  Diode  - Volts 
(at  TZBR) 

TZBR  = Temperature  of  Zener  Diode  - °C 
(at  Breakdown  Voltage) 

Calculate  Zener  Diode  Breakdown  Voltage  at  Reference  Temperature 

Iterate  the  following  equations  until  the  change  in  Zener  breakdown 
voltage  at  reference  temperature  is  less  than  0.1  volts.  The 
number  of  iterations  shall  not  exceed  15. 


VZB3C  = VZBR  * 


1.0  - 


TC  = ZTCOEF  (VZB30) 

Where:  VZB30  = Zener  diode  breakdown  voltage  at  30°C  - Volts 

TC  * Zener  diode  temperature  coefficient  ( %/°C ) as  a 
function  of  VZB30 

ZTCOEF  * Input  table  of  TC  as  a function  of  VZB30 


Step  3 Obtain  Solar  Array  Temperature 

TSAC  = Solar  Array  Temperature  - "Centigrade 
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Step  b 


Step  5 


Step  (> 


Step  i' 


Step  8 


Step  9 


Calculate  Zener  Diode  Operating  Temperature 
TJZ1  = TSAC 

Where:  TJZ1  = Zener  Diode  Operating  Temperature  - °C 

Calculate  Zener  Diode  Breakdown  Voltage 
VZB  = VZB30  • [l.O  + TC-— 3°‘°)j 
TC  = ZTCOEF  {VZB 30} 

Where:  VZB  = Single  Zener  Diode  Breakdown  Voltage  - Volts 
at  TJZI 

Obtain  Number  of  Zener  Diodes  in  a String 
NZS  = Number  of  Zener  Diodes  in  Series 

Calculate  Zener  String  Breakdown  Voltage 
VZSB  = VZB  * NZS 

Where:  VZSB  = Zener  String  Breakdown  Voltage  - Volts 

Calculate  Zener  Diode  Dynamic  Impedance 
ZZ  = ZDIMP  {TJZI,  VZB30} 

Where:  7.D1MP  » Input  Table  of  ZZ  as  a function  of  TJZI  and  VZB30 


Calculate  Zener  Diode  Current-Volatge  Characteristics 

ZI(I)  * Zener  Diode  Current  - Amperes 
at  ZV(I) 

ZV(I)  * Zener  Diode  Voltage  - Volts 
I ■ 1,  20 
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Step  9 (contd) 

as  follows: 

a)  For  I = 1 

ZI(I) 

ZV(I) 

b)  For  I = 2 

ZI(2) 

ZV(2) 

c ) For  1 = 3 

ZI(3) 
ZV(  3) 

d)  For  I = U , 20 

ZI(I) 

ZV(I) 

Step  10 


= 0.0 
= 0.0 

= 0.0 

= VZSB 

= 100.0 

= VZSB  + ZI(3)  * zz  * NZS 

= 0.0 
= 0.0 


RETURN  TO  GENERAL  SHUNT  LIMITER  ROUTINE 
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TEMPERATURE  COMPENSATED  ZENER  DIODE 


For  the  Initial  Load  Line  Analyses  Begin  at  Step  1 
For  Subsequent  Load  Line  Analyses  Begin  at  Step  5 

Step  1 Obtain  Zener  Diode  Operational  Requirements 

NZS  = Number  of  Zener  diodes  in  series  in  a string 
HDZMX  * Maximum  Heat  Dissipation  of  a single  zener  diode  - Watts 
HDER  = Heat  Dissipation  Derating  Factor  for  a single  zener  diode 

Step  2 Calculate  Reference  Zener  Power 
PZRF25  = HDER  * HDZMX 

Where:  PZRF25  x Zener  diode  power  at  25°C  - Watts 

Step  3 Calculate  Reference  Zener  Current 
IZRF25  * CURZ  {HDZMX} 

Where:  IZRF25  * Zener  Diode  Current  at  25°C  - Watts 

CURZ  * Input  Table  of  IZRF25  as  a function  of  HDZMX 

Step  1*  Calculate  Reference  Zener  Voltage 
VZRF25  ■ PZRF25/IZRF25 

Where:  VZRF25  * Zener  Diode  Voltage  at  25°C  - Volts 

Step  5 Obtain  Solar  Array  Temperature 

TSAC  ■ Solar  Array  Temperature  - °C 
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Step.  6 Calculate  Zener  Diode  Operating  Temperature 


TCZ  = TSAC 


Where:  TCZ  = Zener  Diode  Operating  Temperature  - °C 


>tep  T 


Calculate  Zener  Diode  Breakdown  Voltage  Ratio  as  follows: 

a) .  Given  the  formulation: 

RATI  = TCZIV{RATV,  TCZ) 

Where:  RATI  = Zener  Diode  Current  Ratio  at  RATV  and  TCV 

RATV  = Zener  Diode  Voltage  Ratio 

TCZIV  = Input  table  of  Zener  diode  current-voltage 

characteristics  (RATI,  RATV)  as  a function  of  TCV 

b)  Then: 

RATVB  = RATV:  when;  RATI  *0.0 
Where:  RATVB  = Zener  Diode  Breakdown  Voltage  Ratio 


Step  3 


Calculate  Voltage  Ratio  Increment  Size 


RTVINC  * 1.05  - 


Step  9 Calculate  Zener  Diode  Current  Ratio-Voltage  Ratio  Characteristics 
RIZ(J)  ■ Zener  Diode  Current  Ratio  at  RVZ(J) 

RVZ(J)  * Zener  Diode  Voltage  Ratio 
J ■ 1,  20 
as  follows: 


a)  For  J « 1 


RIZ(l)  ■ 0.0 
RVZ(l)  « 0.0 


StepjM  contd) 

b)  For  J = 2 

RIZ(2)  = 0.0 
RVZ(2)  = RATVB 

c)  For  J = 3,  20:  Repeat  Step  80;  18  times 

RVZ(J)  = RVZ( J - 1)  ♦ RTVINC 
RATV  = RVZ(J) 

RATI  = TCZIV{RATV,  TCZ} 

RIZ(J)  * RATI 
J * J ♦ 1 

Step  10  Calculate  Zener  Diode  String  Current  Voltage  Characteristics 
ZI ( J ) * Zener  Diode  Current  at  ZV(J)  - Amperes 
ZV(J)  = Zener  Diode  String  Voltage  - Volts 
J = 1,  20 
as  follows: 

ZI(J)  = IZRF25  * RIZ(J) 

ZV(J)  « IJZS  • VZRF25  * RVZ(J) 

Step  11 


RETURN  TO  CENERAL  SHUNT  LIMITER  ROUTINE 
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ACTIVE  SHUNT  LIMITER 


Stop  1 Obtain  Shunt  Limiter  Operational  Requirements 

TSAC  = Solar  Array  Temperature  - °C 

VSHTOR  = Required  Shunt  Limiter  Turn-On  Voltage  - Volts 
at  TSHREF 

TSHREF  = Shunt-Limiter  Reference  Temperature  - °C 

CSH  = Shunt-Limiter  Turn-On  Voltage  Coefficient  - %/°C 

Step  2 Calculate  Shunt  Limiter  Operating  Temperature 
TSH  = TSAC 

Where:  TSH  = Shunt  Limiter  Operating  Temperature 


Step  7 Calculate  Shunt  Limiter  Turn-On  Voltage 

VSHTO  = VSHTOR  |\.0  + OSH  * (TSH.  - TSHREF) j 

Where:  VSHTO  = Shunt  Limiter  Turn-On  Voltage  - Volts 
(at  TSH) 

Step  h Calculate  Shunt  Limiter  Dynamic  Impedance 
ZSH  = ZSHTAB{TSH} 

Where:  ZSH  = Shunt  Limiter  Dynamic  Impedance  - Ohms 

ZSHTAB  * Input  Table  of  ZSH  as  a function  of  TSH 
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Step  5 Calculate  Shunt  Limiter  Current  Voltage  Characteristics 
ZI(I)  = Shunt  Limiter  Current  - Amperes 
ZV(I)  = Shunt  Limiter  Voltage  - Volts 
I - 1.  3 
as  follows: 

a)  For  I = 1 

ZI(1)  = 0.0 
ZV(1)  = 0.0 

b)  For  I = 2 

ZI(2)  = 0.0 
ZV ( 2 ) ■ VSHT0 

c)  For  I = 3 

ZI(3)  - 100.0 

ZV(  3)  * VSHTO  ♦ ZI( 3)  * ZSH 

d)  For  I = 1*,  20 

ZI(I)  « 0.0 
ZV(I)  ■ 0.0 

Step  6 RETURN  TO  GENERAL  SHUNT  LIMITER  ROUTINE 
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4.2  Kner^r  TU-rage  Group 

The  Energy  Storage  Group  is  made  up  of  the  buoy  batteries,  the  battery 
cables,  the  battery  chargers  and  the  battery  discharge  diodes.  The 

characteristics  of  these  elements  are  consolidated  as  a function  of  the 
battery  states-of-charge,  temperature,  number  of  series  cells,  cable 
resistance,  and  battery  charge  rate  and  are  then  expressed  as  a single 
set  of  current-voltage  characteristics  at  the  unregulated  bus  operating 
point. 

The  Energy  Storage  Group  algorithms  access  a comprehensive  set  of  battery 
charge  and  discharge  data  that  are  ordered  in  specific  battery  operating 
states  and  temperatures.  The  battery  data  are  in  the  form  of  current- 
voltage  curves  at  specified  temperatures  and  depths  of  discharge.  A total 
of  21  curves  are  available  for  each  of  six  temperatures,  126  curves  in  all. 

PROGRAM  ALGORITHMS 

Step  1 Obtain  Battery  Charger  Type 

ICHRT  = Battery  Charger  Type 

0 = Ho  Charger 

1 ■ Constant  Voltage  Charger  with  Current  Limit 

IF;  This  is  Initial  Load  Line  Analysis; 

THEN:  GO  TO  STEP  2 

IF:  Thi3  is  Subsequent  Load  Line  Analysis: 

THEN:  GO  TO  STEP  3 
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Step  2 Obtain  Battery  Characteristics 

DTTESG  * Energy  Storage  Group  Temperature  Rise  - °F 
NBATT  = Number  of  Batteries  in  Parallel 

CB  = Capacity  of  each  Battery  - ampere-hours 
XN  = Number  of  Cells  in  series  in  each  Battery 
RL  = Resistance  of  Cable  Connected  to  each  Battery  - Ohms 
XICHMX  = Maximum  Allowable  Battery  Charge  Current  - Amperes 

Step  3 Obtain  Ambient  Temperature 

TTAMB  = Ambient  Temperature  at  Selected  Location  - °F 

Step  L Calculate  Energy  Storage  Group  Temperature 

TTESQ  - TTAMB  ♦ DTTESG 

Where:  TTESG  = Energy  Storage  Group  Temperature  - °F 

Step  S Obtain  Battery  State-of-Charge 

ow  * ^ t nt of  r***  of  Mph 

Step  6 Calculate  Battery  Current-Voltage  Characteristics 

BRR(J,K,L)  ■ BCQT{VCC(J,K,L),  QBB(K),  TBR(L)} 

VC(J)  ■ VCC(J,K,L)  at  Q£B  ■ QB  and  TBB  ■ TTESG 
BR(J)  « BRR(J,K,L)  at  VCC  » VC(J),  QBB  « OB,  and  TBB  ■ TTESG 
VB( J ) « VC(J)  • XN 
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Step  6 (Ccntd) 

XIB(J)  = BR(J)  • CB 

J = 1,  9 (Data  Points) 

K = 1,  NQBB 
L ■ 1,  NTBB 

Where:  BRR  = Nonnaliaed  Battery  Current  Rate  (expressed  as  ratio  of 
battery  current  (amperes)  to  battery  capacity 
(ampere-hours))  - Hours 

VCC  * Cell  Voltage  - VDC 

QBB  * Battery  State-of-Charge 

TBB  = Battery  Temperature  - 8F 

BCQT  * Input  Table  of  BRR  us  a function  of  VCC,  QBB  and  TBB 

NQBB  * Humber  of  QBB  entries  in  BCQT 

NTBB  = Number  of  TBB  entries  in  BCQT 

VC( J ) « Cell  Voltage  for  each  Battery  - VDC 

* ' 

VB( J ) ■ Battery  Voltage  for  eac's  battery  - VDC 

BH(J)  * Normalized  Battery  Current  for  each  Battery  - Hours”* 

XIB(J)  ■ Battery  Current  for  each  Battery  - Hours-1 

Step  7 Calculate  Effect  of  Parasitic  Losses 

IF:  :<IB(J)  >0.0 

THEN:  VBM(J)  ■ VB(J)  ♦ RL  • XIB(J) 

IF:  XIB(J)  < 0.0 

THEN:  < VBM( J ) ■ VB(J)  - VDI0DE  - RL  • |XIB(J)| 

\VDI0DE  « AD2{ | X I B ( J ) | } 

Where:  VBM  ■ Modified  Battery  Voltage  for  each  Battery  - VDC 

VDI0DE  * Battery  Discharge  Blocking  Diode  Voltage  Drop  - VDC 

AD2  ■ Input  Table  of  VDI0DE  as  a function  of  battery 
discharge  current 

J ■ 1,  9 
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Step  7 (Contd) 

IF:  ICHRT  = 1 

THEN : GO  TO  STEP  17 

Stei:  8 Rearrange  Modified  Battery  Voltage  Data  into  One  Array  in  Ascending 
Order 

TRESLT(J,l)  « F{’^BM( J ) } (in  ascending  order) 

J • 1,  9 

Where:  TRESLT(J,1)  ■ Energy  Storage  Group  Voltage  - VDC 
Step  9 Initialize  Counter  and  Special  Voltage  Array 
LL  « 1 
LY  « 1 

TRESVT(LY)  « TRESLT(LL.l) 

Where:  TRESVT  « Special  Voltage  Array  for  ES  Group  - VDC 

LY  ■ Remaining  Number  of  Modified  Voltage  Points 
Step  10  Index  Counter  and  Compare  Voltage  Differences 
LL  B LL  ♦ 1 

IF:  | (TRESLT(LL.l)  - TRESLT((LL-l),l)|<  0.01 

THEN:  REPEAT  STEP  10 

Step  11  Index  Counter  and  Calculate  Values  in  Special  Voltage  Array 
LY  ■ LY  ♦ 1 

TRESVT (LY)  - TRESLT(LL.l) 

IF:  LL  > 9 

THEN:  CO  TO  STEP  12 

REPEAT  STEPS  10  AND  11 

Step  12  Calculate  Revisions  to  Energy  Storage  Group  Voltage 

TRESLT(LL.l)  -0.0 
LL  • 1,9 

TRESLT(LY.l)  ■ TRESVT(LY) 

LY  - 1,9 
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Step  13  Calculate  Individual  Battery  Voltage 
TRESV(LY)  = TRESLT(LY.l) 

LY  = 1,9 

Where:  TRESV  = Individual  Battery  Voltage  - VDC 
Step  1*4  Calculate  Average  Battery  Cell  Voltages 
VCM(J)  = VCM(J)/XN 
J = 1,9 

Where:  VCM  = Modified  Average  Cell  Voltage  - VDC 

Step  15  Calculate  Individual  Battery  Current 

TRESI(J)  = F{XIB( J) ,VCM( J ) } at  values  of  VCM(J)  = TRESV(J)/XN 
J = 1,9 

Where:  THESI  - Individual  Battery  Current  - Amperes 
Step  16  Calculate  Total  Battery  Current 
TRESLT(LY,2)  = TRESI(LY)  * NBATT 
LY  = 1,9 

Where:  TRESLT(LY,2)  = Energy  Storage  Group  Current  - Amperes 

at  TRESLT(LY ,1 ) 

TRKSLT( LY , 1 ) = Energy  Storage  Group  Voltage  - VDC 
RETURN  TO  PERFORMANCE  ANALYSIS  ROUTINE 
Step  17  Enlarge  Modified  Current-Voltage  Characteristic  Array 
XIBBM(JJ)  ■ XIB(J) 

VBUM(JJ)  » VBM(J) 

JJ  ■ J 
J • 1,9 

But,  allow  for  an  extra  location  so  that  (JJ)  . ■ 10 

maximum 
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Step  1 6 Insert  Battery  Voltage  Corresponding  to  Maximum  Allowable  Battery 

Current  (XICHMX)  into  Modified  Current-Voltage  Characteristic  Array 
such  that  the  array  contains  I-V  data  for  that  point. 

Thus, 

VBBM(JJ)  = Modified  Battery  Voltage  for  each  Battery  - VDC 

XIBRM(JJ)  - Battery  Current  for  each  Battery  - Ampere 
Corresponding  to  VBBM(JJ) 

JJ  = 1,10 

Step  19  Obtain  Battery  Charger  Reference  Voltages 
VCHIO  = VCHI0T{TTESG} 

VCHISA  = VCHIST{TTESG} 


Where:  VCHIO  = Battery  Charger  Input  Voltage  at  Turn-On  - VDC 

(Minimum  Voltage  Drop  at  Zero  Current  Level) 

VCHISA  = Battery  Charger  Input  Voltage  - VDC 

at  which  operation  changes  from  "saturated"  to 
"active"  conditions 


VCHIOT  * Input  table  of  VCHIO  as  a function  of  TTESG 

VCHIST  * Input  table  of  VCHISA  as  a function  of  TTESG 
Step  20  Obtain  Battery  Charger  Impedance 

ZCHRS  = ZCHRST{ TTESG} 


ZCHRA  = ZCHRSA{ TTESG) 

Where:  ZCHRS  ■ Output  Impedance  of  Battery  Charger  in  "saturated" 
condition  - Ohms 

ZCHRA  * Output  Impedance  of  Battery  Charger  in  "active" 
condition  - Ohms 


ZCHRST  * Input  Table  of  ZCHRS  as  a function  of  TTESG 
ZCHRAT  * Input  Table  of  ZCHRA  as  a function  of  TTESG 
Step  21  Obtain  Battery  Charger  Impedance 
ZCHRS  ■ ZCHRST{ TTESG) 

ZCHRA  « ZCHRSA{ TTESG) 
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Step  21  (Cont) 

Where:  ZCHRS  = Output  Impedance  of  Battery  Charger  in  "saturated" 
condition  - Ohms 

ZCHRA  = Output  Impedance  of  Battery  Cnarger  ir:  "active" 
condition  - Ohms 

ZCHRST  = Input  Table  of  ZCHRS  as  a function  of  TTESG 
ZCHRAT  = Input  Table  of  ZCHRA  as  a function  of  TTESG 
Step  2j  Initialize  Counter 
JJ  = 1 

Step  2h  Calculate  Energy  Storage  Unit  Discharge  I-V  Chracteristics  and  Then 
Increment  Counter 

IF:  XIBBM(JJ)  _>  0.0 

THEM : GO  TO  STEP  25 

VESI(JJ)  = VBBM(JJ) 

XIEST(JJ)  = XIBBM(JJ) 

JJ  = JJ  + 1 

Where:  VESI(JJ)  = Energy  Storage  Unit  Input  Voltage  - VDC 

XIESI(JJ)  = Energy  Storage  Unit  Input  Current 

(Corresponding  to  VESI(JJ))  - Amperes 

Repeat  Step  2h 

Step  <?5  Compare  Power  Sources  Group  Maximum  Voltage  with  Reference  Voltages 

IF:  (VPSGMX ) > (VCHIO  ♦ VBBM(JJ)) 

THEN:  GO  TO  STEP  28 

Step  P6  Calculate  Energy  Storage  Unit  Charge  I-V  Characteristics  and 
Increment  Counter 

IF:  JJ  > 10 

THEM:  GO  TO  STEP  27 

VESI ( JJ ; ■ VBBM(JJ) 

XIESI(JJ)  ■ 0.0 
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REPEAT  STEP  26 
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Step  28 


Step  29 


Step  30 


Step  31 


Calculate  Estimate  of  Charger  Input  Voltage 

VCHOOS  = VBBM(JJ)  + XIBBM(JJ ) * ZCHRS 

VESS  = VCHOOS  + VCHIO 

IF:  VESS  >_  VCHISA 

THEN:  CO  TO  STEP  30 

Where:  VCHOOS  = Battery  Charger  Output  Voltage  (in  "saturated" 
condition)  at  zero  current  - VDC 

VESS  = Estimate  of  Battery  Charger  Input  Voltage  in 
"saturated"  condition  - VDC 

Calculate  Energy  Storage  Unit  Charge  I-V  Characteristics  and 
Increment  Counter 

IF:  JJ  > 10 

THEN:  CO  TO  STEP  27 

XIESI(JJ)  = XIBHM(JJ) 

VESI(JJ)  = VESS 

JJ  * JJ  ♦ 1 

RETURN  TO  STEP  28 

Calculate  Estimate  of  Charger  Input  Voltage 
VCHOOA  * VBBM(JJ)  + XIBBM(JJ)  » ZCHRA 
VESA  ■ VCHIT{ VCHOOA, TTESG) 

Wiiere:  VCHOOA  * Battery  Charger  Output  Voltage  (in  "active" 
condition)  at  zero  current  - VDC 

VESA  * Estimate  of  Battery  Charger  Input  Voltage  in 
"active"  condition  - VDC 

VCHIT  * Input  Table  of  VESA  as  a function  of  VCHOOA  and 
TTESG 

Calculate  Energy  Storage  Unit  Charge  I-V  Characteristics  and 
Increment  Counter 

IF:  JJ  > 10 

THEN : CO  TO  STEP  27 
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Step  31  (Cuntd) 

XIESI(JJ)  = XIBBM(JJ) 

VESI(JJ)  = VESA 

JJ  = JJ  + 1 
RETURN  TO  STEP  30 

Ste£  32  Compare  Battery  Charger  Input  Current  with  Reference  Current  Limit 

IF:  XIESI(JJ)  > XIC121X 

THEM : XIESI(JJ)  = XICHMX 

JJ  =1,10 

Jt-pp  33  Rearrange  Energy  Storage  Unit  Voltage  Data  into  One  Array  in 
Ascending  Order 

TRESLT(JJ,i)  = F{VESI ( J J ) } (in  ascending  order) 

JJ  = 1,  10 

Step  3l*  Initialize  Counter  and  Special  Voltage  Array 
LL  = 1 
LY  = 1 

TRESVT(LY)  « TRESLT{LL,1 ) 

Step  35  index  Counter  and  Compare  Voltage  Differences 
LL  * LL  ♦ 1 

JF:  | (TRESLT(LL,1 ) - TRESLT((LL  - 1),1)|  < 0.01 

THEM:  REPEAT  STEP  35 

Step  jb  Index  Counter  and  Calculate  Values  in  Special  Voltage  Array 
LY  » LY  ♦ 1 

TRESVT(LY)  = TRESLT(LL.l) 

IF:  LL  > 10 

n'HEM : GO  TO  STEP  37 
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REPEAT  STEPS  35  AND  36 
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Step  3J  Calculate  Revisions  to  Energy  Storage  Group  Voltage 
TRESLT(LL.l)  =0.0 
LL  = 1,10 

TRESLT(LY.l)  = TRESVT(LY) 

LY  = 1,10 

Step  38  Calculate  Individual  Energy  Storage  Unit  Voltage 
TRESV(LY)  = TRESLT(LY,l) 

LY  = 1,10 

Step  39  Calculate  Energy  Storage  Unit  Equivalent  Single  Cell  Voltage 
VESIC(JJ)  * VESI ( JJ )/XN 
JJ  = 1,10 

Where:  VESIC  = Equivalent  Single  Cell  Voltage  of  an  Energy  Storage 
Unit  - VDC 

Step  ^0  Calculate  Energy  Storage  Unit  Current 
TRESI(LY)  * F{XIESI(JJ ), VESIC (JJ )} 

At  values  of:  VESIC (JJ)  ■ TRESV(LY }/XH 

JJ  ■ 1,10 
LY  • 1,10 

Step  kl  Calculate  Total  Energy  Storage  Group  Current 
TRESLT(LY.P)  ■ TRESI(LY)  • HBATT 
LY  » 1,10 

RETURN  TO  PERFORMANCE  ANALYSIS  ROUTINE 
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U . 3 Power  Conditioning  and  Distribution  Group 


The  Power  Conditioning  and  Distribution  Group  is  made  up  of  two  sub- 
assemblies:  the  Lamp  Flasher  and  the  Housekeeping  Regulator.  The  char- 
acteristics of  the  subassemblies  are  computed  as  a function  of  the  lamp 
flasher  pattern  and  the  lamp  flasher  condition  (on,  off,  or  flashing). 
These  characteristics  are  then  shifted  for  the  combined  effects  of 
wiring  and  connector  series  resistance  to  give  a single  set  of  c irre.nt- 
voltage  curves  at  the  unregulated  bus. 


PROGRAM  ALGORITHMS 


Step  1 


Jtep  2 


Step  3 


For  the  Initial  Load  Line  Analysis,  GO  TO  STEP  1 
For  Subsequent  Load  Line  Analysis,  GO  TO  STEP  13 

Obtain  Flasher  Pattern  Type 

IF:  IFTYPE  = 0,  GO  TO  STEP  3 

Where:  IFTYPE  * Type  of  flasher  pattern 
IFTYPE  * 0:  Non-Standard  Pattern 
IFTYPE  > 0:  Standard  Pattern 

Calculate  Standard  Flasher  Pattern 
TL1(J)  = TL0( IFTYPE ,J ) 

(1  ^ IFTYPE  <_  15)  (15  standard  pattern  types) 

(1  < J < 16) 


i Up  to  16  steps  per  pattern i 
l Alternate  On/Off  f 


Where:  TL0( IFTYPE, J ) * Input  table  containing  the  patterns 

exhibited  by  the  Standard  Lamp  Flashers 


TL1  ■ Selected  Lamp  Flasher  Pattern 


GO  TO  STEP  4 


Calculate  Non-Standard  Flasher  Pattern 

TL1(J)  ■ TLL1(J)  (1  < J < 16) 

Where:  TLL1(J)  ■ Input  data  containing  up  to  16  alternate  on-off 
steps  for  the  Non-Standard  Flasher  Pattern 


Step  I4  Calculate  Total  Duration  of  Lamp  Illumination  and  Lamp  Shut-Off 

15 

TLON  = ^ TL1(J) 

J=1 ,3,5. . . 

16 

TLOFF  = ^ TL1(J) 

J=2 . i*  .6 . . • 

Where:  TLON  = Total  duration  of  lamp  illumination  | in  a single 

< flasher 

TLOFF  = Total  duration  of  lamp  shut-off  | period 

IF:  TLON  0 and  TLOFF  <_  0 stop  program  and 
Print:  "No  flasher  pattern  entries" 

Step  5 Calculate  Lamp  Duty  Cycle 
DL  * TLON/ (TLON  ♦ TLOFF) 

Where:  DL  x Lamp  Duty  Cycle 

Step  6 Obtain  Lamp  Characteristics 

VLR  ■ Lamp  Voltage  Rating  - VDC 

CLR  = Lamp  Current  Rating  - Amperes 

CL3  * Cold-Filament  Lamp  Surge  Coefficient 

Step  7 Calculate  Actual  Lamp  Current 

IL  - CLS  *CLR 

Where:  IL  * Actual  Lamp  Current  -Amperes 


Note:  IF: 

DL  ■ 1.0 

THEN: 

CLS  ■ 1.0 

IF: 

DL  < 1.0 

THEN: 

CLS  > 1.0 
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3 top  8 Calculate  Actual  Lamp  Resistance 

PI  _ VLR 
IL 

Where:  RL  = Actual  Lamp  Resistance  - Ohms 

Otep  9 Calculate  Average  Lamp  Current 
Tl  = IL  * DL 

Where:  IL  = Average  Lamp  Current  - Amperes 

Step  10  Calculate  Effective  Lamp  Resistance 
RL  = VLR/IL 

Where:  RL  = Effective  Lamp  Resistance  - Ohms 

Otep  11  Obtain  Raw  Rower  Bu3  Voltage  Limits  and  User  Load  Cable  Resistance 
VMINIV  = Minimum  Raw  Power  Bus  Voltage  - VDC 
VMAXIV  = Maximum  Raw  Power  Bus  Voltage  - VDC 
RLL  * User  Load  Cable  Resistance  - Ohms 

Step  1?  Calculate  PCD  Group  Voltage  Increment 
VINCIV  * (VMAXIV  - VMINIV)/50.0 

Where:  VINCIV  « PCD  Group  Voltage  Increment  - VDC 

Step  13  Obtain  PCD  Equipment  Temperature  Characteristics 
TTAMB  ■ Ambient  Temperature  - #F 
DTTPCD  ■ PCD  Equipment  Temperature  Rise  - °F 
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Step  ll*  Calculate  PCD  Equipment  Temperature 

rrrcD  = ttamb  + dttpcd 

Where:  TTPCD  = PCD  Equipment  Temperature  - °F 

Step  i5  Compare  Raw  Power  Bus  Minimum  Voltage  With  Reference 

IF:  VMINIV  < VRIO 

THEN:  00  TO  STEP  16 

IF:  (VMINIV  > VRIO)  and  (VMINIV  <_  VRISA) 

THEN:  GO  TO  STEP  24 

IF:  VMINIV  > VRISA 

THEN:  GO  TO  STEP  29 

Where:  VRIO  = Minimum  (No  Current)  Voltage  Drop  - VDC 

Across  Lamp  Regulator  in  "Saturated"  Condition 

VRISA  ■ Voltage  level  at  which  lamp  regulator  - VDC 
changes  from  "Saturated"  condition  operation 
to  "Active"  operation 

VRIO  ■ VRIOTfTTPCD) 

VRISA  « VRISATfTTPCD) 

Step  l6  Initialize  Counter  and  Lamp  Regulator  Voltage 
J * 1 

VRI(J)  * VMINIV 

Step  17  Calculate  Lamp  Regulator  Current 
IRI(J,1 } ■ 0.0 
IFI(J,P)  - 0.0 
IRI(J,3)  » 0.0 

Where:  VRI(J)  * Lamp  Regulator  Input  Voltage  - VDC 

IRI(J,K)  ■ Lamp  Regulator  Input  Current  - Volts 
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Stop  17  (contd) 


Step  18 


Step  19 


Step  20 


When:  K = 1 - Lamp  Off 

K = 2 - Lamp  Flashing  - Effective 
K = 3 - Lamp  On 


Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare 
With  Reference 

J = J + 1 

VRI(J)  * VRI ( J - 1)  ♦ VINCIV 


IF:  (VRI(J)  > VRIO)  and 
IF:  (VRI(J)  < VMAXIVT" 
THEN:  CO  TO  STEP  20 

IF:  VRI(J)  > VMAXIV 
THEN:  CO  TO  STEP  32 


Calculate  Lamp  Regulator  Currents 

IRI ( J, 1 ) ■ 0.0 

IRI(J,2)  * 0.0 

IRI(J,3)  ■ 0.0 

REPEAT  STEPS  18  AND  19 


Calculate  Lamp  Regulator  Current 
IRI(J,1)  - 0.0 

IRI{ J,2)  ■ (VRI(J)  - VRIO)/ {RL  ♦ ZRS) 

IRI(J,3)  ■ (VRI(J)  - VRI0)/(RL  ♦ ZRS) 

Where:  ZRS  ■ Regulator  Impedance  in  "Saturated"  Condition  - Ohms 


ZRS  ■ ZRST{TTPCD) 

ZRST  ■ Input  Table  of  ZRS  as  a function  of  TTPCD 


Step  21 


Step  22 


Step  23 


Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare  with 
Reference  ' 

J = J ♦ 1 

VRI(J)  = VRI ( J - 1)  ♦ VINCIV 

IF:  (VRI { J ) > VRISA)  and 

IF:  (VRI(J)  < VMAXIV ) 

THEN:  Go  to  Step  22 

IF:  VRI ( J ) > VMAXIV 

THEN:  GO  TO  STEP  32 

REPEAT  STEPS  20  AND  21 


Calculate  Lamp  Regulator  Currents 

IRI(J,1)  * 0.0 

IRI(J,2)  * Vlb/(RL  ♦ ZRA) 

IRI(J,3)  ■ Vlb/(rl  ♦ ZRA) 

Where:  V^B  * Regulator  Output  Voltage  at  Zero  Current  - Volts 
ZRA  * Regulator  impedance  in  "Active"  region  - Ohms 
VLB  ■ VLBTfVRI,  TTPCD) 

ZRA  - ZRAT{  TTPCD) 

Where:  VLBT  ■ Input  Table  of  VLB  a3  a function  of  VRI  and  TTPCD 
ZRAT  ■ Input  Table  of  ZRA  as  a function  of  TTPCD 

Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare  With 
Reference 

J ■ J ♦ 1 

VRI(J)  ■ VRI( J - 1)  ♦ VINCIV 

IF:  VRI(J)  > VMAXIV 

THEN:  CO  TO  STEP  32 


REPEAT  STFPS  22  AND  23 
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Step  21*  Initialize  Counter  and  Lamp  Regulator  Voltage 
J = 1 

VRI(J)  = VMIIIIV 

Step  25  Calculate  Lamp  Regulator  Currents 

I R I ( J , 1 ) = 0.0 

I R I ( J , 2 ) = (VRI(J)  - VRIO  )/(PL  + ZRS) 

IRI(.I,3)  = (VRI(J)  - VRI0)/(RL  + ZRS) 

Step  Cti  Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare  with 
Reference 

J = J ♦ 1 

VRI(J)  = VRI ( J - 1)  ♦ VINCIV 

IF:  (VRI(J)  > VRISA)  and 

JT:  ( VRI (J ) < VMAXTV) 

THEN : CO  TO  STEF  27 

IF:  VRI(J)  > VMAXIV 

THEIi  CO  TO  STEP  32 

REPEAT  STEPS  25  AN!)  26 

Step  27  Calculate  Lamp  Regulator  Currents 
IHI(J.l)  » 0.0 
IR1(J,2)  « VLB/ (RL  ♦ ZRA) 

1HI(J,D  = VLB/ (HI  ♦ ZRA) 


Step  28  Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare  with 
Reference 

J ■ J ♦ 1 

VRI(J)  = VRI ( J - 1)  + VINCIV 

IF:  VRI(J)  > V MAX IV 

THEN:  GO  TO  STEP  32 

REPEAT  STEPS  27  AND  28 

S*;ep  29  Initialize  Counter  and  Lamp  Regulator  Voltage 
J * 1 

VRI(J)  * VMINIV 

Step  30  Calculate  Lamp  Regulator  Currents 
IRI(J,1)  * 0.0 
IRI(J,2)  ■ VLB/ (rIT  ♦ ZRA ) 

IRI(J,3)  * VLB/ ( RL  ♦ ZRA) 

Step  31  Increment  Counter  and  Lamp  Regulator  Voltage  and  Compare  with 
Reference 

J » J ♦ 1 

VRI(J)  « VRI ( J - 1)  ♦ VINCIV 

IF:  VRI(J)  > VMAXIV 

THEN:  GO  TO  STEP  32 

REPEAT  STEPS  30  AND  31 
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Step  32  Calculate  PCD  Group  Current 
XI ( J ,K ) = IHI(J)  + IRI  ( J ,K ) 

IHI(J)  = IHIT{VRI(J),  TTPCD} 

J = l,  50 

K = 1,  3 

Where:  XI(J,K)  = PCD  Group  Current  - Amperes 

IHI(J)  = Housekeeping  Load  Regulator  Input 
Current  - Amperes 

IHIT  = Input  Table  of  IHI(J)  as  a function  of  VRI(J) 
and  TTPCD 

Step  33  Calculate  PCD  Group  Voltage 

X X ( J , K ) = VRI(J)  + XI(J,K)  * RLL 
J * 1,  50 
K * 1,  3 

Where:  XX(J,K)  = PCD  Group  Voltage  - VDC 
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5.  MERGE 


The  MERGE  program  set  is  a group  of  computer  programs  (TDFlb,  DECK280, 
and  LISTMERGE)  which  provide  a support  function  to  the  DSPA  program.  The 
purpose  of  the  MERGE  package  is  to  provide  a means  for  extracting  actual 
weather  information  from  NOAA  data  tapes  and  using  it  as  input  to  the 
DSPA  program.  The  three  MERGE  programs  and  their  functions  are  described 
below. 


i > 

J 


r- 


11  ST  MERGE 

( MERGE  ( 

PROGRAM 

J LISTING  j 

FIGURE  5-1.  MERGE  COMPUTER  PROGRAMS  OVERVIEW 
5.1.  Creation  of  a MERGE  File  (TDF14) 


The  building  of  the  skeletal  MERGE  file  is  controlled  by  the  TDF14 
program.  The  user  requests  the  creation  of  a MERGE  file  to  span  a parti- 
cular period  of  years  for  a selected  location  from  a NOAA  TDF-ll  weather 
tape.  The  TDFll*  program  then  extracts  the  date,  temperature,  and  wind 
velocity  data  from  the  NOAA  tape,  builds  a one-day  record  consisting  of 
2lt  hourly  observations  of  temperature  and  wind  velocity  and  space  for 
solar  insolation,  and  sequentially  writes  the  day's  information  to  a 
MERGE  file. 
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PROGRAM  ALGORITHMS 


2 

Obtain  one  day's  weather  data  from  TDF-ll*  input  tape 

Step 

ry 

Convert  TDF-lk  data  to  numerics  (DECODE) 

Step, 

2. 

Write  MERGE  file  record 

Step 

2 

If:  end  of  TDF-ll*  tape. 

Then:  Stop  Program 

Otherwise:  Go  TO  STEP  1 

5.2  Addition  of  Solar  Insolation  Data  to  MERGE  File  (DECK280) 


The  addition  of  solar  insolation  data  to  a MERGE  file  created  by  TDFll* 
is  performed  by  the  DECK20O  program.  The  user  requests  that,  for  a 
particular  MERGE  location,  NOAA  DECK-280  tape  data  from  a specified  loca- 
tion be  inserted  into  the  file.  The  DECK280  program  extracts  the  solar 
radiation  data  (in  Langleys)  from  the  NOAA  tape,  converts  the  data  to 
vatts/Bquare  meter,  and  adds  the  data  to  the  appropriate  day  and  hour 
position  in  the  MERGE  file. 
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FIGURE  5-3.  DECK280  COMPUTER  PROGRAM 
PROGRAM  ALGORITHMS 


Step  1 Obtain  solar  radiation  data  from  DECK-280  input  tape 

Step  2 Convert  DECK-280  data  to  numerics  (DECODE) 

Step  3 Convert  DECK-280  data  to  watts/square  meter 

Q = S * 41.82/3.6 

where:  S ■ Solar  Insolation  - Langleys 

Q * Solar  Insolation  - watts/meter^ 

Step  4 Add  solar  insolation  data  to  appropriate  day  and  hour  in 
MERGE  file 

Step  5 If:  end  of  Deck-280  tape. 

Then:  Stop  Program 
Otherwise:  GO  TO  STEP  1 


5.3  D ■ -iplaying  MERGE  File  Data  (LISTMERGE) 

The  averge  MERGE  file  consists  of  from  10  to  12  years  of  hourly  tempera- 
ture, wind  velocity,  and  solar  insolation  data  recordings.  The 
LISTMERGE  program  permits  the  user  to  randomly  view  any  number  of 
sequential  days  within  the  file  beginning  at  any  date  contained  in  the 
file. 
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6. 


STAT 


The  STAT  program  set  is  a group  of  computer  programs  (STATS  and  PROFILE) 
which  provide  a support  function  to  the  DSPA  program.  The  purpose  of  the 
STAT  package  is  to  provide  a means  of  producing  a single  year  of  statis- 
tical data  from  a MERGE  file,  and  using  it  as  input  to  the  DSPA  program. 
The  two  STAT  programs  and  their  functions  are  described  below. 


FIGURE  6-1.  STAT  COMPUTER  PROGRAMS  OVERVIEW 


6.1  Statistical  Analysis  of  MERGE  File  Data  (STATS) 

The  STATS  program  uses  the  10  to  12  years  of  MERGE  file  weather  data  to 
produce  a one-year  statistical  file.  The  procedure  involves  the  averaging 
of  data  for  a given  hour  of  a given  day  of  each  of  the  years  contained  in 
the  MERGE  file.  The  statistical  data  is  then  written  to  a STATS  file 
which  is  used  as  input  to  the  PROFILE  program  (see  Section  6.2  below). 
Specifically,  the  STATS  program  computes  and  outputs  (both  to  the  STATS 
file  and  to  a printer)  the  following  statistics: 

la)  Average  temperature  for  each  hour  of  one  year. 

lb)  Average  wind  velocity  for  each  hour  of  one  year. 

lc)  Average  solar  insolation  for  each  hour  of  one  year. 

2a)  Average  wind  velocity  for  each  day  of  each  data  year. 

2b)  Average  solar  insolation  for  each  day  of  each  data  year. 

3a)  Average  wind  velocity  for  each  day  of  one  year. 

3b)  Average  solar  insolation  for  each  day  of  one  year. 
l*a)  Average  temperature  for  each  month  of  each  data  year. 
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4b  Average  wind  velocity  for  each  month  of  each  data  year. 

4c)  Average  solar  insolation  for  each  month  of  each  data  year. 
5a)  Average  temperature  for  each  month  of  one  year. 

5b)  Average  wind  velocity  for  each  month  of  one  year. 

5c)  Average  solar  insolation  for  each  month  of  one  year. 

6a)  Standard  deviation  of  statistics  gathered  in  4a. 

6b)  Standard  deviation  of  statistics  gathered  in  4b. 

6c)  Standard  deviation  of  statistics  gathered  in  4c. 

7a)  Maximum  temperature  for  each  year. 

?b)  Minimum  temperature  for  each  year. 

8a)  Mean  and  standard  deviation  of  statistics  gathered  in  7a. 

8b)  Mean  and  standard  deviation  of  statistics  gathered  in  7b. 


PROGRAM  ALGORITHMS 

Step  1 Initialize  yearly  minimum  and  maximum  temperatures 

TMin(J)  = 1000.0 
TMAX(J)  = -1000.0 
for  J*l,  12 

Step  2 Initialize  monthly  sums  for  each  year 

MTSUM(J)  = 0.0 
MVSUM(J)  * 0.0 
MQSUM(J)  =0.0 
for  J*l,  12 

where:  MTSUM(J)  ■ monthly  temperature  total  for  year  J 

MVSUM(J)  ■ monthly  wind  velocity  total  for  year  J 
MQSUM(J)  ■ monthly  solar  insolation  total  for  year  J 
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FIGURE  6-2.  STATS  COMPUTER  PROGRAM 


>v 


5040-27  (Change  1) 


Step  3 

Initialize  daily  sums  for  each  year 

DVSUM(J)  = 0.0 

DQSUM(J)  = 0.0 

for  J=l,  12 

where:  DVSUM(J)  = daily  wind  velocity  total  for  vear  J 

DQSUM(J)  ■ daily  solar  insolation  total  for  year  J 

Step  1* 

Obtain  hourly  temperature,  wind  velocity,  and  solar  insolation 

data  for  same  day  of  each  MERGE  year  (liYRS  < 12) 

T(I,J)  = temperature  for  hour  I of  year  J 
V(I,J)  = wind  velocity  for  hour  I of  year  J 
Q(I,J)  = solar  insolation  for  hour  I of  year 

J 

Step  5 

Initialize  hour  counter 

I = 1 

Step  6 

Initialize  hourly  sums 

TSUM  * 0.0 

VSUM  ■ 0.0 

QSUM  ■ 0.0 

where:  TSUM  ■ hourly  temperature  total  over  all  years 

VSUM  * hourly  wind  velocity  total  over  all  years 

QSUM  ■ hourly  solar  insolation  total  over  all  years 

Step  7 

Initialize  year  counter 

J ■ 1 

Step  8 

Add  hour's  data  to  sums 

TSUM  ■ TSUM  ♦ T(I,J) 

VSUM  - VSUM  ♦ V(I,J) 

QSUM  - QSUM  ♦ (Q(I,J) 

DVSUM{ J ) ■ DVSUM(J)  ♦ V{7,J) 

DQSUM(J)  ■ DQSUM(J)  ♦ Q'I,J) 

MTSUM(J)  ■ MTSUM(J)  ♦ T(I,J) 

MVSUM{ J ) « MVSUM(J)  ♦ V(I,J) 

MQSUM(J)  - MQSUM(J)  ♦ Q(I,J) 
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Steg  9 Select  minimum  and  maximum  temperatures 

TMIN(J)  = AM  INI  ( TMIIi  { J ) ,T(  I ,J ) ) 

TMAX { J ) = AMAX1 ( TMAX ( J ) , T ( I , J ) ) 

Step  10  If : J = NY RS 

Then:  GO  TO  STEP  11 
Otherwise:  J=J+i 
And:  GO  TO  STEP  8 

Step  11  Compute  hourly  statistics 

TOUT(I)  = TSUM/NYRS 
VOUT(I)  = VSUM/NYRS 
QOUT(I)  = QSUM/HYRS 

where:  TOUT(I)  = average  temperature  for  hour  I 

VOUT(I)  * average  wind  velocity  for  hour  I 
QQUT  (I)  = average  solar  insolation  for  hour  I 

Step  IP  If:  I = 2U 

Then:  GO  TO  i 3 
Otherwise:  1=1+1 
And:  TO  TO  STEP  6 

Step  13  Compute  daily  statistics  for  each  year 

VOAY(J)  = DVSUM(J)/(2fc  • NYR5) 

QDAY(J)  = DQSUM(J)/(2li  * NYRS) 
for  J=1,NYRS 

where:  VDAY(J)  * average  wind  velocity  for  given  day  of  year  J 

viDAY(J)  ■ average  solar  insolation  for  given  day  of 
year  J 

Step  ll  Write  statistical  day  record  to  STATS  file 

Step  15  If:  last  day  of  month 
Then:  GO  TO  STEP  l6 
Otherwise:  GO  TO  STEP  3 
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Tte;  +6  Compute  monthly  statistics  for  each  year 

TOUT ( J ) = MTSUM(J)/(2U*NYRS*NDYS) 

YOUT(J)  = MVSUM(J)/(2U*NYRS*NDYS) 

QOUT(J)  = MQSUM( J )/( 2U*NYRS*KDYS ) 
for  J=1  ,tiYR5 

where:  TOUT(J)  = average  temperature  for  given  month  of  year  J 
VOUT(J)  = average  wind  velocity  for  given  month  of 
year  J 

QO'uT(J)  = average  solar  insolation  for  given  month  of 
year  J 

NDYS  = n -amber  of  days  in  current  month 
Step  l7  Compute  mean  monthly  statistics  and  standard  deviations 


NYRS 

XQUT(13)=^  XOUT(J)/IIYRS 

J-l 


/ i 

i riYHS  \ 

\ 1 

XOUTdM  * SQRT 

( £ X0UT(J)»»2  - NYRS*  1 

\ J-l  ' 

£ XOUT(J) 1**2 

U«i  / 

J / (NYRS- ') 

where:  XOUT  is  used  to  represent  each  of  the  TOUT,  VQUT,  and 
(<QUT  variables  since  the  equations  are  of  the  same 
form  for  each 


■'tep  18  Write  statistical  month  record  to  STATS  file 


Step  19  If:  last  day  of  year 
Then:  GO  TO  STEP  20 
Otherwise:  GO  To  STEP  2 
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Step  20 


TOUT ( 2 ) 


TOUT  ( 1* ) = 


Step  21  Write  yearly  statistical  data  to  STATS  file 
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6.2  Environmental  Profiling  of  STATS  Data  (PROFILE) 


The  PROFILE  program  uses  the  one  year  of  statistically  prepared  STATS 
data  to  produce  a modified  statistical  file  for  use  with  the  DSPA 
program.  The  procedure  involves  the  scaling  of  the  STATS  file  data,  on  a 
monthly  basis,  by  factors  computed  from  user-specified  proportions  and 
confidence  levels.  In  addition,  the  PROFILE  program  will,  upon  request, 
perform  a worst  case  analysis  for  low  solar  insolation,  low  wind,  or 
high  wind  periods.  The  revised  weather  data  is  written  to  a STAT  file 
for  subsequent  use  as  input  to  the  DSPA  program. 


FIGURE  6-3.  PROFILE  COMPUTER  PROGRAM 
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PROGRAM  ALGORITHMS 


Step  1 


Read  monthly  user  input  request 


Step  2 


ALPHAQ  = 
ALPHAT  = 
ALPHAV  = 
ALPHHV  = 
ALPHLQ  = 


confidence  level  (0  to  l)  for  solar  insolation  data 

confidence  level  (0  to  1)  for  temperature  data 

confidence  level  (0  to  1)  for  wind  velocity  data 

confidence  level  (0  to  l)  for  high  wind  worst  case 

confidence  level  (0  to  l)  for  low  insolation  worst 


ALPHLV  = confidence  level  (0  to  1)  for  low  wind  worst  case 
LH(l)  = temperature  flag:  -1  ■ low  profile 

0 = means  profile 
+1  = high  profile 

LH(2)  = wind  velocity  flag  (-1,  0,  or  +l) 

LH(3)  = solar  insolation  flag  (-1,  0,  or  +l) 

LH(4)  = low  insolation  flag:  0 » no  worst  case  analysis 

+1  = worst  case  analysis 
LH(5)  * low  wind  flag  (0  or  +l) 

LH(6)  = high  wind  flag  (0  or  +l) 

PHV  = scale  factor  (>l)  for  high  wind  worst  case 
PLQ  * scale  factor  (0  to  l)  for  low  insolation  worst  case 
PLV  = scale  factor  (0  to  1)  for  low  wind  worst  case 
PQ  = proportion  (0  to  l)  for  solar  insolation  data 
PT  = proportion  (0  to  l)  for  temperature  data 
PV  = proportion  (0  to  l)  for  wind  velocity  data 

Compute  tolerance  limit  factors 
ZA  * F{ A, STABLE} 

ZP  = F{P,ZTABLE} 

AL  = 1.0  - ZA**2/(2.0#(NYRS  -1)) 

BL  « ZP«»2  - ZA**2/NYRS 

CL  » (ZP  + SQRT(ZP<*»2  - AL#BL ) ) / AL 

for  each  of  CLT,  CLV,  and  CLQ 
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where : A = specified  confidence  level 
P = specified  proportion 

ZTABLE  = table  of  inverse  error  function  vs.  percent 
CLT  = temperature  tolerance  limit  factor 
CLV  = wind  velocity  tolerance  limit  factor 
CLQ  = solar  insolation  tolerance  limit  factor 


Step  3 Compute  low/high  delta/soale  factors 

DLHT  = LH(l ) * CLT  • TM(2) 

RLHV  = 1.0  + LH(2)  * CLV  • VM(2)/VM(l) 

RLHQ  = 1.0  + LH(3)  * CLQ  * QM(2)/QM(l) 

where : DLHT  = temperature  delta  factor 

RLHV  = wind  velocity  scale  factor 
RLHQ  * solar  insolation  scale  factor 
TM(2)  * monthly  temperature  standard  deviation 
VM(l)  *=  average  monthly  wind  velocity 
VM(2)  = monthly  wind  velocity  standard  deviation 
QM(l)  * average  monthly  solar  insolation 
QM(2)  « monthly  solar  insolation  standard  deviation 


Step  4 


Read  next  STATS  file  day 

T(I)  * temperature  for  hour  I 
V(I)  * wind  velocity  for  hour  I 
Q(l)  ■ solar  insolation  for  hour  I 
for  I - 1,24 


Step  5 Modify  STATS  data 

T( I)  « DLHT  ♦ T ( I ) 

V( I)  ■ RLHV  * V(I) 

Q(I)  « RLHQ  * Q(I) 

Step  6 Write  modified  day's  data  to  STAT  file 
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Step  7 If:  last  day  of  month 
Then:  GO  TO  STEP  8 
Otherwise:  GO  TO  STEP  4 

Step  8 If:  worst  case  analysis  requested 
Then:  GO  TO  STEP  9 
Otherwise:  GO  TO  STEP  13 

Step  9 Compute  Poisson  maximum  likelihood  estimate 
LAMBDA  = (SUM  - CNT)/CNT 

for  each  of  high  wind,  low  wind,  and  low  insolation 

where : LAMBDA  = maximum  likelihood  estimate 

SUM  * total  number  of  bad*  days  for  current  month 
CNT  a number  of  strings  of  bad  days  for  current 
month 

Step  10  Compute  number  of  sequential  worst  case  days  for  current  month 


NR  such  that 


1 

J«1 


LAMBDA** J/J! 


> confidence  level 


NDYS 

£ 

j-i 


LMBDA**I/I! 


for  each  NRHV,  NRLQ,  NRLV 


where:  NDYS  ■ number  of  days  in  current  month 

NRHV  ■ number  of  sequential  high  wind  worst  case  days 
to  be  centered  about  the  20th  day 
NRLQ  « number  of  sequential  low  insolation  worst 

case  days  to  be  centered  about  the  15th  day 
NRLV  ■ number  of  sequential  low  wind  worst  case 
days  to  be  centered  about  the  10th  day 

•A  bad  day  is  a day  for  which  (Q  < PLQ*  Q „ ) 

— avg.  for  month 
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i 

Step  11  Modify  worst  case  day's  data  for  appropriate  days  of 
current  month. 

V(I)  = PLV  * V(l)  for  low  wind  period 

Q( I ) = PLQ  * Q(l)  for  low  insolation  period 

V(I)  PHV  * V(l)  for  high  wind  period 

Step  12  Re-write  worst  case  days  to  STAT  file 

Step  13  If:  last  day  of  year 

Then:  GO  TO  STEP  14 
Otherwise:  TO  TO  STEP  1 

Step  lU  Stop  Program 


